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NOTICES 


Readers are invited to submit written discussions on any of the papers contained in 
this issue. Discussions received within two months from the publication date of this 
issue will be published in a following issue. Correspondence should be addressed to 
Society Headquarters. 


ASLE is constantly looking for new and useful papers on every aspect of lubrication 
including research, design, application and maintenance. Authors having papers which 
they believe would be of interest to the readers of either ASLE Transactions or Lubri- 
cation Engineering are encouraged to submit them to the Society for publication con- 
sideration. If an author desires, he may request that a paper also be considered for 
presentation at one of the Society’s meetings. Presentation before a meeting, however, 
is not a prerequisite for publication. An Author’s Guide covering the preparation of 
papers for publication is available from Society Headquarters on request. 
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Contact Fatigue under Oscillatory Normal Load 


By J.C. TYLER,! R. A. BURTON,” and P. M. KU (aste)* 


Experimental studies were made wherein contact fatigue was brought about by a cyclic variation 
of normal load applied to the contact between a ball and a flat specimen. Curves of flat-specimen 
temperature versus time showed that a rapid temperature rise occurred in the initial stage of crack 
formation rather than after the appearance of fully developed spalls; thus the temperature criterion 
appeared to provide an early indication of fatigue. Fatigue failure was observed in unlubricated 
as well as lubricated contact, but with much longer life when unlubricated. 

Toroidal rings of hardened and softened material were noted both on and below the contact 
surface. These hardness variations bore symmetrical relationship to the contact geometry and 
were apparently related to the stress pattern. 

It was found that the fatigue cracks did not initiate in the most severely hardened material but 
in the material adjacent to it. Further, the slope of the cracks tended to follow the boundary of 
the hardened region in general direction. Propagation of microcracks in this general direction 
appeared to be influenced by the prior austenitic grain boundaries. These microcracks at times 


appeared to propagate along carbide stringers. 


Introduction 


In recent years, a great deal of interest has been shown 
by the jet engine industry, the bearing industry, and the 
lubricant industry in the problem of rolling-contact bear- 
ing fatigue. Because of its complex nature, particularly 
when full bearings are considered, this problem has not 
lent itself to a quick solution. Contributions to the under- 
standing of the failure process have been made at a 
number of levels: phenomenological studies of full bear- 
ings, simplified bearing element studies, elasto-hydro- 
dynamic lubrication studies, stress analysis along with 
the application of failure criteria, dislocation theories, 
and metallographic analysis of materials undergoing the 
characteristic stresses. Through such concerted studies, 
a basic understanding has been slowly emerging. 

The research reported here has been devoted to ob- 
taining information on several aspects of the problem by 
experimentally producing contact fatigue under delib- 
erately simple circumstances. No attempt was made to 
simulate actual bearing operation. However, it is believed 
that the test conditions were sufficiently realistic to 
permit insights gained to be applied to rolling-contact 
bearing fatigue. To accomplish this, an apparatus was 
developed to produce contact fatigue at a single point, 


Contributed by the ASLE Technical Committee on Bearings and 
Bearing Lubrication and presented at the Annual Meeting of the 
American Society of Lubrication Engineers held in New York, 
April, 1963. 

1 Associate Research Engineer, Southwest Research Institute, 
San Antonio, Texas. 

2 Staff Scientist, Southwest Research Institute, San Antonio, 
Texas. 

3 Director, Department of Aerospace Propulsion Research, 
Southwest Research Institute, San Antonio, Texas. 


through the application of an oscillatory normal load 
between a ball and a flat surface. This technique was 
also found to lend itself to early detection of fatigue 
failure through observation of changes in the flat speci- 
men temperature during testing. 

It should be noted that Kennedy (7) had earlier used 
oscillatory normal loading in contact fatigue studies. The 
present work is felt to represent a considerable extension 
of that initial effort. In particular, the present work has 
been facilitated by the thermal damage-detection tech- 
nique. This has proven much more sensitive to the early 
stages of damage than the ultrasonic method reported by 
Kennedy. 


Apparatus 


The principal features of the oscillatory fatigue ap- 
paratus are shown in Fig. 1. The test specimens con- 
sisted of a ball (attached to a tool holder by an epoxy 
resin) in contact with a flat, disc-shaped specimen at D. 
The flat specimen was held by a resonant support, B, 
and rested on a piezoelectric crystal. The specimens were 
statically loaded by means of pneumatic pressure applied 
at C and acting on a bellows from below. This load was 
transmitted through fastenings at nodes of the resonant 
support, so as not to affect the vibration of the rig. In 
addition, a dynamic load was also applied to the speci- 
mens by a magnetostrictive transducer G, which caused 
the tool holder A to vibrate in the longitudinal direction 
at the driver frequency, 26 kc. The entire rig was O-ring 
sealed to permit control of the atmosphere inside the 
apparatus, 

Figure 2 illustrates the load variation on the specimens 
during operation. The instantaneous load varied between 
a minimum, a, and a maximum, b. The contact area 
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Fic. 1. Assembly of fatigue apparatus 


varied with load, with diameter d, and d, corresponding 
respectively to loads a and b. 

Figure 3 shows the mounting arrangement of the 
barium titanate crystal. The ball specimen held in the 
tool holder was forced against the flat-faced specimen 
which set on the suitably insulated, wafer-shaped crystal. 
Insulation between the lower specimen and the crystal 
was provided by a thin layer of plastic. Contrary to 
intuitive expectations, the compliance of the plastic was 
not high since it could deform appreciably only in com- 
pression under the rapid stressing. The accelerometer 
effect, due to the vibrating mass of the crystal, was in- 
vestigated and found to be insignificant for this partic- 
ular application. It is also of interest that the wave form 
of the output signal from the crystal was not perfectly 
sinusoidal, as shown in Fig. 2, but was asymmetrical, 
since Hertzian contact acted as a nonlinear spring. 

The output of the BaTiO, crystal was calibrated by 
subjecting the specimens and support system to selected 


STATIC 
LOAD 


a 
n 


Fic. 2. Variation of load and contact area on fatigue speci- 
mens. 
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static step-loads and reading the piezoelectric output 
with an electrometer. This was further corrected for the 
alternating current losses due to the “loading” effect of 
the panel instruments during oscillatory operation. 

A ;,-inch stainless steel bayonet thermocouple was 
inserted into the interior of the flat specimen as shown 
in Fig. 3. The thermoelectric potential from the thermo- 
couple was recorded, thus giving a permanent flat-speci- 
men temperature history for each test run. 

Because the area of interest on the specimens was rela- 
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Fic. 3. Crystal load cell and thermocouple arrangement 
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Fic. 4. Surface profile of flat specimen after one load cycle without lubricant 


tively small, it could be studied intensively by many 
other physical techniques. These included sectioning, 
etching, microhardness measurement, indentation profile 
measurement, and electron beam microanalysis. 


Description of specimen damage 


When a ball and a flat specimen, in contact, are sub- 
jected to oscillatory normal loading, there are two dis- 
tinct contact circles formed on both of the specimens 
(see Fig. 2). The smaller or inner circle corresponds to 
the minimum load since the specimens do not separate 
during testing. The larger or outer circle corresponds to 
the maximum load. The area between these two circles 
represents that portion of the surface which is alternately 
brought into contact with the mating surface and pulled 
away from it. 

As a consequence of the oscillatory loading, three dis- 
tinguishable types of specimen damage were observable 
following testing, these being fretting, spalling, and 
plastic deformation. 

The fretting exhibited itself as a light erosion and 
staining of the region between the inner and outer con- 
tact circles. It appeared very early after cyclic operation 
was started but remained secondary to the other types 
of damage as to the depth of penetration, even for ex- 
tended test runs of several million cycles duration. It 
could not be eliminated by the use of oil in the contact 
zone and by the presence of an inert gas atmosphere. It 
was worsened when the elastic modulus of the ball dif- 
fered greatly from that of the flat, such as when a ceramic 
ball was used against a steel flat. 

Spalling was clearly identifiable from fretting both by 
the depth affected and by its appearance under the micro- 
scope. Though the term “spalling” is used for con- 


venience of discussion, this type of damage progressed 
as a succession of events: crack formation, crack growth 
and interweaving, and flaking out of the crack-surrounded 
areas. The cracks and the circumferential groove result- 
ing from a succession of spalls tended to start just outside 
the periphery of the inner contact circle and ultimately 
to extend both inward and outward as the damage ad- 
vanced. 

Plastic deformation or depression of the center of the 
inner contact circle could be clearly distinguished from 
the other types of damage since the surface within the 
inner circle remained undamaged during testing. Even 
after very long periods of operation, the mirrorlike finish 
of this surface was preserved, and light polishing scratches 
remained on it undisturbed. Thus, a measure of plastic 
deformation of the contact zone was provided by the 
depth to which the center of the inner circle was de- 
pressed below the level of the surrounding surface. 

Figure 4 shows the surface profile of a 52100 steel flat 
specimen, subjected to one normal load application 
against a 52100 steel ball at a calculated maximum Hertz 
stress of 750,000 psi, without any lubricant present. The 
maximum plastic deformation, occurring at the center of 
the contact region, was approximately 15 pinches. In an- 
other test under identical conditions, where the speci- 
mens experienced 750,000,000 load cycles without fatigue 
failure, the surface profile of the flat specimen was simi- 
lar but the maximum depth of the depression was ap- 
proximately 40 pinches. In this case of dry contact, the 
maximum amount of fretting was of the order of 10 
winches after 750,000,000 load cycles. 

It is interesting to compare the accumulation of plastic 
deformation with repeated normal loading, as observed 
here, with that for rolling contact. Moyar and Sinclair 
(2) presented data for a flat plate grooved by a rolling 


257 

peci- 
the 
t of 
own 
‘mo- 
eci- 
rela- 
ED 


J. C. Tyzer, R. A. Burton, P. M. Ku 


CALCULATED OUTER 


CONTACT DIAMETER 


CALCULATED 
INNER CONTAC 
DIAMETER 


lOO MICRO IN. 


| 


Fic. 5. Surface profile of fatigued flat specimen after 162,000,000 load cycles with lubricant 


ball at a maximum contact stress of 730,000 psi, where 
the depth of the groove in the flat plate increased from 
about 53 winches at 6000 stress cycles to 84 winches at 
1,340,000 stress cycles. These observations are at vari- 
ance with Tabor’s findings (3) which indicated no ac- 
cumulation of plastic deformation with repeated loading. 

Turning now to a fatigue damaged specimen, Fig. 5 
shows the surface profile of a flat specimen subjected to 
162,000,000 load cycles under the same test conditions 
with an SAE 10 white mineral oil in the contact zone. 
There was a deep circumferential groove between the 
inner and outer contact circle, and a fatigue crack ex- 
tended circumferentially for approximately 100 degrees 
along the bottom of the groove. The true depth of the 
crack is not shown on the profile because the acute angle of 
the crack prohibited full penetration of the traversing 
stylus. The figure further shows that the maximum depth 
of the depression at the center of the contact zone was ap- 
proximately 80 winches. This was considerably more than 
that for the case where no fatigue failure occurred, 
possibly because the presence of the circumferential 
groove reduced the actual contact area and therefore 
increased the stress level in the central area. 

Figure 6 shows a typical flat-specimen temperature 
trace versus time. Upon the application of oscillatory 
load (Regime I), the temperature rose steadily at first, 
then gradually levelled off, and thereafter remained con- 
stant so long as no fatigue was observed. The abrupt rise 
in temperature (Regime II) was found to be related to 
the appearance of fatigue cracks. It has been shown (4) 
that the energy dissipation prior to crack formation (i.e., 
that associated with contact hysteresis and wear) pro- 
ceeded in an orderly and essentially predictable manner; 
however, the magnitude of the initial temperature rise 


bore no apparent relationship to the inception of fatigue. 
On the other hand, whenever fatigue cracks were devel- 
oped, a “break” in the temperature rise curve invariably 
resulted, and this was taken as evidence of increased 
energy dissipation associated with the cracks. The tem- 
perature break was found to be a very sensitive means of 
fatigue detection. In tests stopped immediately after the 
temperature break, very small cracks were at times found 
in post-test inspections. Advanced spalling and grooving 
were found only when the test was allowed to continue 
long after the temperature break. 
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Fic. 6. Variation of flat-specimen internal temperature 
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Metallurgical observations 


Many of the test specimens were subjected to various 
metallurgical investigations such as cross sectioning, etch- 
ing, and extensive inspection under the microscope in 
order to gain information associated with specimen 
damage. Standard techniques were employed for prepar- 


Fic. 7. Details of fatigue crack in flat 52100 steel specimen. 
(A) crack relative to austenitic grain structure; (B) cross section 
of crack. 


ing metallographic specimens, but extreme care was taken 
during grinding to avoid overheating. 

Figure 7A illustrates a typical fatigue crack in one of 
the flat 52100 steel specimens run with an SAE 10 white 
mineral oil in the contact area. The principal crack had 
small cracks propagating from it along curved paths as 
shown by the arrows | and 2. If the testing had continued 
on this specimen, it is very likely that the small cracks 
would have extended until they again joined the principal 
crack, then the material completely surrounded by cracks 
would have been removed from the specimen giving rise 
to a spall or pit. An area where this had already hap- 


pened is shown by arrow 3. As the principal crack con- 
tinued to propagate and the smaller cracks continued to 
grow from it in this manner, a circumferential groove 
would form just outside the inner contact circle and 
then spread both inward and outward, as explained previ- 
ously. Figure 7B shows the crack as it appeared on the 
sectioned surface of the test specimen. The crack extended 
under the contact area at an angle of approximately 30° 
and reached a maximum depth of approximately 0.0015 


Fic. 8. Details of fatigue crack in M-10 steel ball. (A) crack 
in surface; (B) cross section of crack. 


inch below the contact surface, It is interesting to note, 
though not surprising, that the surface material was 
much softer (refer to the Knoop hardness indentations 
in Fig. 7A) in the area directly over the extension of the 
crack into the subsurface material. 


Figure 8 illustrates a typical fatigue crack in an M-10 
steel ball run with an SAE 10 white mineral oil in the 
contact area. The crack on the surface (A) of the ball 
did not have the small circular cracks propagating from 
it, as in the previous case. But there were definite exten- 
sions from this principal crack, and evidence relating 
these to the prior austenitic grain boundaries will be 


- 
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presented later. Again the crack, as shown on the sec- 
tioned surface (B), extended under the contact area and 
had a depth of approximately 0.0012 inch below the 
contact surface. As before, the crack formed an angle of 
approximately 30° with the contact surface, sloping 
downward toward the contact center. 


Fic. 9. Details of failed area in flat 52100 steel specimen run 
with dry contact. (A) surface of flat specimen after sectioning ; 
(B) cross section of specimen through failed region. 


Figure 9 illustrates both the contact surface (A) and 
sectioned surface (B) of a failed flat 52100 steel speci- 
men that was run with a dry contact. The pit was rela- 
tively large but only extended to a depth of approxi- 
mately 0.0008 inch. There was a very small crack in the 
material diametrically across from the pitted area. This 
crack was visible in the contact surface but was so small 
that it did not appear in the sectioned view. Location of 
the failure relative to the contact area was in the same 
zone as those observed on specimens tested with a lubri- 
cant in the contact area. After a light nital etch, the 
sectioned surface appeared as shown in Fig. 10. The dark 
streaks in the material below the pitted area have been 
identified as small cracks, and a definite tendency can 
be observed for these cracks to follow a path of carbides 
in the material. There has been other evidence indicating 
that the fatigue cracks tended to follow carbide stringers 
in both the M-10 and 52100 steel specimens. Unfortu- 


nately, these details have not been clearly shown in the 
photomicrographs. This is because the carbide stringer 
is most apparent at low magnification where the carbide 
particles themselves cannot be identified. At higher mag- 
nification, the stringer loses identity and one sees only 
individual particles. This offers no special problem for 


Fic. 10. Appearance of subsurface cracks following nital etch 
of specimen shown in Fig. 9. 


visual observations when magnification may be changed 
easily, but prohibits clear representation of the situation 
in photographs. 

Many of the test specimens were etched in a solution* 
to bring out the prior austenitic grain structure. The 
purpose of the etch was to determine if any abnormalities 
observed in the specimen material were related to the 
austenitic grain structure. Photographs that were typical 
of these observations are shown in Figs. 7A, 11A, and 12. 
There is a strong indication that the initial fatigue crack 
did not follow the austenitic grain boundaries. The fact 
that the cracks did coincide with the grain boundaries 
in some areas appeared to be incidental rather than the 
rule. On the other hand, outgrowths (shown by the 
arrows) from the pitted or spalled area in Figs. 11A and 
12 tended to propagate along the austenitic grain bound- 
aries. It is conceivable that as these outgrowths increase 
in number and extend farther along the grain boundaries, 
they might either consume or undermine the essentially 
undamaged material between them, and by this process 
the pitted or spalled areas would continue to grow. An- 
other interesting observation displayed by Fig. 12 was 
the enormous amount of slip, as shown by the numerous 
slip lines, in the specimen material immediately sur- 
rounding the fatigued area. 

Figure 13 illustrates the structural alteration observed 
beneath the contact area in a failed M-10 steel ball. 
There was a crack in the surface of the ball after testing. 
Material was removed from the surface of the ball until 
only the deepest portion of the crack remained, then the 
ball was sectioned through the cracked area and sub- 
jected to a nital etch. The material directly below the 
contact zone was etched more rapidly, thus creating the 
darkened region shown in the photograph (Fig. 13A). A 


4 For 52100 steel, chromic acid and caustic soda in water was 
used; for M-10 steel, 50% HCl was used. 
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Fic. 11. Details of pits in flat 52100 steel specimen. (A) pits 
relative to austenitic grain structure; (B) cross section cf pits. 


Fic. 12. Details of fatigue failure in M-10 steel ball relative 
to austenitic grain structure. 


Fic. 13. Details of transformed region in M-10 steel ball 
(A) cross section of specimen through failed region; (B) enlarged 
view of split carbide. 


split carbide particle (arrows) was found in the specimen 
material directly below the region occupied by the previ- 
ously mentioned crack. This carbide particle was at the 
edge of the transformed zone and actually laid in a 
region where the darkened material extended into the 
material that had not been transformed. No carbides in 
this or other specimens were observed with tongues or 
intrusions. 

One flat specimen that had exhibited a “break” in the 
specimen temperature curve during testing, but did not 
appear to have any surface cracks when cleaned by super- 
ficial polishing of the contact surface, was etched to bring 
out the austenitic grain structure. Microscopic observa- 
tion of the contact area, after the etch, revealed a very 
small crack in the surface of the specimen. Subsequent 
removal of additional material from the surface of the 
specimen provided evidence that the crack was increasing 
in length and width as it extended into the subsurface 
layers. After approximately 0.0001 inch of material had 
been removed in this manner, the specimen was sectioned 
to determine the depth and appearance of the crack in 
the subsurface region. Figure 14 shows that the maximum 
depth of the crack was approximately 0.0013 inch and, 
as before, the angle between contact surface and crack 
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Fic. 14. Crack bifurcation in subsurface material of flat 52100 
steel specimen. 


was approximately 30°. There was a bifurcation of the 
crack at a depth of approximately 0.001 inch. Metallo- 
graphic examination did not reveal any definite non- 
homogeneity associated with this bifurcation. 


Electron beam microanalysis 


Two metallographically prepared samples were ana- 
lyzed with an electron beam microanalyzer. One sample 
was a failed 52100 steel flat specimen and the other a 
failed M-10 steel ball. The specimens were analyzed for 
second phases, austenitic grain boundary segregation, and 
segregation in the area of the cracks. Paths were traversed 
on the specimens, and, in addition, separate discrete 
points in the material were analyzed. 


52100 STEEL SPECIMEN 


Of interest was segregation at the austenitic grain 
boundaries and segregation in the area of the crack. The 
grain boundaries were found to be enriched in chromium 
and manganese, but the silicon distribution was relatively 
uniform. One point in the crack was noted to be slightly 
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RADIAL DISTANCE FROM CENTER OF CONTACT, IN. 


rich in silicon. A visible fluorescence was noted when 
the electron beam was positioned on this point implying 
the presence of a discrete particle of Si0.. There were 
indications, though not conclusive, that a chromium en- 
richment surrounded the crack, while in the crack itself, 
the concentration was normal. 


M-10 STEEL SPECIMEN 


Of interest was the carbide chemistry and segregation 
at the grain boundaries and in the matrix surrounding the 
carbides. The carbides were found to be rich in vanadium 
and molybdenum with small amounts of chromium. Al- 
though no analyses were performed for iron, it seemed 
unlikely, based on the determined V-Mo-Cr values, that 
more than a few per cent of iron was present in the 
carbides. The matrix concentration levels for areas within 
5-10 u distance of the carbides agreed with the normal 
matrix levels within the experimental error. The grain 
boundaries were noted to be enriched in vanadium and 
molybdenum. 


Microhardness variations 


Numerous microhardness traverses were made across 
the contact areas of the flat specimens. Through a process 
of removal of material layers and repeated hardness 
measurements, both surface and subsurface regions were 
examined. Failed and unfailed specimens were found to 
exhibit the same pattern of hardness variation, but this 
variation was more severe in the failed specimens. 

Figure 15 illustrates details of the hardness variation 
on a typical failed 52100 steel flat specimen. The speci- 
men had an average Knoop hardness of 780 before test- 
ing. All of the hardness indentations were made with a 
500-gm load. There was an extremely hard zone on the 
specimen surface immediately inside the inner contact 
circle. The hardness decreased with depth into the 
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Fic. 15. Knoop hardness contours in cross section of flat specimen 
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specimen. The center of the contact area was much softer 
than this annular ring, even on the surface of the speci- 
men, and the hardness of this center zone decreased with 
depth until it reached a minimum at a depth of approxi- 
mately 0.002 inch. There were also other clearly identi- 
fied soft zones in the specimen subsurface material. Due 
to the nearness of symmetry of these soft zones, they 
appeared to form toroidal rings in the material. It is 
highly unlikely that these soft zones were in the material 
prior to testing because they were symmetrical relative 
to the contact area, because numerous tests outside of 
the contact areas have not given such a variation in 
hardness, and because additional but less detailed studies 
of other specimens have confirmed many of the major 
features of the distribution. 

Regions of hardness variations were noted by Jones 
(5), Carter (6), and Bush e¢ al. (7) in rolling contact. 
There are some differences in the observed hardness and 
softness distributions with those reported here. These 
differences can be explained by the fact that rolling must, 
of necessity, smear the zones and eliminate the possibility 
of a clear resolution of such regions. 

To determine if the soft spots were the result of 
thermal annealing, estimates were made of contact tem- 
perature using a simple experimental setup where tem- 
perature drops between the specimens were measured 
for different rates of heat flow from one to the other (4). 
These indicated that the time-mean contact temperature 
was only about 175 F during operation without failure. 
In view of this, it was felt that cyclic variations could 
scarcely be sufficient to bring this periodically up to the 
temperature range necessary for thermal annealing of the 
contact zone. Of course, if a crack forms and dissipation 
rises locally, it is possible that annealing can occur and 
influence propagation. 

Another item of interest is the fact that fatigue cracks 
in the surface material were always in the zone immedi- 
ately outside of the inner contact circle and were circum- 
ferential in nature. This is an indication that the fatigue 
cracks did not initiate in the most severely hardened 
material but in the material adjacent to it. Further, the 
slope of the cracks inward tended to follow the boundary 
of the hardened region in general direction. 


Stress analysis 


When the hardness data and observations on crack 
formation were checked against predictions of fatigue 
failure from existing theories (8), contradictions were 
found. Whereas such theories would call for deeper 
sources of failure initiation, experimental results indicated 
that the region of interest might be near the surface of 
contact. In the hope of resolving this question, calcula- 
tions of surface shear stress were undertaken. While the 
results have been found to be provocative, they still do 
not provide an answer as to what aspect of stress varia- 
tion actually causes crack formation. 


Fic. 16. Calculated surface shear-stress variations 
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Surface shear stresses were calculated from the principal 
stresses of classic Hertzian theory (9), for a ;,-inch 
steel ball in contact with a steel flat. The results are 
shown in Fig. 16. The load range which served as a basis 
for these calculations produced a variation of maximum 
Hertz stress from 398,000 to 695,000 psi. These Hertz 
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Fic. 17. Variation of calculated octahedral shear stress on 
specimen surface. 


stresses give a contact radius of 0.006 and 0.0104 inch, 
respectively. The shear stresses, Tz, and were 
calculated for both the minimum load and the maximum 
load, thus giving representative changes of the maximum 
shearing stresses for oscillatory normal loading. Values 
for the shear stresses t,g and t,:, shown in Figs. 16A 
and 16B, change in magnitude only during a stress 
cycle, whereas the values for t.4, shown in Fig. 16C, 
change in both magnitude and direction. Due to this 
reversal of direction, the change in t,4 at the minimum 
contact radius (0.006 inch) is greater than the change 
in Tt, OF Ty: at any contact radius during a stress cycle. 
This change in t.4 of 76,000 psi is, however, only slightly 
greater than the change in 1, of 74,500 psi at a radius of 
0.0104 inch. The fact that the direction of t., does 
change may, of course, be of significance. An observation 
that seems to verify this significance is that fatigue 
cracks in the surface material have always appeared 
in the same general vicinity of the specimen as this 
maximum change in shear stress. 

The octahedral shear stresses of the surface material 
at different contact radii were calculated in terms of 
the principal stresses. These calculations were also made 
for a variation of maximum Hertz stress from 398,000 
to 695,000 psi. The variation of the octahedral shear 
stress across the specimen contact area, for four different 
Hertz stresses, is shown in Fig. 17. Since the octahedral 
shear stress is very closely related to the distortion 
energy theory of failure, it may be logical to assume 
that the change in octahedral shear stress in the speci- 
men failure zone would be a maximum, but the curves 
in Fig. 17 contradict this assumption; in fact, the change 
in octahedral shear stress is maximal at the outer rather 
than inner contact radius. 

A more informative presentation of surface stress 
behavior is shown in Fig. 18. Here octahedral stress is 


given as a function of time for selected. radii from the 
center of contact. Note that the stress variation just out- 
side the inner contact (r = 0.007 inch) passes through 
three peaks per cycle while the much higher stress at 
r = 0.0104 inch passes through only one peak per load 
cycle. This observation in itself is not sufficient to 


80 x 107 
RADIUS FROM CENTER OF 
CONTACT, IN. xX 103 
10.4-—— 
f 
60+ / 4 
|/ WA \ 
30 J 
4 
ONE STRESS CYCLE——————+} 
° 


TIME 


Fic. 18. Time variation of octahedral shear stress on specimen 
surface. 


explain the observed location of fatigue cracks. On the 
other hand, it does suggest that factors are at work which, 
until now, appear to have been neglected. Totally aside 
from considerations of maximum shear reversal, questions 
regarding the number of stress peaks per load cycle are 
now raised and must ultimately be dealt with. In addition, 
better criteria of failure are needed in the presence of 
triaxial stress variations when the hydrostatic stress 
component is significantly large. 


Fatigue life observations 


Defining “fatigue life” as the number of stress cycles 
before the flat-specimen temperature break, Weibull 
plots have been prepared as shown in Figs. 19A and 
19B. Figure 19A shows data for 52100 steel balls in 
contact with 52100 steel flats, and is based on eight 
failures in the twelve tests performed. The maximum 
Hertz stress averaged 734,000 psi and SAE 10 white 
mineral oil was applied as lubricant. 

Figure 19B shows a corresponding Weibull plot 
for M-10 steel balls operating on 52100 steel flat speci- 
mens, with five failures in a total of eight tests. The 
lubricant was again SAE 10 white mineral oil and the 
maximum Hertz stress averaged 730,000 psi. 
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In analyzing the data of the twenty tests in the 
above two Weibull plots, a definite correlation was 
found between the probability of specimen failure and 
applied deformation amplitude as determined by the 
power setting of the driving generator. Since an initial 
constant dynamic load (crystal output) had served as 
a criterion for setting the power control of the generator, 
variation of this setting was an indication of variations 
of local elastic compliance of the specimens. 
Proceeding under the assumption that power setting 
variation had contributed to the scatter of the first test 
series, it was decided to run a series of tests wherein the 
power to the transducer would be held constant, irrespec- 
tive of indicated dynamic load. Twelve such tests were per- 
formed and a Weibull plot of the ten failures obtained 
is shown in Fig. 19C. The lubricant was SAE 10 white 
mineral oil and the maximum Hertz stress was 750,000 
psi. Since the plot does not appear to be a continuous 
function, two straight lines (dashed) may be drawn. The 
best single line (solid) to fit the data points (calculated 
by the method of least squares) appears to have a slope 
almost identical to that shown in Fig. 19A. If the best 
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Fic. 19. Weibull plots of experimental data 
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fitted line from this test series is corrected for the stress 
level by using the inverse tenth-power relation of fatigue 
life with maximum Hertz stress as proposed by Carter 
(6), it is found to be very nearly the same as that in 
Fig. 19A in both slope and position. This indicates that 
the peculiar arrangement of ‘the data in Fig. 19C may 
simply be scatter and that a sufficient amount of the 
data has not been collected to straighten the line out. 
If this be so, constant dynamic load tests and constant 
amplitude tests give closely comparable results. 
Figure 19D is a Weibull plot for 52100 steel balls 
in contact with 52100 steel flat specimens with a 
dry contact. There were five failures in eleven test 
runs. The maximum Hertz stress was 750,000 psi. 
The solid line was again drawn from calculations 
using the method of least squares. The Bio life 
for this series of dry tests was approximately 40 
times as great as the value obtained with mineral oil 
in the contact area (Fig. 19C). This large increase in 
fatigue life may very well be significant in explaining 
why Way (10) concluded that pitting due to rolling 
contact would not occur without a lubricant. The maxi- 
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mum duration of any dry tests performed by Way was 
18,000,000 stress cycles, where 50% of his lubricated 
rollers pitted between 670,000 and 1,450,000 stress 
cycles. Confirmation of longer life in the absence of 
lubricant must not be taken as confirmation of Way’s 
theory of mechanical participation of the lubricant in 
crack growth, however. 


Conclusions 


The oscillatory normal contact between a ball and a 
flat, as employed in this investigation, was not intended 
as a direct simulation of rolling contact. However, it 
has been found to be uniquely suited for studying some 
aspects of the rolling-contact fatigue problem. Cracking 
and spalling of the specimens could be produced in the 
confined regions subjected to high stress. Furthermore, 
the concentration of damage permitted intensive analysis 
of the contact regions. 

Curves of flat specimen temperature versus time showed 
that an abrupt temperature rise occurred in the initial 
stage of crack formation rather than after the appearance 
of fully developed spalls; thus the temperature criterion 
appeared to provide an early indication of fatigue. 
Weibull plots of experimental data for both lubricated 
and unlubricated oscillatory normal contact based on 
the temperature criterion exhibited a similar trend 
to those for rolling contact based on the incipience of 
spalls. The fatigue life for unlubricated normal contact 
was found to be much longer than that for lubricated 
contact. 

Measurements on the permanent deformation of the 
flat test specimen for both single and multiple stress 
cycles were made. These measurements led to the con- 
clusion that there was an accumulation of plastic defor- 
mation with successive load application. 

Toroidal rings of hardened and softened material 
were observed both on and below the contact surface. 
These hardness variations bore symmetrical relationship 
to the contact geometry and were apparently related 
to the stress pattern. It was found that the fatigue 
cracks did not initiate in the most severely hardened 
material, but in the material adjacent to it. The slope 
of the cracks tended to follow the boundary of the 
hardened region in general direction. 

The occurrence of hardened and softened rings in 
clearly defined zones argued against a general annealing 
process as a result of energy dissipation in the contact. 
Furthermore, thermal measurements have established 
the general magnitude of temperature rise to be con- 
siderably below that necessary for annealing at an 
appreciable rate. In view of these, it was concluded 
that the softened and hardened regions were the result 
of peculiarities in the response of the material to the 
three-dimensional aspects of the stress variation. If this 


be the case, refinement of existing criteria for work 
hardening and failure appears necessary before an ex- 
planation of these zones can be found. 


Possible participation of the prior austenitic grain 
boundaries in fatigue failures was of interest. These 
boundaries, as were the carbides, were shown to be en- 
riched in vanadium and molybdenum in the M-10 steel. 
Photomicrographs also displayed an abundance of car- 
bides along the austenitic grain boundaries in the M-10 
steel; for the same steel, cracks were randomly observed 
to follow carbide stringers. A chromium and manganese 
enrichment at the austenitic grain boundaries in 52100 
steel specimens was observed. There was evidence, 
although not conclusive, that the matrix surrounding a 
fatigue crack in this material was also enriched in 
chromium. The propagation of previously formed pits 
or spalls was clearly related to these boundaries. These 
observations taken together point toward the conclusion 
that, in the steels used, the prior austenitic grain 
boundaries participated in fatigue failures. 
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Trpor TALLIAN (SKF Industries, Inc., Philadelphia, Pennsylvania) : 


The authors have presented very substantial work in the rela- 
tively poorly explored field of Hertzian point loading in cyclic 
compression. Aside from Kennedy (Ref. / in the paper), only one 
other paper comes to mind in this field (A1). The value of this 
technique is, of course, the ability to concentrate all test occur- 
rences in a very small volume of material. For the metallurgist 
desiring to examine structural changes in fatigue, no valid objec- 
tion seems to arise to this experimental arrangement. 

There are, however, several objections from a mechanics point 
of view to using results obtained by this technique in interpreting 
rolling contact fatigue. The authors are aware of these, but it may 
be helpful to summarize them as follows: 

(a) In the cyclic compression test, no elastohydrodynamic film 
exists between the contacting surfaces. None would be expected 
from hydrodynamics and none is shown by the results which 
show unavoidable fretting and wear. As a consequence, high 
values of surface traction must exist in the contact area (A2) 
since the contact is between bodies of different curvature. The 
authors observe, accordingly, cracks initiating at, or near, the 
surface, in the peripheral zone of the contact, where slip is high. 
They do not observe fatigue in the layer of maximum shear 
stress, which would be much deeper below the surface (A3). 

Contrary to this situation, rolling contacts are partially or fully 
separated by elastohydrodynamic films under most fatigue test 
conditions (A4), and there is, accordingly, no fretting, little or no 
progressive wear, and fatigue cracks initiate classically in the zone 
of maximum shear stresses. Conditions when cracks initiate at 
the surface do exist, but even then they penetrate generally to the 
maximum shear zone. 

(b) The critical subsurface shear stress in cyclic compression is 
the maximum shear stress in the axis of symmetry, whereas it is 
the maximum alternating shear stress in the case of rolling con- 
tact. These two stresses operate on different planes which again 
may lead to differences in crack formation. 

Once these differences are recognized, the authors’ results become 
highly interesting because they may well be symptomatic of the 
surface-layer phenomena occurring in inadequately lubricated 
rolling contacts. Remarkable in this respect are the hardness 
changes reported. Their detailed explanation in terms of strain 
history seems a worthwhile project for the future. 

The results of the electron-beam microanalysis are probably a 
disappointment to the authors. It is surmised that they expected 
more salient features around the cracks. It is suggested that the 
lack of such features may be due to the fact that the type of 
stressing used causes cracks to originate at surface defects rather 
than defective structures in the material. 

The authors give surface shear stress variations calculated from 
“classic Hertzian theory.” The analysis is not included, but it 
appears that it was based on a frictionless contact. If friction 
(surface traction) is considered, shear stresses of much higher 
order of magnitude are obtained (A5, A6). (On the scale of the 
surface asperities worn away, these stresses obviously exceeded 
the shear strength of the metal.) It is believed that the actual 
shear stresses existing in the surface are quite sufficient to explain 
the surface failures observed. Again, it is noted that such high 
shear stresses cannot be transmitted by the elastohydrodynamic 
films often existing in rolling contact which is why these contacts 
do not fail at the surface. 

By far the most interesting observation of the authors is the 
existence of a discontinuity in the heat generation rate during the 
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tests. It is stated that this occurs when a “failure” takes place. 
The existence of a correlation between these two phenomena 
would suggest that the extra heat generated after failure stems 
from the vicinity of the “failure” or crack. The observed cracks 
are quite small and one is led to wonder about the high rate of 
heat dissipation per unit volume that this concept implies. If it 
could be shown by further study that a fatigue crack does, indeed, 
generate heat at a high rate once it is formed, this would be a 
valuable basic finding. 
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Joun S. Tawresey (The Franklin Institute Laboratories for Re- 
search and Development, Philadelphia, Pennsylvania) : 


Those of us associated with the progressive development of 
ball and roller bearings have been compelled to live with numerous 
vagaries influencing their fatigue life performance. Accordingly 
we welcome any well-considered effort directed to study of the 
fatigue process. The authors have adopted a new approach which 
has the advantage of simplification, in dealing with a single com- 
ponent element, rather than the complex system comprising a full 
scale bearing. However, it is understood that in their investigation 
no attempt has been made to actually simulate the conditions that 
apply in a rolling contact bearing. 

During the past year the writer’s laboratories have been en- 
gaged in an extensive test of the life performance of needle roller 
bearings of the type employed on helicopter rotor horizontal and 
vertical hinges. Such bearings do not rotate, but oscillate back 
and forth about a mean position, in cyclic relation to the aero- 
dynamic loads imposed upon the blade. The load contact between 
roller and raceways maintain a fixed indexing and extends over a 
localized area of about 5°. Test runs are made under controlled 
conditions of load, speed, angle of oscillation, and type of lubri- 
cant. The tests are continued until fatigue failure is made evident 
by spalling of the load surfaces. Some 400 bearings have been 
life tested in this program under various operating conditions. 

This extensive test of full scale bearings, with limited and 
localized load areas, subjected to repetitive contact stress cycles, 
would appear to afford an opportunity for study in connection 
with the author’s work and conclusions. It would indeed be 
gratifying to establish a basis of correlation between full scale 
bearing tests and a streamlined fatigue investigation such as that 
developed by the authors. This would appear to be a first and 
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logical step in generalizing the application of the method to the 
problems of rolling contact. 

Particular interest is taken in Fig. 19 of the paper, wherein 
Weibull plots are shown for the experimental data. Plate “C”, for 
52100 steel and Hertz stress of 750,000 psi, is in striking agree- 
ment with similar data obtained in the bearing tests referred to. 
The break in the Weibull distribution line, at about 50% statisti- 
cal percent failed, has been found to be characteristic of many 
tests with oscillating bearings. 

We believe that the author’s method merits further study as a 
means for evaluating fatigue performance, in the presence of high 
surface stresses, and possibly in the case of rolling contact. 


K. L. Jounson (Engineering Laboratory, Cambridge, England): 


1. Stress cycles in rolling contact. There is now much evidence 
that the loads at which contact fatigue becomes a serious prob- 
lem are sufficiently high to strain the material locally above its 
elastic limit. The present paper supports this view. Now in a 
loading cycle in which the material is deformed plastically, the 
stress and deformation are dependent upon the history of loading. 
The stress cycles experienced by material elements due to an 
oscillatory normal load are likely to be very different, therefore, 
from the stress cycles in rolling contact. (See J. E. Merwin and 
K. L. Jounson, “An analysis of plastic deformation in rolling 
contact,” Proc. Inst. Mech. Eng. London, 1963.) 

Keeping this point in mind, however, the investigation has given 
rise to a number of interesting questions concerning repeated 
contact loading. 

2. Plastic deformation during unloading. The paper reports a 
progressive increase in plastic deformation with cycles of load 
and shows bands of softened material beneath the surface which 
suggests a repeated cycle of reversing plastic strain (see the work 
of L. E. Coffin). These observations are contrasted with the work 
of Tabor (1948) which describes the unloading and reloading 
processes as perfectly elastic. 

I have given some thought to this question and have come to 
the conclusion that a heavily loaded contact will undergo reversed 
plastic strain on unloading and thus experience a continuous cycle 
of plastic strain with repeated loading. This may be shown as 
follows, following Tabor “Hardness of Metals” (Oxford, 1951). 

During first loading the surface deforms to give an indentation 
with a contact radius a. The exact pressure distribution is un- 
known but is not likely to differ greatly from the Hertz form of 
distribution! with a maximum pressure pf, given by: 


Po = % Pay. = [C1] 


where Y is the yield stress in tension and c is a factor which has 
the value 1.1 at initial yield and rises to a maximum of about 3 
for a fully plastic indentation. At full load the material beneath 
the surface on the z-axis will be in a stress condition given by 
the yield criterion: 

o,—0,=Y [C2] 


We now consider unloading as the superposition of a Hertzian 
distribution of surface tension whose maximum value is given by 
Eq. [C1] (with the sign reversed). If this process is elastic, it will 


1 The calculation has been repeated assuming a uniform dis- 
tribution of pressure at the end of first loading. This modifies 
Eq. [C5] to give c > 3.05 for reversed yield. While such a value 
is just in excess of the possible values of c, the actual pressure 
distribution is likely to be closer to Hertz than to a constant 
pressure. 
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produce an additional state of stress at the point r = 0, s = 0.47a, 
given by 


6, — 6, = —0.62p) = —0.93cY [C3] 


The final state of stress is therefore given by the superposition of 
[C2] and [C3] with the result: 


6,— 6, = (1— 0.93c) [C4] 


Now reversed yielding will occur if the value of 6,—«,. in [C4] 
exceeds — Y, i.e., if 


0.93¢>2 
c>215 [C3] 


Theoretical and experimental values of c as a function of load 
are given by Tabor (loc. cit.) in Figs. 28 and 29. For heavy loads 
it is clear that c exceeds the value given by Eq. [C5]. 

The above calculation is based upon an ideally plastic material 
or one which strain hardens isotropically. For a real metal, the 
Bauchinger effect would make reversed flow on unleading even 
more probable. 

It will be appreciated that the region of reversed flow will, in 
all probability, be small and contained below the surface. Its 
influence on the geometry of the indentation will be very slight 
which explains why it was not observed in Tabor’s experiments. 
Its importance is likely to arise only in repeated loading for a 
large number of cycles. 

3. Stress analysis in repeated loading. The authors have drawn 
attention to the probable significance of the complex cycles of 
stress in arriving at an explanation of the hardened and softened 
zones, the mechanism of cumulative plastic deformation and also 
of the nature of the failure. I would add that such a stress analy- 
sis should take into account plastic deformation and _ residual 
stresses introduced in the first loading. Though the analysis of 
the elastic-plastic problem presents considerable difficulties and 
can only be approximate, recent investigation (e.g., Merwin and 
Johnson, 1963) suggest possible methods of approach. 

4. Fretting and wear. From the work of Mindlin and Dere- 
siewicz [J. Appl. Mech. 327 (1953)] confirmed experimentally by 
Johnson [J. Mech. Eng. Sci. 3, 362 (1961)] fretting would not 
be expected at the interface of two similar metals loaded within 
the cone of friction. One is tempted to ask how closely the con- 
tact force in the present apparatus can be controlled to be axial 
under dynamic conditions. The case of dissimilar elastic solids 
leaded normally has been analyzed recently by Goodman [J. Appl. 
Mech. 515 (1962) and developed in a more recent unpublished 
paper]. 


AuTuHors’ CLOSURE: 


The authors are indebted to Mr. Tallian, Mr. Tawresey, and 
Dr. Johnson for their interesting and helpful comments. All 
three discussors have remarked on the relation between normal- 
contact fatigue and rolling-contact fatigue. We have stated in the 
paper that our objective was to obtain basic information on 
several aspects of contact fatigue under deliberately simple cir- 
cumstances, and no attempt was made to simulate actual rolling- 
contact bearing operation. The application of basic data to an 
actual piece of hardware must certainly take into consideration 
the differences in geometry, kinematics, etc. It is in this light that 
the parallel cited by Mr. Tawresey between our results and his 
results on oscillating bearings is particularly noteworthy. 

Mr. Tallian has said that no hydrodynamic action would be 
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expected in the normal oscillating contact. We feel, however, that 
hydrodynamic action is possible. A fluid with a pressure- 
sensitive viscosity can provide load support under oscillatory 
motion. Further, Dr. F. W. Smith, in his discussion of our com- 
panion paper (4), has pointed out that “the fluctuation of load 
at 26kc causes the contact zone to sweep from the diameter d, 
to d, with what is in effect a rolling motion with a peak velocity 
of the order of meters per second. Under these conditions one 
would expect a hydrodynamic film to form in the shaded annulus 
and to be absent from the center of the zone.” 

As to the traction resulting from differences in specimen curva- 
tures, we would like to note first that the differences in curva- 
ture are not requisite to the experimental method. Furthermore, 
we feel that one can be misled in cases of elastic deformation by 
looking at the differences in arc length and assuming that there 
must be a corresponding slip. 

He states that “the critical subsurface shear stress is the maxi- 
mum shear stress in the axis of symmetry, . . .” We feel that our 
microhardness data do not bear this out, and that further exam- 
ination of failure criteria under tri-axial stress is called for by 
this evidence. 


Other than to point out the above qualifying arguments, we 
concur with Mr. Tallian’s general theme. 


Turning now to the discussion by Dr. Johnson, we wish first to 
say that we have followed his analyses involving plasticity effects 
with great interest. As he has noted, both the bands of micro- 
hardness alteration and the continuing deformation under cyclic 
loading raise important questions as to the nature of plastic flow 
in the contact. We can only hope to make a closer tie between 
improved analysis of such effects and observed material behavior 
in the future. 


As to the question of how closely the contact force cculd be 
controlled to be axial under dynamic conditions, there was 
normally very little relative motion between the specimens in the 
lateral direction. We must point out that there were cases in which 
a slight relative motion occurred, and these were readily identified 
by examining the test specimens under a microscope. In the 
majority of the tests, no perceptible relative motion was noted, 
and the center area of the specimens, that is, the area that 
remained in contact at all times during testing, was very clean 
and undamaged even after numerous stress cycles. 
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Crater Wear and Failure of High Speed Steel Tools 
in Dry and Lubricated Cutting 


By A. DORINSON (asLE)! 


Detailed measurements were made of the growth of the wear crater in high speed tools when 
cutting AISI 1141 steel, dry and in the presence of a compounded cutting oil, at two speeds. The 
time-dependence of wear was studied by plotting the maximum cross sectional area of the crater. 
After a break-in period, wear assumed a constant rate which depended on the speed and the 
absence or presence of cutting oil. For a given speed, the wear rate for lubricated cutting was half 
the rate dry. Tool failure was found not to be an abrupt transition from orderly wear to burn- 
out, but rather a progression from orderly wear to catastrophically high wear without burn-out, 
and eventually to failure. These phenomena are discussed in their relation to the author’s previ- 


ously proposed theory of cutting tool wear (1). 


Introduction 


One of the means used to evaluate cutting oils is the 
lathe-turning tool-life test, where single-point tools are 
run to complete failure over a range of cutting speeds 
and the data so obtained are fitted to the logarithmic 
form of the Taylor equation 


log V + n log L = log C [1] 


In this equation, V is the cutting speed, L is tool life, 
n is an empirical constant whose value depends on the 
particular cutting conditions (a given workpiece, given 
tool geometry and material, given chip width, and given 
cutting fluid), and C is the cutting speed for unit tool 
life. The Taylor equation has been found very useful 
for quantitative comparison of cutting oils in spite of 
the fact that quite often the repeatability of the tool 
life determinations is only mediocre. 

The overt indications of tool failure are unmistakable; 
sparks are generated at the tool point, power consumption 
increases sharply, the surface finish of the workpiece 
deteriorates, and the cutting edge and the nose of the 
tool are completely eroded away. These are familiar 


Contributed by the ASLE Technical Committee on Fluids for 
Metalworking and presented at the Annual Meeting of the 
American Society of Lubrication Engineers held in New York, 
April, 1963. 

1 Section Leader, Fundamentals of Lubrication Section, Sinclair 
Research, Inc., 400 East Sibley Boulevard, Harvey, Illinois. 


phenomena. It is also well known that the motion of the 
chip across the face of the cutting tool quickly wears 
a crater in it. A reasonable presumption is that the 
ultimate end of crater growth is tool failure. The in- 
vestigation reported here concerns crater growth, its 
transformation into tool failure, and the effect of cutting 
oil thereon. 


Experimental 
LATHE-TURNING TESTS 


The turning tests were carried out on a single log 
of normalized AISI 1141 steel with high speed tools of 
the following geometry: 8-14-6-6-6-15-3/64. The depth 
of cut was 0.100 inch, the feed was 0.012 inch per 
revolution, and the cutting speeds were 138 and 150 fpm. 

To obtain a progressive picture of the wear at a 
given speed, a series of tools was employed, each tool 
cutting for a longer time than the preceding one. Some 
of the tools were run to failure by burn-out. The number 
of tools for each experiment and the spacing of the 
cutting times depended on the cutting speed and whether 
or not cutting oil was used. 

One set of experiments was run dry at the two cutting 
speeds; in the other, a cutting fluid was used. This was 
a compounded oil containing dissolved elemental sulfur, 
a sulfurized additive, a chlorinated additive, and a fatty 
additive, in a petroleum base oil. Some pertinent char- 
acteristics of this fluid are given in Table 1. 
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TABLE 1 
Characteristics of Cutting Oil 


Analysis: Sulfur (%) 1.15 
Chlorine (%) 0.29 
Viscosity (cs/100 F) 33.51 
Constituents: Elemental sulfur, sulfurized additive, chlo- 


rinated additive, fatty additive, and pe- 
troleum oil. 


EXAMINATION OF THE TOOLS 


When a tool was disengaged from the cut, the crater 
was found to be filled with adherent chip metal and a 
large built-up edge of chip metal was present. The rake 
face had a polished zone at the trailing border of the 
crater. The grinding marks on the tool were not obscured 
in this zone, indicating that the chip had rubbed on 
this portion of the face but had not worn it significantly. 
Bands of temper colors surrounded the crater. A brightly 
polished zone which did not obscure the grind marks was 
also present on the flank and the nose of the tool. Evi- 
dently there had been considerable rubbing contact with 
the shoulder of the cut and the newly generated workpiece 
surface, but this contact had not resulted in measurable 
wear. A band of temper colors was found on the borders 
of the bright zone. 

The first step in the examination of the crater was to 
hone away the adherent chip metal to the level of the 
tool face and then to etch strongly with 2% nital 
solution. Figure 1 shows the typical appearance of the 
crater thus exposed. The tool was mounted in a measuring 
microscope with its face perpendicular to the optical 
axis and the outline of the scar was mapped by means 
of rectangular coordinates. 

The next step was to section the tool perpendicular to 
the cutting edge so that the cross sectional area of the 
crater could be mapped. At least three cross sections were 
measured for each tool. Figure 2 shows a typical cross 


Fic. 1. View of honed and etched rake face. Original mag- 
nification 17 x. 
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Fic. 2. Cross section of crater. Original magnification 50 x. 
Note adhering chip metal in crater. 


section after etching with nital reagent to delineate 
the boundary between the crater and the adherent chip 
metal. Figure 3 shows plots of three typical cross 
sections of a tool, the depth being exaggerated with 
respect to the width by a factor of five. The true pro- 
portions of these sections are indicated by the dashed 
lines. 

The area of a cross section was evaluated by plotting 
on ruled paper, counting the squares, and multiplying 
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Fic. 3. Cross sections of crater at different locations along its 
axis. Vertical scale 5 times horizontal scale. Dotted lines show 
true proportions of sections. A at the heel; B in the middle; C 
near the nose. 
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by the scale factor. As can be seen from Fig. 3, the 
cross sectional area depended on where the crater was 
sectioned. In the present work, the greatest significance 
is attached to the region where the crater is widest and 
deepest. The cross sectional area in this region, here- 
inafter designated as the maximum area, was found by 
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Fic. 4. Cross sectional area of crater as a function of axial 
location. 


plotting the measured cross sectional areas as a function 
of their location along the crater, as shown in Fig. 4, and 
reading off the maximum value of the area from such a 
plot. 

EXPERIMENTAL RESULTS 


The time-dependent development of the crater was 
examined in a number of ways. Figure 5 shows two 
examples of the cumulative mapping of the outlines of 
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Fic. 5. Cumulative growth of rake face scar at 150 fpm. Dry: 
1, 0.25 minutes; 2, 0.50; 3, 1.00; 4, 2.00; 5, 3.00. Lubricated: 1, 
2.00 minutes; 2, 4.00; 3, 7.00. 


the craters on the rake face of the tool. When cutting 
dry at 150 fpm, at first the scar grows in all directions, 
but eventually the progress of the scar border towards the 
cutting edge and the nose of the tool comes to a halt 
and further enlargement of the scar occurs by progress 
of the trailing border and the border farthest from the 
nose. There has been a certain amount of simplification 
in the outlines of the scars shown in Fig. 5. The scar in 
Fig. 1 exhibits a winglike feature in the region farthest 
from the nose of the tool which has been omitted 
in Fig. 5. This winglike portion of the scar is very shallow 
and does not seem to play an important part in the 
growth of the crater or the failure of the tool by 
burn-out. Another aspect of the scar evident in Fig. 1 
but not shown in Fig. 5 is the groove worn by the edge 
of the chip. This, too, seems to exert little influence on 
the progressive development of the crater or the break- 
down of the tool. 

The second half of Fig. 5 shows the cumulative 
growth of the rake face scar at 150 fpm in the presence 
of cutting oil. On comparing the growth of the two scars 
for equivalent fractions of the cutting time (from 2 to 
7 minutes lubricated, 1-3 minutes dry), it is observed 
that the rate of growth is perceptibly less for the lubri- 
cated case. The major influence of the cutting oil seems 
to be not so much a uniform decrease in the rate of scar 
growth as it is an inhibition of further growth once the 
scar has passed a critical size. It will also be noticed 
that in the later stages of dry cutting the scar develops 
a protuberance in the nose region which is absent when 
cutting oil is used. 

Cumulative growth of crater cross sections is shown 
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Fic. 6. Cumulative growth of cross section at 138 fpm. Dry: 
1, 1.00 minute; 2, 2.50; 3, 3.00; 4, 4.50; 5, 6.00. Lubricated: 1, 
5.00 minutes; 2, 10.00; 3, 15.00. 
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Fic. 7. Course of crater wear. Dry: A, 150 fpm; C, 138 fpm. Lubricated: B, 150 fpm; D, 138 fpm. 4 : Failure by burn-out; Vv : 


catastrophic wear. 


in Fig. 6. The upper half of this figure shows the develop- 
ment of the crater when cutting dry at 138 fpm. At 
first the crater enlarges by movement of both the leading 
and the trailing borders of the scar on the rake face and 
by increase in depth. Eventually the advance of the 
leading border ceases, and further enlargement of the 
crater is by movement of the trailing border and by 
increase in depth. As wear progresses, the deepest portion 
of the crater moves toward the trailing border. The lower 
half of the figure indicates that for lubricated cutting 
at 138 fpm the increase of the crater depth becomes 
relatively more important than the movement of the 
trailing border of the rake face scar. 

Figure 7 shows the variation of the maximum areas of 
the craters with cutting time. Although each curve 
terminates in a linear portion, the over-all variation with 
cutting time is not the linear one predicted by Trigger 
and Chao (2) on theoretical grounds and actually 
observed by them for cemented carbide tools. On com- 
puting the wear rates from the slopes of the linear 
portions of the curves, it is found that the cutting oil 
employed halved the wear rate at both the cutting 
speeds studied, as shown in Table 2. 

The arrows pointing upward in Fig. 7 indicate the 


TABLE 2 
Effect of Cutting Oil on Terminal Wear Rate 


Cutting Ratio 
speed Wear rate (dry / 
(fpm) (10-5 in?/min) lubricated) 
150, dry 3.35 
150, lubricated 1.60 
2.09 
138, dry 1.35 
138, lubricated 0.68 
1.99 


times at which tools failed by burn-out; the numbered 
arrows pointing downward show the times at which 
cutting was stopped with tools which were subsequently 
found to have damaged cutting edges. Figure 8 shows 
views of the honed and etched rake faces and flanks 
of the tools corresponding to arrows 1 and 2 in Fig. 7. 
Tool 1 is characteristic of the beginning of breakdown; 
the land between the crater and the cutting edge is only 
barely eroded, and although the cutting edge appears 
badly damaged in the flank view, its original outlines 
can easily be reconstructed in a cross sectional diagram. 
The small bulge at the trailing border of the wear scar 
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Fic. 8. Catastropic wear at 150fpm, lubricated. Left: 7 minutes; (top) flank; (bottom) rake face. Right: 10 minutes; (top) 
flank; (bottom) rake face. Original magnification: rake face 4 x, flank 25 x. Arrows on rake face views show location of flank 


views. 


grows with continued cutting until eventually the wear 
crater and the cutting edge assume an appearance 
typified by tool 2. When the tool is in this condition, 
both the width and the depth of the crater are signifi- 
cantly greater than extrapolation from the smooth wear 
rate would indicate. Yet even with the tool in this 
state, the cutting process does not exhibit any of the 
gross phenomena associated with burn-out. 


The writer has published a treatment of cutting tool 
wear intended as a theoretical basis for the Taylor 
equation (1). One of the postulates in this treatment was 
that failure of the tool by burn-out is associated with 
a critical wear volume, or, in the two-dimensional ap- 
proach used in the present paper, with a critical cross 
sectional area. If the transition from orderly crater wear 
to burn-out were abrupt and catastrophic, then the 
upward-pointing arrows in Fig. 7 would indicate this 
critical cross sectional area. An average value of 21.1 
(+ 2.5) & 10~-° sq. in. is computed for the assumed 
critical area from the location of the arrows in Fig. 7. 


At first glance the assumption seems reasonably well 
confirmed by the evidence. 


However, one must also consider the state of affairs 
represented by the downward-pointing arrows in Fig. 7: 
breakdown of the cutting edge and rapid enlargement of 
the central portion of the crater. For instance, for 
lubricated cutting at 150 fpm, the first signs of break- 
down appear at approximately 54% of the time to burn- 
out; the area associated with this condition is 61% 
of the assumed critical area. The transition from orderly 
crater wear to burn-out is more complex than postulated 
in the writer’s previous treatment (1) and the proposed 
critical area criterion is beset with considerable doubt. 


If the critical area concept is tenable, the following 
equation applies: 
W/W =? /t [2] 


where W and W’ are the wear rates at two different 
cutting speeds and ¢ and ?¢’ are the corresponding failure 
times.” Failure time data for dry cutting (the only 


2 This relation follows from Eq. [18] of Ref. 1. 
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TABLE 3 
Tool Life, Dry Cutting 


Cutting 
speed Tool life 
(fpm) (min) Average 
150 4.78 
5.10 4.94 
138 11.41 
15.22 13.32 


data of sufficient scope available) are shown in Table 3. 
The ratio ¢’/t computed from these data is 2.69. The 
ratio W/W’ computed from the corresponding wear 
data in Table 2 is 2.48. The discrepancy between the 
two ratios is 9%. This suggests that the effect of cutting 
speed on the wear rate of the transition from orderly 
wear to failure is the same as for orderly wear. If this be 
true, then comparisons among various cutting fluids 
obtained in the orderly wear domain would also be 
valid in the catastrophic wear domain. 


Another indication of the complexity of the wear 
process is the fact that the outline of the scar on the 
rake face grows in length as well as in breadth. There 
must, therefore, be some curl across the width of the 
chip. Contact pressures and contact temperatures would 
be expected to be greater in the early stages of cutting 
than when the crater has grown somewhat. This, perhaps, 
accounts for the higher wear rate shown in Fig. 7 by the 
initial part of each curve. The nonuniform pressure and 
temperature distribution along the length of the crater 
probably persists throughout the cutting process, as in- 
dicated by the fact that breakdown of the cutting edge 


begins in the middle instead of at the nose or the 
heel of the scar (see Fig. 8). 

A clue to the nature of the transition from orderly 
wear to catastrophic wear can be found in the cumulative 
growth of the crater cross section as depicted in Fig. 6. 
Here it can be seen that after the location of the leading 
border of the crater has been stabilized, further growth 
of the crater by increase in depth means that the slope 
of the leading side becomes steeper, thereby decreasing 
the volume of metal in the land between the crater 
and the cutting edge. Conceivably this could alter the 
heat dissipation in the land sufficiently to cause break- 
down of the cutting edge. There would then be two 
separate domains of action for the cutting oil: one, 
orderly wear; the other, catastrophic wear leading to 
failure of the tool by burn-out. The influence of cutting 
oil on the rate of orderly wear has been demonstrated 
in this investigation. Elucidation of the influence of 
cutting oil on the transition from orderly wear to 
failure would require a close study of the details of 
this transition. 
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Wear Studies on Moldings of Polytetrafluoroethylene 
Resin. Considerations of Crystallinity 
and Graphite Content 


By JOHN F. LONTZ! and MILES C. KUMNICK* 


The wear of molded bearings of polytetrafluoroethylene resin has been investigated in relation to 
crystallinity and filler content (graphite) as factors effecting changes in selected mechanical prop- 
erties. In a homogeneous, unfilled resin molding, the rate of wear increases in direct proportion to 
the modulus in flexure and inversely with yield strain. A useful, simple expression is proposed 
to describe this relationship. In a heterogeneous, graphite-filled molding at a low filler content, 
the rate of wear initially approximates that of the homogeneous unfilled system and then decreases 
with time; at moderate filler content, the wear rate is constant at levels below that of the unfilled 
resin. The rate of wear is also proportional to the applied load at low loading under conditions 
where bearing temperature remains reasonably constant; the wear data indicate the importance 
of removing frictional heat especially at high bearing loads. Practical suggestions are included for 


design of efficient bearings. 


Introduction 


POLYTETRAFLUOROETHYLENE (PTFE) resins possess 
unique combinations of low coefficient of friction along 
with such important attributes as chemical, thermal, and 
electrical resistance, and mechanical endurance at low 
and high temperatures. Molded forms of the resin in 
various compositions are assuming recognized utility as 
bearings or solid lubricants in mechanical components 
of various designs (1-6). Although the literature is 
replete with studies on friction (7-12) and to some 
extent on wear (1, 6, 8, 13), an understanding of physical 
or mechanical properties and structural characteristics 
of molded PTFE in relation to wear has yet to be 
developed. From the lubrication engineering standpoint, 
wear or damage to the molded resin in bearing applica- 
tions may often be more important than frictional 
resistance itself. The fact that molded resins undergo 
deformation and wear under bearing loads has prompted 
numerous modifications with various fillers for increasing 
resistance to mechanical stress. It has been proposed 
(14-16) that wear is dependent upon the bulk mech- 
anical properties of the wearing surfaces. Since the 
mechanical properties of the molded or fabricated forms 
of PTFE resin are dependent upon (a) crystallinity 
(17) and (b) filler content (18), the present study is 
a preliminary attempt to develop appropriate correla- 
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tion of wear to these two factors so as to provide 
guidance on resin composition and bearing design. 

As semicrystalline materials, PTFE resins in molded 
form assume a range of crystallinity depending upon 
fabrication, especially the applied postsintering and an- 
nealing conditions. Figure 1 illustrates the effect of 
crystallinity (17) from 50 to 82% level of the basic 
resin, on a selected, limited set of mechanical properties, 
namely modulus in flexure, yield strain, and hardness. Of 
these the first two have been chosen for consideration of 
wear since they reflect the greatest fluctuation in meas- 
ured values and can be readily associated with deforma- 
tion under bearing load and with the yielding or tearing 
away of the resin as it forms the wear debris. There 
may be other physical and mechanical constants, in- 
cluding those involved in cyclic or repeated stressing, 
which may prove more amenable to correlation with 
wear with continued study. 

The marked effect of crystallinity on the mechanical 
properties indicated in Fig. 1 is in turn related to the 
organization of the submicroscopic crystalline structure. 
It has been recently shown (19) that semicrystalline 
PTFE gradually changes from a completely disordered 
structure into a well-defined banded structure with 
increasing crystallinity level, as indicated in Fig. 2. 
In effect, the amorphous portion gradually fills up with 
highly organized lamellae as the crystallinity increases 
to the indicated 89% level. The difference between the 
53 and 89% crystallinity levels is that the former is 
quenched or cooled rapidly during postsintering whereas 
the latter involves slow cooling or prolonged annealing at 
near melt (300C) temperature. Although deferred for 
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Fic. 1. Effect of crystallinity on selected mechanical properties 
[Thomas et al. (17)] (23.5C). 


later study, the manner in which the crystalline, lamellar, 
and noncrystalline or amorphous entities undergo rupture 
into wear debris should provide valuable clue on filler 
modification of resin molding for maximum wear 
resistance. 


Another practical means for modifying the resin is 
by incorporating fillers. In this study graphite was 
selected for several reasons, namely, it too possesses a 
lamellar structure with well-defined crystalline order, is 
nonabrasive, and is chemically inert compared to the 
innumerable fillers that are useful in enhancing various 
physical properties of PTFE moldings. Figure 3 shows 
the changes in flexural modulus, yield strain, and hard- 
ness in relation to graphite content. As with increasing 
crystallinity level, the modulus increases and the yield 
strain decreases with increasing graphite content. 

The experimental approach in this study involved first 
determining the wear rate using techniques applied by 
previous investigators (1, 8, 9) and secondly extending 
this to a more practical journal bearing type of test. This 
work is concerned principally with the tetrafluoroethylene 
(TFE) homopolymer, identified usually as TFE-fluoro- 
carbon resin, relegating that on the melt extrudable 
hexafluoropropylene copolymer, or FEP-fluorocarbon 
resin to a later consideration. 


Experimental details 
MATERIALS 


(1) Polytetrafluoroethylene. Test moldings were pre- 
pared principally from type III TFE resin described 
in the ASTM Designation D °1457-62-T. Commercial 
forms of type III in the form of molding powder and 
aqueous dispersions were used, the latter particularly for 
compounding with graphite. 

(2) Graphite. Fine-ground powder, made to 99% 
minus-200 mesh screen, was employed in compounding to 
10, 15, and 25 volume % filler based on combined 
fluorocarbon resin and graphite content. 


Fic. 2. Structures in PTFE resin at three levels of crystallinity [Starkweather (19) ] 
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Fic. 3. Effect of graphite content on selected mechanical 
properties (23.5 C). 


MOo.LpING PROCEDURE 


Table 1 summarizes the molding conditions for pre- 
paring test sheets (4 x 5 X 0.060 inch thick) for me- 
chanical properties and stock cylinders (2-inch diameter 
X 2-inch height) from which machined rings and journal 
bearings were made. Crystallinity levels were determined 
from gravimetrically measured density (17) and checked 
by infrared using the characteristic amorphous band 
absorption (20) calibrated by reference to X-ray data. 
The molded density from ambient cooling, indicated in 
Table 1, averaged 54 + 2% even though the cooling 
rate varied from 3- to 5-fold. To attain high (82%) 
crystallinity, the cooling rate was adjusted to 0.02-0.03 
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C per minute. The intermediate density (68%) was 
attained by annealing at 300-305 C for 96 hours. 


SAMPLE PREPARATION 


The test specimens were machined with liberal air 
cooling to high polish and cleaned with several rinsings 
of carbon tetrachloride to remove machine oil and other 
atmospheric contaminants. Prior to use, the machined 
specimens were dried in vacuum desiccator for a minimum 
of 24 hours and stored in sealed containers to avoid 
chance pick-up of atmospheric dust or oil owing to the 
highly electrostatic character of the resin. 


Test EQUIPMENT 


(1) Stationary pin and rotating ring. This test as- 
sembly is an adaptation of the arrangements used by 
Archard and Hirst (21) and Flom and Porile (8). The 
principal dimensions are listed in Table 2 and shown 
schematically in Fig. 4. The test involves measuring 
the wear at the resin to resin interface. The wear volume 
and wear rate were measured from weight loss of the 
stationary pin and converted to volume loss (removed) 
by dividing by the density at the respective crystallinity 
level and graphite content. No systematic attempt was 
made to recover the debris for identification and analysis. 


(2) Journal bearing. This test equipment is shown 
schematically in Fig. 5. The dimension and testing 
features are summarized in Table 3 in comparison with 
those used by Ricklin and Miller (1) and White (9). 
This test involves wear at the resin to metal interface. 
The rotating shaft is made of AISI 410 stainless steel 
machined to 0.7500 inch with a tolerance of 0.0005-inch 
maximum oversize and polished to a 3-5 winch finish. 
The PTFE bearings, machined to 0.090 wall, are encased 
in a mild steel ring using slot and key to fix the bearing 
to the housing. A thermocouple well, using a ;'; by 
4-inch deep boring, was drilled in the center of the 
shaft for measuring temperatures developed during the 
test. As with the pin and ring test, the extent of 
wear is determined from the weight loss of the resin 


TABLE 1 
Molding Procedures for Test Specimens 
TFE-Fluorocarbon Resin: ASTM-1457-62T, Type III Class 1 


Cooling conditions 


Preform Rate 
Resin Mold pressure Sintering Pressure (360-250 C) 
form Composition form (Ib/sq. in.) conditions® (Ib/sq. in.) (°C/min) 
Powder Homogeneous (unfilled) Cylinder 1000 In mold 700 10-15 
Homogeneous (unfilled) Sheet 1000 Free (open) None 3 
Aqueous dispersion? Heterogeneous (graphite) Cylinder 5000—10,000 In mold 2500-5000 10-15 
Heterogeneous (graphite) Sheet 5000-10,000 In mold 2500-5000 10-15 


@ 360C for 1 hour of resin temperature (thermocouple insert). 
» Type III, class 1 equivalent admixed with graphite and co-coagulated by stirring. 
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TABLE 2 
Rotating Ring and Stationary Pin Wear Dimensions and Test 
Conditions 
References 
Descriptive features (21) (8) This study 
Shaft (in.) 0.50-1.250 
Ring (in.) 0.9375 2.0 2.0 
Pin (in.) 0.2500 0.150 0.150 
Loads, 
Min. (lb) 0.1 0.04 0.2 
Max. (Ib) 22 0.46 2.2 
Speeds, 
Min. (it min) 370 2 374 
Max. (ft/min) 370 196 393 
Distance 
Min. (10° cm) 
Max. (10% cm) 6 0.9 16 
Temp. (° C) undetd. undetd. ambient 
@ See Fig. 12. 


LOAD (W) 
STATIONARY PIN 


ROTATING RING 


ROTATING SHAFT 


Fic. +. Stationary pin and rotating ring for wear testing 


HOUSING (CASE) 

TEST BEARING 
- STEEL SHAFT 
CLEARANCE 
SLOT AND KEY 


TORQUE 
_w DIRECTION 


LOAD 


Fic. 5. Journal bearing wear test assembly 


TABLE 3 
Journal Bearing Wear Dimensions and Testing Conditions 


References 

Descriptive features (1) (9) This study 
Bearing, 

OD (in.) 0.920 

ID (in.) 0.50 a 0.740 

Wall (in.) 0.090 

Clear. (in.) — 0.005 0.005 
Speeds, 

Min. (ft/min) 20 1 60 

Max. (ft/min) 111 118 215 
Loads, 

Min. (Ib) 120 44 8 
Max. (Ib) 44 70 
Temp. (° C) undetd. undetd. 80-1307 

@ See Fig. 12. 


bearing converted to volume loss by the appropriate 
density. The wear debris at first assumes the form of 
isolated, microscopic flakes such as shown in Fig. 6, 
in this case under a high load amounting to 7000 pounds 
per sq. in. applied against polished stainless steel bushing 
for short distance of travel. On prolonged wear the 
isolated flakes appear to be smeared out and overloaded 
with continuous flaky films. It will be seen that the flakes 
in the case of the resin debris range from 2 to 8p in 
diameter with an estimated thickness as low as 0.01 u 
based on the interference rings. With graphite com- 
position, a similar flaking can also be seen but with 
considerable carbon in predominance. 


Results and discussion 
PIN AND RinG TEST 


With increasing crystallinity of the fluorocarbon resin, 
the extent of wear as measured by the volume (V) 
removed from the pin increases as indicated in Fig. 7. 
A fairly linear dependence of wear volume on the sliding 
distance (S) is evident for all three of the crystallinity 
levels employed at a constant load. With varying 
load (W), the extent of wear as measured in 
terms of wear volume per unit of distance travelled 
(V/S) is proportional to the applied load, as shown 
in Fig. 8, where the high crystallinity level (82%) data 
indicate a slope of m = 1, or approximately thereof, in 
the low stress up to 100-gm load. At higher stresses 
there is a marked departure from this unit slope (82% 
crystallinity) suggesting that some undetermined factor 
is affecting the load effect through changes in some 
feature of resin strength or structure under high stress. 
It would appear that the same factors prevailed for 
the high loads (> 10? gm) in the tests by previous 
authors whose data points are also indicated in Fig. 8 
but with no identity as to crystallinity level or bulk 
mechanical properties. Since normal fabrication pro- 
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SLIDING RESIN + 
UNFILLED DIRECTION 25% GRAPHITE 
Fic. 6. Wear debris deposited on polished steel. Resin crystallinity, 54%; magnification, 500 x 
10 x 1073 1 T 
PIN AND RING Test | a 
| CRYSTALLINITY 
> 4 
7 
WwW 
54% < 62% 
ae = 
© 64% 
; [PIN AND RING TEST | © 56% 
CRYSTALLINIT Y 
10-0 at l 
10 10? 103" 


5 10 15 20 


SLIDING DISTANCE -cm 


Fic. 7. Wear vs. sliding distance at various crystallinity levels. 
Load, 100 gm; speed, 190 cm per second. 


cedures do not usually lead to such high levels as 82%, 
it is surmised that the high loads (> 107 gm) reflect 
thermal effects, as measured in succeeding journal tests, 
and attendant decrease in the mechanical strength. This 
suggests that for a discriminating study of wear of the 
fluorocarbon resin it is important to consider ambient 
conditions and structure at the wearing interface. 
The pronounced changes in mechanical constants 
of the fluorocarbon resin moldings and observed changes 
in wear rates in relation to crystallinity provide a start 


LOAD-gm 
Fic. 8. Wear rate vs. load at various crystallinity levels. Load, 
as indicated; speed, 190cm per second. @ Archard and Hirst, 
crystallinity undetermined; p Flom and Dorile, crystallinity un- 
determined. 


in developing some one form of an empirical relationship. 
Though any such relationship may be far from perfect 
in recounting all the necessary factors related to wear, 
it is nonetheless important from mechanical design stand- 
point to employ mechanical or elastic constants to monitor 
the molding bearing for minimum wear. Just as elastic 
constants, both in static and dynamic sense, are proving 
most enlightening in understanding for instance frictional 
effects (22, 23), so the mechanical constants quite 
properly should be considered in accounting differences in 
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the extent of wear. The two most prominent constants 
chosen in this consideration are (a) the flexural modulus 
(E£), since it is reasonable to conceive of circumferential 
flexing as a mode of stressing of the pin where the 
load position represents the fixed end of a short flexing 
cylinder and the contacting wearing surface represents 
the deformed end, and (b) the yield strain as this 
provides the onset on the initial rupture from the bulk 
molding to form the wear debris. From these two material 
constants, the wear rate, expressed in terms of the 
volume removed per unit of sliding distance (V/S) 
from qualitative inspection of the data in Table 4, 


TABLE 4 
Selected Mechanical Properties and Wear Rates vs. Crystallinity 
Pin and ring test 


Crystallinity Modulus Yield Wear/distance (V/S) 
level flexural strain 10° cm3/cm 
(%) (Ib/in*¢) (in/in) Exptl. Calc.? 
54 72,000 0.36 0.377 0.378 
64 118,000 0.24 0.752 0.770 
82 170,000 0.09 1.786 1.786 


® Converted to gm/cm*? by X 70.31 for Eq. [1]. 
» From Eq. [1], x = 0.961 > 1; y = 0486 > 0.5; K = 4.503 X 
10—-9, MKS units, 


appears to be directly proportional to the modulus and 
indirectly proportional to the yield strain. Applying 
separate exponents for these two material characteristics, 
the following general expression can, with discretion or 
limitations, be considered: 


V/S = K W (E*/e") [1] 


where K is a proportionality constant, W the applied 
load, and x and y are the appropriate exponents, which 
like the constant (K) assume values characteristic of 
the experimental conditions, the undefined but constant 
resin characteristics, and the nature of the mechanism 
of wear. Numerical solution of Eq. [1] for the three 
levels of crystallinity listed in Table 4 provides values 
of x = 1 and y= 0.5; these in turn lead to fairly good 
agreement between the experimental and calculated values 
for the wear rate, expressed as the V/S quotient. A 
suitable test of Eq. [1] would consist of extending the 
experiment to one or two orders of magnitude higher 
loads (10*-10' gm) while at the same time measuring 
the increasing interfacial temperature at the wearing 
surface and the corresponding flexural modulus and yield 
strengths. 

Effect of graphite. The admixture of graphite, thus 
giving rise to a heterogeneous system, markedly reduces 
the wear in general agreement with observations of 
previous investigators (1). The results are indicated 
graphically in Fig. 9, in which two interesting effects 
are evident. First, the homogeneous, nonfilled resin in- 
itially undergoes rapid wear, up to 1 & 10° cm of sliding 
distance, followed by a markedly lower wear rate. 


T T T T 
PIN AND RING TEST | 
or CONTROL (NO GRAPHITE) 
z 
10% GRAPHITE 
[=) 
24 
3 
25% GRAPHITE 
L 
10 15 20 


SLIDING DISTANCE -—cm 


Fic. 9. Wear vs. sliding distance at various graphite contents. 
Crystallinity, 54%; test load, 1000 gm; speed, 200 cm per second. 


Second, the 10 and 25% graphite content decrease the 
extent of wear but with a noteworthy difference, namely, 
the 10% graphite content initially approximates the 
extent of wear observed with the nonfilled control, 
whereas the 25% content eliminates this initial, high 
wear completely. From the practical standpoint, this 
observation suggests that bearing dimensions for graphite 
content less than 25 volume % should allow for greater 
wear-in than that applied, if at all, to the 25 volume 
% graphite content. The factor contributing to the 
initial high rate in the case of the nonfilled and the 
10% graphite content followed by the marked change 
in slope is not altogether clear. It is conjectured that 
some temperature-dependent effects, including changes 
in the E*/e” ratios and possibly hysteresis response, play 
an undefined role in offsetting the initial particulate 
rupture into wear debris. 

The application of Eq. [1] to the graphite content 
based on experimental and calculated values for the 
wear volume/distance (V/S) ratio is indicated in Table 
5. In this series only the homogeneous, nonfilled resin 
tested at a tenfold higher load provides a good agreement 
with the calculated V/S value, while those of the 10 and 


TABLE 5 
Selected Mechanical Properties and Wear Rates vs. Graphite 
Pin and ring test 


Wear/Distance (V/S) 


Graphite Modulus Yield 
content flexural strain 10° cm#/cm 
(%) (Ib/in?¢) (in/in) Exptl. Calc. 
0 72,000 0.36 3.85 3.78 
10 96,000 0.20 0.280 (6.88) 
25 140,000 0.04 0.140 (2.24) 


@ Converted to gm/cm? by xX 70.31 for Eq. [1]. 
» From Eq. [1]. 


25% graphite content, following the initial wear-in period, 
are at least an order of magnitude lower than demanded 
by Eq. [1]. Thus, the wear-resistant effect imparted by 
the lamellar graphite is emphasized quite significantly. 
It would appear that the supporting or protective matrix 
of the graphite overrides the effect of the selected 
mechanical constants. The data in Table 5 reduce to a 
general empirical relationship 


or «(Vp [2] 


where V, represents the filler volume and V, the resin 
volume. Thus a 2.5-fold increase in graphite content 
imparts an approximately equivalent (2.8) decrease 
in wear as the data in Table 5 indicate. 


JOURNAL BEARING WEAR TESTS 


The wear-resistant effect of graphite provided the 
basis for examining the 25% content in the more practical 
journal bearing tests with a view to ascertaining the 
effect of load and especially temperature on bearing wear 
and design. Figure 10 indicates the extent of wear (V/S) 
observed at various loads involving two test speeds. The 
corresponding apparent increase in radius, calculated 
from the wear volume removed from the starting bearing 
dimensions, is indicated at right portion of graph to 
provide an estimate of radial changes that could develop. 
The load effect on wear rate for the journal bearing 
tests on the 25% graphite composition, plotted in Fig. 11 
as an extension of the pin and ring test, appears to 
provide a slope of unity, with some deviation to be sure, 
in which V/S«W" with n= 1. Despite the deviation, 
the load effect in this case does not appear to provide 
too serious departure from the pin and ring to the 
journal test conditions. The load does in fact cause a 
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Fic. 10. Wear vs. sliding distance; 25% graphite content, 
54% crystallinity. Test speed: © 30cm per second; @ 110cm 
per second. 
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PIN AND RING TEST| | | JOURNAL BEARING 
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~ 
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Fic. 11. Wear rate vs. load comparing pin/ring and journal 
bearing tests. 


JOURNAL BEARING TEST | 
80 a 
4 12 16 20x 105 
LOAD - gm 


Fic. 12. Effect of load on bearing wall temperature; 25% 
graphite; speed, 110 cm per second. 


marked rise in temperature as indicated in Fig. 12. 
This in turn affects markedly the mechanical properties 
(18) and should account for the deviation. 


THERMAL CONSIDERATIONS IN DESIGN 


The design of the journal bearing in this test was 
based on a 0.090-inch wall molded component. This 
raises the question as to whether this thickness is 
acceptable for solid lubricant bearing surfaces. One 
approach to this query is to consider, as shown graphi- 
cally in Fig. 13, the thermal gradient and thermal flow 
and design the bearing-case to assure operating tempera- 
tures consistent with Fig. 12. Another would be to con- 
sider redesigning the interfacial wearing surface so as 
to provide cooling channels, helically grooved at an 
optimum pitch and angle to draw off frictional heat. 
Both approaches are shown graphically in Fig. 13. 
Other modifications with borings in the housing case, 
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(A) REDESIGNED TO: 


SHAFT 

BEARING 

GROOVES 

CASE OR HOUSING 


Fic. 13. Bearing design for ambient cooling. (A) Decreasing 
temperature with decreasing bearing wall. (B) Cooling of bearing- 
shaft intersurface with groovings. 


shafts, and other stationary and rotating components 
for efficient removal of heat can provide a variety of 
useful design features. Some additional study needs to 
be given toward determining thermal gradients and 
dynamic friction developed in the course of the wearing 
process. With this type of information at hand, the 
selection of appropriate composition for general me- 
chanical services involving wear can be approached in 
a manner less dependent upon costly and indeterminate 
trial-and-error methods. 


Conclusions 


The extent of wear with moldings of PTFE resins in 
homogeneous, unfilled moldings is directly proportional 
to increasing crystallinity content. With a heterogeneous 
system using graphite as filler, the wear is reduced in 
proportion to increasing graphite content. An approxi- 
mate relationship can be applied to equate wear with 
resin crystallinity and graphite content based on selected 
mechanical constants, in particular flexural modulus and 
yield strain, on the one hand and graphite content on 
the other. Temperatures developed during the wear 
period may have to be considered in relation to changes 
in the mechanical constants and in turn for developing 
most efficient bearing design. 
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DISCUSSION 


C. J. SPEERSCHNEIDER AND C. H. Li (Honeywell Research Center, 
Hopkins, Minnesota) : 


The authors’ equation relating the wear rate to the modulus 
in flexure and yield strain provides a basis for predicting wear 
of PTFE. Their relationship, kw E*/e¥, introduces three 
unknowns k, x, and y, and the experimental observations of 
Fig. 4 provide three equations by which the unknowns can 
be determined. The experimental and calculated values of V/S 
must agree as shown in Table 5 by virtue of the simultaneous 
solution of the three equations. It seems advisable to have at 
least one more set of data to verify this important relationship. 
This might be one of the reasons why their equation does not 
apply to filled PTFE as indicated in Table 5. 

The authors suggested that microscopic observations of the 
surface were most useful in studying wear of filled PTFE. In the 
case of spherical Al,O, filled material, we also found that 
microscopic observations were very useful, and we were able 
to observe the increasing sphere concentration as steady state 
wear was approached. Once the steady state rate was reached, 
the surface had become saturated with filler, the initial concen- 
tration having little effect on the final concentrations. In addition, 
using the electron microscope, it was possible to determine whether 
the filler particles were coated with PTFE. 

Work of this nature, as conducted by these authors, is ex- 
tremely important, both for the prediction and comparison 
of wear. Possibly the next step in an analysis such as the 
authors have presented would be to relate the properties (E 
and ¢) to microscopic fracture involved in wear. 


R. C. Bowers 


The authors are to be commended for aceomplishing a 
difficult assignment, the reproducible measurement of wear and 
the reduction of some of the data to a simple equation which 
describes the rate of wear in terms of certain physical properties 
of the material. 

The equation derived by the authors, in which the wear rate 
of unfilled TFE at low loads was found to be dependent merely 
upon the modulus in flexure and yield strain, is quite remarkable 
when the complexity of the wear process is considered. It 
is of interest to speculate on the nature of wear which would 
be dependent only upon these properties. Four types of wear 
are generally recognized: (a) adhesive, (b) corrosive, (c) 
abrasive, and (d) fatigue. With most sliding solids including 
plastics, wear is most commonly of the adhesive type which 
is caused by the “welding” of contacting asperities and the subse- 
quent subsurface shearing of the weaker material. However, 
TFE, because of its low specific adhesion, is an exception, as 
evidenced by its anomalously low coefficient of friction, at least 
at low sliding speeds and loads. The chemical inertness of TFE 
precludes corrosive wear. Abrasive wear rapidly diminishes as 
the surfaces are “run in” if the system is kept free of extraneous 
abrasive particles. Wear by fatigue is usually associated with 
rolling friction but with the apparatus used the rhythmic loading 
and unloading of the TFE could induce subsurface cracks which 
eventually result in wear particles. It is because the wear of TFE 
in these experiments is predominantly of the fatigue type that 
it can be related to the modulus in flexure and yield strain. 
However, with most sliding surfaces wear due to fatigue is 
usually small compared to adhesive, abrasive, or corrosive wear. 
TFE is exceptional in this respect. 


The observed increase in wear rate of the unfilled TFE with 
increasing crystallinity is of particular interest. We have observed 
that the friction of many polymers decreases with increasing 
crystallinity. These facts again indicate that the wear at low loads 
was not adhesive since this type of wear would increase with 
increasing friction. 

The nonlinear increase in wear rate at loads in excess of 
100 gm has been attributed to a decrease in yield elongation as 
a consequence of the higher surface temperatures encountered 
at these loads. This temperature rise could also cause sufficient 
thermal softening to permit more intimate contact and greater 
adhesion. At the higher loads the total wear may be the sum 
of both adhesive and fatigue wear. It would be of interest 
to know if the surfaces appeared more pitted at the lower loads 
than at the higher loads. 

The authors stated that most of the results are concerned 
with TFE. Apparently at least one measurement was made 
with a filled composition of FEP (a copolymer of TFE and 
hexafluoropropylene). It cannot be assumed that these two 
resins will have the same wear rate. We have found that 
at low sliding velocities the kinetic coefficient of friction with 
FEP is greater than that of TFE by a factor of three. Have 
any wear measurements been made for both resins under 
comparable conditions ? 

The wear rate of graphite-filled compositions (Fig. 10) was 
greatest at the lowest load and least at the highest load when 
the sliding velocity was 110cm per second. Wear rate was also 
much larger when the load was 7300gm and the speed 30cm 
per second than it was for a 7700 gm load at 110cm per second. 
This contradicts the argument that where experimental conditions 
promote greater surface temperature, greater wear results. Would 
the authors comment on this? 

Many fillers have been successfully used to decrease the wear 
of TFE. These include such diverse materials as glass fibers, 
asbestos, molybdenum metal, bronze, iron oxide, and a synthetic, 
colloidal, hydrated aluminum oxide. The resulting compositions 
are much harder than the unfilled resin. It is rather curious that 
graphite was the one material chosen in this study as the 
filler. Graphite is most generally used as a friction reducing 
additive, such a material is not needed with TFE. Was there 
a particular reason for selecting graphite? Have any other 
fillers been examined ? 


AuTuHorS’ CLOSURE 


The suggestion by Speerschneider and Li on the possible 
extension of Eq. [1] to heterogeneous systems, such as that 
containing graphite described in the present paper or aluminum 
oxide to which they refer, provides an attractive sequel on 
fundamental aspects of wear which has yet to be pursued by 
the authors. 

It is conceivable that some more general and definitive expres- 
sion than Eq. [2] can be developed within limited range of 
filler volume and only under fairly narrow conditions. This 
is because the wear rate, for instance at low filler loading, 
undergoes a breaking-in period from which it levels off to 
steady-state wear. Hence we are confronted with a transition 
effect. Then we have the case with continued increase in the 
filler content beyond equal volume ratios where the filler exceeds 
the resin. Here we should expect a different and possibly higher 
wear rate as the bonding of the filler to the resin is notoriously 
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poor in contrast to the low filler content (25 volume %) where 
the filler is substantially contained in the resin. 

The comments by R. C. Bowers are indeed quite illuminating 
and definitive in delineating the four specific types of wear 
with implications much broader than our specific area with 
TFE-fluorocarbon resin. The suggestion that at higher loads the 
total wear may be the sum of both adhesive and fatigue wear 
must remain unanswered at the moment because no close check 
was made on surface appearance of the worn resin for pitted 
structure. 

While no systematic examination of the wear of FEP-fluoro- 
carbon resin (copolymer of tetrafluoroethylene and hexafluoro- 
propylene) was made in relation to crystallinity and graphite 
content, there is evidence of higher wear rate (approximately 
3- to 4-fold) from a single journal bearing test with a 10 
volume % graphite composition compared to an equivalent 
TFE-graphite composition. This would be consistent with the 
commenter’s observation of a 3-fold difference between the 
kinetic coefficient of friction of FEP and TFE resins. 


The marked difference in the wear rates shown in Fig. 10 
at different speeds (30 vs. 110cm per second), especially as one 
compares the data at nearly constant load (7300 and 7700 gm 
respectively), presents an anomaly as the commenter indicates. 
The authors have no phenomenological explanation, such as could 
be used as a rationale, based either on physical inspection and 
characterization of the worn surface or more importantly on 
the removal or entrapment of the wear debris in relation to the 
applied journal speed. The authors would look to differences either 
in elastic response at small strains or hysteresis effects which 
would initiate a chain of mechanisms leading to higher wear 
at the lower speed, as a pure speculation. 

Graphite was selected in this study because of its lamellar 
structure similar to that of the lamellae of the highly crystalline 
TFE resin and more particularly because of its chemical inertness 
and nonabrasive character compared to several fillers listed by 
the commenters. However, their implied suggestion warrants 
continued studies with different fillers that can impart a new 
set of mechanical properties. 
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Lubrication by Low-Melting-Point Metals at 
Elevated Temperatures 


By M. IMAI! and E. RABINOWICZ (AsLeE)* 


Measurements of the friction coefficient as a function of temperature have been carried out 
using surfaces of stainless steel covered by lead, tin, zinc, cadmium, and Wood’s metal applied 
in powder form. Some work has also been done with a few other metal and nonmetal surfaces. 
In cases where the interaction between the low melting substance and sliding surfaces is high, as 
revealed by the occurrence of wetting, the friction reaches a peak just below the melting 
temperature of the substance, and then drops to considerably lower values just above the 
melting point. The peak below the melting temperature is associated with the formation of 
large adhering fragments of the low melting substance on the sliding surface. When there is no 
wetting, the low-melting substance has, either below or above its melting point, essentially no 


effect on the friction. 


Introduction 


A STANDARD method of studying the mode of action of 
lubricants is to measure the friction of lubricated surfaces 
while the temperature is raised. This method dates back 
to 1922, when Hardy and Doubleday (1) found that 
when organic materials were used as lubricants, there 
was a sudden increase of friction at the melting point of 
the materials and that the friction remained high there- 
after. This work has been often repeated and extended, 
most prominently by Bowden et al. (2), who showed 
that, if the lubricant could react with the surfaces, then 
it was at the melting point of the products of reaction that 
the rise in friction was observed. 

A second class of materials which has been studied 
in detail consists of glass-forming materials. It has been 
shown by Peterson e¢ al. (3) and by Rabinowicz and 
Imai (4) that with these materials there is a pronounced 
peak in the friction at the softening point of the glass, 
with the friction much lower on either side of the soften- 
ing temperature. 

In view of the fact that these two types of materials 
behave so differently, it was decided to study other 
materials which could possibly be used as boundary 
lubricants, to see how their frictional behavior varied 
at their melting point. 

Most of the authors’ work has been done using various 
low melting metals as lubricants, which by virtue of 


Contributed by the ASLE Technical Committee on Lubrication 
Fundamentals and presented at the Annual Meeting of the 
American Society of Lubrication Engineers held in New York, 
April, 1963. 

1 Research Assistant in Mechanical Engineering, Massachusetts 
Institute of Technology, Cambridge, Massachusetts. 

2 Associate Professor of Mechanical Engineering, Massachusetts 
Institute of Technology, Cambridge, Massachusetts. 


their simple chemical structure are not subject to thermal 
decomposition or degradation. Unfortunately, low-melting 
metals do react with an air environment, and this does 
produce some difficulty. 

There have been but a limited number of studies 
reported of the lubricating properties of metals. Thus, 
Vinogradov et al. (5) have made observations on the 
use of liquid Wood’s metal and mercury as lubricants, 
while Campbell (6) has investigated liquid sodium, and 
Coffin (7) has studied liquid NaK, the sodium-potassium 
eutectic. In none of the cases, apparently, were observa- 
tions made at the melting point of the metals. In addi- 
tion, Umeda and Nakano (8) and Kingsbury and 
Rabinowicz (9) have studied the frictional properties 
of unlubricated uniike metals which form eutectics. In 
this case, a striking drop in the friction was observed as 
the temperature was raised to the eutectic melting 
temperature. On the other hand, a few observations 
reported by the authors (4) were not so clear cut. 
Used as lubricants on stainless steel, tin and cadmium 
gave an increase in friction at their melting temperature, 
while lead and zinc showed a decrease. Clearly, this 
situation deserved further investigation. 


Apparatus 


The friction apparatus used in this study utilizes the 
geometry of three pins of one material contacting a 
rotating plate of the other material, the region of contact 
being inside a furnace. The plate is mounted on a shaft 
which extends outside the furnace and is connected 
through a pulley system to a variable-speed electric 
motor. The pins are mounted in a self-aligning holder 
which is connected to another shaft, the far end of 
which is outside the furnace and is restrained from rota- 
tion by means of strain gage rings. The loading is by 
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Fic. 1. Schematic illustration of high-temperature friction ap- 
paratus. 


means of weights applied to this shaft. The furnace is 
made of welded steel plate and is air-tight except for 
the front door, which is bolted on and supplied with a 
cooling system to protect the polychloroprene gasket, 
and the two holes for the shaits, which are covered by 
rotating seals and also cooled. The chamber is electrically 
heated, using resistance wires. A schematic drawing of 
the apparatus is shown in Fig. 1. Further details are 
given by Mamin (10). 

For the runs described in this paper, the temperature 
was gradually raised from room temperature, during 
which time sliding was continuous and the friction was 
monitored by an electrical recorder. The temperature 
was measured by a thermocouple connected to one of 
the pins. The pins were of 4-inch diameter and were 
normally given a slightly rounded end (radius about 2 
inches). The flat specimen was of dimensions 2.25 
2.25 0.5 inches with a recess of 2-inch diameter and 4- 
inch depth (Fig. 2). Low-melting metal powder of commer- 
cidl purity was applied to this recess in excess quantity so 
that the effect of the oxidation of the metal powder could 
be minimized. This preparation differs from that in 
previous tests, in which limited amount of the soft 
metal, in powder form, was applied to the flat surface. 
The electric contact resistance of the contacts was also 
monitored in some of the tests in order to detect both 
physical and chemical changes occurring on the sliding 
surface. This measurement was made using a Wheatstone- 
bridge circuit, in which the electrical contact resistance 
was shunted across a 10 ohm resistor, and the signal was 
fed into an electrical recorder of time constant 14 second. 
Atmospheric control of the furnace was carried out 
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Fic. 2. Schematic illustration of the test geometry and the 
test specimens. 


by replacing the air with argon of commercial purity. 
In this process, the argon was slowly introduced from 
its cylinder into the furnace while the air was pumped 
out by a vacuum pump. After 15 minutes, the pump was 
shut off and the friction run was started. During the 
test, argon was continuously fed into the furnace at a 
slight positive pressure in order to prevent back flow 
of air into the furnace. 

In all friction tests, there was a break-in period of 
sliding at room temperature for a few minutes until the 
friction reached a steady-state value. Then the tempera- 


. ture was raised. 


In almost all cases, the repeatability of the friction 
data was quite good. The level of the friction coefficient 
on repeat testing showed a variation of perhaps + 10%, 
and the main features (such as peaks, dips, etc.) occurred 
at temperatures within + 5C at low temperatures, 
+ 25C at high temperatures. 


Results 


FRICTIONAL PROPERTIES OF LOoW-MELTING PURE METALS 


Friction data, using stainless steel riders on a stainless 
steel flat covered with powder of the low-melting metal, 
are shown in Fig. 3. It is seen that atmosphere has some 
effect on the friction of these metals. This indicates that 
the results were affected by the formation of oxides during 
the tests. Hence, for purposes of comparison, tests were 
also made using the oxides of the metals. These results 
are also plotted in Fig. 3. 

The friction-temperature plots of lead and of lead 
monoxide are quite similar, giving a friction coefficient 
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Fic. 3. Friction-temperature plot for stainless steel on stain- 
less steel (grade 303) with metal powders. Load, 1700 gm; speed, 
0.1 cm per second. 


of about 0.5 at room temperature and of 0.15—0.20 
at higher temperatures. A possible explanation is that 
lead was oxidized rapidly during the test, so that the 
plot obtained in the test using lead was actually that 
characteristic of lead monoxide (PbO). It is considered 
likely that lead was oxidized even in argon, owing to 
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residual traces of oxygen impurity in our apparatus. 
The other oxides of lead, namely the red oxide (Pb3;0,) 
and black oxide (PbO2) gave somewhat higher friction, 
as shown in Fig. 4. However, the friction of these oxides 
also showed a gradual decrease as the temperature was 
raised. The decrease is probably due to the formation of 
lead monoxide (PbO) since the latter is the main form 
stable at higher temperatures. Lead monoxide melts at 
888 C, and all curves show an increase in friction and a 
decrease in electric resistance as this temperature is 
approached. 

The friction plots of tin and stannic oxide (SnO2) 
sliding in air are also similar, and this suggests that 
tin was oxidized during this test. However, the friction 
plot of tin in argon is quite different from that in 
air. It showed a peak below the melting point and a 
very sharp decrease at the melting point. The test was 
stopped at 300 C, and a photograph of the sliding surface 
was taken after the specimen was brought down to room 
temperature (Fig. 5). It may be seen that material in 
metallic form is adhering to the sliding surface, although 
a portion of the metal powder had become oxidized during 
the test. 
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Fic. 5. Sliding surface of stainless steel flat specimen covered 
with tin, at 300 C. 


The friction plots of zinc and cadmium also gave a 
peak below the melting point, while the oxides of these 
metals showed little change at this temperature. The 
examination of the sliding surface after the tests in 
argon revealed that zinc and cadmium remained in 
the metallic form, as in the case of tin, although again 
a portion of the metal powder had become oxidized during 
the test. 


There is a distinct difference between the friction- 
temperature curve of lead and that of the other metals 
shown in Fig. 3. As was discussed above, lead was 
presumably oxidized rapidly, and the curves using lead 
as a lubricant merely show the frictional properties of the 
oxides. On the other hand, the other metals remained 
metallic above the melting point, at least in argon. 
The friction-temperature curves of these metals show 
a peak below the melting point and a sharp decrease at 
that temperature. The friction above the melting point 
becomes lower than that of the oxides (or of unlubricated 
stainless steel), suggesting that the molten metal has 
some lubricating capacity. The oxides of these metals, 
as might have been expected, did not produce any change 
in frictional properties at or near the melting point of the 
metals. 


In order to generalize the above results, tests were 
also made using titanium surfaces. Figure 6 shows the 
friction-temperature plots for titanium surfaces covered 
with cadmium and tin. The plots again give a peak and 
a dip at the melting point. 


The above results show that metal powders give a peak 
in friction just below their melting point, if the metals 
have not become oxidized. Above the melting point, the 
molten metal has some lubricating capacity. In order to 
study these effects under conditions which could be more 
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Fic. 6. Friction-temperature plot for titanium on titanium 


with metal powders. Load, 1700 gm; speed, 0.1 cm per second, in 
argon. 


readily controlled, it was decided to carry out another 
series of tests using a metal of still lower melting point, 
namely, Wood’s metal. 


FRICTIONAL PROPERTIES OF Woop’s METAL 


The eutectic temperature of Wood’s metal (an alloy 
of cadmium, tin, lead, and bismuth) is 69-72 C, and this 
makes it possible to study phenomena associated with 
its melting in a temperature range near room tempera- 
ture, before the effects of oxide formation or of thermal 
degradation of sliding surfaces become serious. Tests 
were carried out using Wood’s metal powder applied to 
metallic and nonmetallic systems; the results are shown 
in Fig. 7. The solid lines are the plots obtained simply 
by applying Wood’s metal powder to the sliding surface. 
However, in some cases, Wood’s metal formed droplets 
on the sliding surface and did not cover the surface 
completely. In order to ensure the presence of enough 
Wood’s metal to fully cover the surface, tests were made 
in which the surface was flooded with molten Wood’s 
metal to a depth of 4-inch through the addition of a 
chunk of Wood’s metal just above the melting point. 
These flooded tests were made while the specimen was 
heated and then cooled, and mean values of the results, 
obtained respectively during heating and cooling, are 
plotted in dotted lines. Also, the frictional properties of 
the unlubricated sliding surfaces without Wood’s metal 
are shown in the same figure for comparison. 

Copper on copper with Wood’s metal gave a peak 
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below the melting point. Above the melting point, the 
friction dropped to 0.2 and the low friction persisted up 
to 200C, at which temperature the test was stopped. 
Stainless steel on stainless steel with Wood’s metal 
also gave an increase and decrease in friction when 
passing through the melting point, but the friction went 
up again as the temperature was raised further. The 
phenol formaldehyde resin system with Wood’s metal 
gave a friction curve similar to that of stainless steel but 
with the maximum and minimum less pronounced. It 
is noted that the friction curve for these three sliding 
surfaces is quite different in the absence of Wood’s metal. 
On the other hand, the nylon on nylon system was 
quite unaffected by the presence of Wood’s metal, showing 
neither the increase nor the decrease at the melting 
point of the latter. This result is important, because it 
shows that the peak just below the melting point is not 
a property only of the soft metal, but that it is strongly 
dependent on the nature of the sliding surfaces. 


INTERACTION OF Woop’s METAL AND THE SLIDING 
SURFACES 


Interesting changes were noted in the physical form of 
the Wood’s metal film on the sliding surface during the 
course of the tests. Figure 8 shows the sliding surface 
of phenol formaldehyde resin and nylon at 100C and 
200 C, respectively. It is seen that the sliding surface of 
phenol formaldehyde resin is covered by a film of Wood’s 
metal at 100C, but the film is broken at 200C. A 
similar change in the film formation was also observed in 
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Fic. 7. Friction-temperature plot for metal and nonmetal sliding systems with 
Wood’s metal. Load, 1700 gm; speed, 0.1 cm per second, in air. 


the test using stainless steel systems. However, the 
copper surface was covered by the film both at 100C 
and at 200C, while as is shown in Fig. 8, the nylon 
surface was never covered by the molten Wood’s metal 
film. 

These changes in the appearance of the Wood’s metal 
film were paralleled by changes in the friction. Thus, 
the friction of the stainless steel and phenol formaldehyde 
resin systems was low at 100C but increased at 200C. 
The friction of copper was always low above the melting 
point of Wood’s metal and that of nylon system was 
not affected by the presence of Wood’s metal at all. 
Hence, it would seem that there is a very good correlation 
between the friction and the film formation. 

It might be argued that the change in friction is 
simply due to the physical disappearance of the Wood’s 
metal film and the consequent exposure of unprotected slid- 
ing surfaces. However, the tests with the flooded surfaces 
show that the same friction changes occur under flooded 
and unflooded conditions, and from this we may deduce 
that the variation in friction is caused not simply by 
the presence or absence of the molten Wood’s metal 
film but by the change in wettability of the liquid-solid 
system. This change can produce two independent effects, 
namely a change in the physical form of the film, and 
a variation in the friction coefficient. 

It follows that there is a very good correlation between 
frictional behavior and wettability; when the molten 
metal shows poor wettability on the solid surface, the 
friction curve shows no peak and no drop (e.g., nylon) 
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Fic. 8. Sliding surface of phenol formaldehyde resin and nylon flat specimens covered 


with Wood’s metal, at 100 and 200 C. 


but when the wettability is better, as with the other 
materials, the friction curve shows both a peak and a 
drop. 

FRICTIONAL PROPERTIES OF GLASS 


In order to generalize these conclusions, it was decided 
to carry out tests with other liquid-forming solids. 
Earlier work (4) has shown that glass-forming materials 
also show a peak and a drop in the friction at a tempera- 
ture close to the softening point of the glass. However, 
no tests to investigate the dependence of this phenomenon 
on wettability between the glass and the sliding surface 
appear to have been carried out. Therefore, tests were 
carried out using glass-forming powders; the results 
are shown in Fig. 9. 

Sodium silicate glass applied to the stainless steel 
surface gave a friction-temperature plot similar to that 
of most of the low-melting metals, namely it showed the 
peak near the softening point (but much more pro- 


nounced), and the decrease above that temperature. 
However, these variations in friction disappeared com- 
pletely when sodium silicate glass was applied to 
a compacted boron nitride surface. 

Boron nitride was used since it is a material which is 
known to have poor wettability with respect to almost 
all liquids, but which can be wetted by boric oxide. 
Accordingly, the friction plot for boron nitride covered 
by boric oxide was determined, and this system does 
indeed show the characteristic peak and dip in the 
friction. 

Thus, there is evidence that wettability plays a very 
important role in the frictional behavior, not only of 
metals, but also of glass-forming materials. 


MATERIAL TRANSFER 


So far, the frictional properties of the various sliding 
systems have been related to the wettability of the 
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lubricant, which is a property that can be determined 
only above the melting point. Investigation of the sliding 
systems below this temperature has also been attempted, 
with interesting results. 

Friction tests were made using Wood’s metal riders 
on copper, stainless steel, phenol formaldehyde resin, 
and nylon surfaces. The temperature was kept constant 
at 60C during each test, and the duration of each 
test was 20 minutes. Figure 10 shows the photographs of 
the surface taken after the tests. 

It is seen that a considerable amount of material 
transfer of Wood’s metal has taken place on the copper 
and stainless steel surfaces, and some transfer to the 
phenol formaldehyde resin. The nylon surface does not 
show any trace of material transfer. These results are 
in striking similarity to the friction results of Fig. 7, 
where the friction curves of copper and_ stainless 
steel show a prominent peak at 60C, while phenol 
formaldehyde resin shows a less pronounced peak, and 
the curve for nylon does not show a maximum. 


The first feature of these data which warrants com- 
ment is the fact that, in general, soft metals such as lead 
and cadmium, in solid form, are found to be poor 
lubricants for steel, while these materials have generally 
been regarded as good lubricants. This distinction can 
be explained by the fact that other investigators (11, 12) 
used a thin film of the soft metal on the hard sliding 
surface, whereas this was not the case in the present 
experiments. It is known that the low friction of a soft 
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Fic. 9. Friction-temperature plot for glass-forming materials. Load, 1700 gm; speed 
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metal film is due to the ease with which the soft metal 
may be sheared, while the load is supported by the 
hard substrate. In the present tests the soft metal film 
tended to be much thicker, so that this lubricating effect 
disappeared. 

As to the results themselves, it has been shown that 
friction data obtained with low-melting metals at elevated 
temperatures are affected by oxide formation. Among 
the metals tested lead was most readily oxidized and 
this is in agreement with the oxidation data given in 
the literature (13). It is shown that the parabolic 
constant of the oxidation rate of lead is of the order 
of 10-1! gm?/cm* sec., while that of zinc and cadmium 
is of the order of 10—1* gm?/cm* sec. 

When tests are carried out in an argon atmosphere 
so that the friction test is but little affected by oxide 
formation, the friction plots fall into a general pattern, 
in which there is an increase in friction near the melting 
point and a subsequent decrease above that temperature. 
Only surfaces of low wettability, such as nylon (with 
Wood’s metal), or boron nitride (with silicate glass), 
show neither the peak nor the decrease. 

A likely explanation for the friction peak below 
the melting point is that, when the lubricating material 
becomes softer, adhesive effects between the material 
adhering to each surface become more pronounced, and 
consequently the area that has to be sheared, but it not 
supporting the load, becomes larger. It is clear that 
such an increase in shear area will result in an increase 
in friction. 

According to recent discussions of the adhesion 
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woop'S METAL ON 
COPPER 60C 


WOODS METAL ON 
PHENOL FORMALDEHYDE 


RESIN 60C 


WOOD'S METAL ON 
NYLON 60 C 


Fic. 10. Surfaces of various materials slid against solid Wood’s metal at 60 C 


process, a high level of adhesion would be expected to 
result from a high value of the ratio W,,/p, where W., 
is the surface energy of adhesion of the materials and 
p the hardness (14). For materials near their melting 
point, W,,/p is high (since p is low); consequently, 
strong adhesive effects are to be expected. 

Tests have been carried out to show the increase, 
with temperature, of the amount of cadmium adhering 
to a steel surface when cadmium riders were slid over 
a steel surface for one hour at constant temperature. 
Results are shown in Fig. 11. This system is rather 
similar to that of steel sliding on steel in the presence 
of cadmium powder, and the increase, with temperature, 
in the mass of adherent cadmium particles seems to be 
related to the increase, with temperature, of the W/p 
ratio. 

As regards the decrease in friction above the melting 
point, it is clear that this is due to the formation of a 


lubricating film of the molten metals. The sliding velocity 
in our tests was 0.1cm per second, and the viscosity 
of molten metals is about 2-5 centipoise, so that it is 
80 
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Fic. 11. Material transfer and friction for cadmium on stain- 


less steel in air. Load, 1700 gm; speed, 0.1 cm per second, duration 
of test 1 hour. 


a 
t 
20 3 
| 


most unlikely that a hydrodynamic lubrication film 
could have formed. Also, the results using Wood’s metal, 
which have shown the dependence of the friction on the 
nature of the bulk sliding surfaces, are consistent with 
boundary lubrication rather than hydrodynamic lubrica- 
tion. 

An adsorbed film of a molten metal is only a few Ang- 
stroms thick, compared to values of about 30 A for long- 
chain organic materials; hence it is natural that molten 
metals are much poorer lubricants than are long-chain 
organic liquids. The effectiveness of the molten metals is 
probably due to the fact that, when adsorbing on a solid 
surface, they lower its energy of adhesion. This explana- 
tion is borne out by the fact that those surfaces on which 
the molten material did not adsorb effectively, as indi- 
cated by lack of wetting ability, were those surfaces 
which the liquid did not lubricate effectively, even when 
the liquid was allowed to flood the surface. 

In fact, one of the most noteworthy results of these 
studies is that they demonstrate the importance of strong 
interaction between the interposed material and the slid- 
ing surface. Below the melting point, this interaction 
leads to the adhesion of large welded particles on the 
sliding surfaces, and consequently an increase in the 
friction, while above the melting point this interaction 
leads to the adhesion of liquid molecules on the solid 
surfaces, with the consequent occurrence of boundary 
lubrication. Neither the peak nor the drop in the friction 
was found without the other. 

As a final point, it is of interest to consider the rele- 
vance of these results to the problem of selecting good 
compatible materials for use under sliding conditions. 
One well-known criterion is that of miscibility (15). Of 
the four metals of Fig. 3, it is known that iron and 
cadmium are immiscible; that iron and lead as well as 
iron and zinc are miscible, while iron and tin form inter- 
metallic compounds. However, the frictional behavior 
was similar in all cases. 

The criterion of miscibility may be taken a step 
further. It does not seem at all likely that Wood’s metal 
and phenol formaldehyde resin are miscible. But yet this 
system gave a friction curve similar to that of miscible 
pairs. It would seem that we cannot predict the frictional 
behavior of a sliding system only on the basis of 
miscibility. 

However, the present tests suggest that wetting ability 
is a more reliable guide to the tendency towards the 
occurrence of high friction and high material adhesion 
in solids, and the tendency for lubrication to occur in 
systems with liquid lubricants present. Probably, the 
wettability which matters in this connection is not that 
measured in static tests, but rather the tendency for 
lubricant to wet the friction track during actual sliding. 
This can be considered a dynamic test. 

Bondi (16) suggests that wetting ability and misci- 
bility are highly correlated, so that a selection of material 
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pairs according to wettability might not differ greatly 
from one based on miscibility. However, even given this 
correlation, the criterion of wettability still seems to be 
a better one. During sliding, junctions are made and 
broken many times per second, so that there is little 
opportunity for the inter-diffusion of metals to occur. 
However, wettability indicates a high energy of adhesion, 
and it has been shown that this can influence friction and 
wear processes directly (14). 
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Antiwear and Extreme Pressure Additives 
for High-Temperature Lubricants 


By JOHN O. SMITH,' GLENN R. WILSON,” and JOHN R. STEMNISKI® 


New aircraft engines require new lubricants that will function under extreme conditions. Syn- 
thetic lubricants such as polyphenyl ethers are being explored as functiona! fluids in this new 
area. Although the thermal and oxidative properties of polyphenyl ethers are quite good, the 
lubricity properties are not outstanding. This paper summarizes research done in our laboratory 
on lubricity additives for polyphenyl ethers. A number of additives have been found to be effec- 
tive. In general they contained a hetero atom (phosphorus, halogen, sulfur) incorporated into an 
organic structure. The trichloroacetoxy functional group in an organometallic structure gave 
compounds with very high extreme pressure (EP) weld points. Although no pattern has been 
established for antiwear additives, several additives were found that reduced wear. These were not 
limited to one chemical class. Some additives gave good EP we'd points and also reduced wear, 
but no one compound was extremely effective in both areas. 


Introduction 


SEVERAL hundred compounds were screened as additives 
in a variety of lubricant base stocks (mineral oils, syn- 
thetic esters, organosilanes, silicones, and polyphenyl 
ethers) for improving their antiwear and extreme pres- 
sure (EP) properties at temperatures up to 600 F. The 
polyphenyl ether (PPE) fluids, of all those investigated, 
represent the most thermally (800-850 F) and oxida- 
tively (500-550 F) stable of any completely organic 
liquids developed to date. In this paper discussion is 
restricted to the PPE base stock. In particular, the PPE 
base stock used in all experiments was a blend of meta- 
and para-isomers of five-phenyl-containing polyphenyl 
ethers in proportions that provided the greatest liquid 
range. Although these PPE fluids have high thermal and 
oxidative stabilities, they are quite viscous (370 centi- 
stokes at 100 F and 13.2 centistokes at 210F), have 
poor low temperature properties (45 F pour point), 
average load-carrying capacity (140-150kg), and av- 
erage antiwear properties. 

For screening purposes, the conventional Shell four 
ball wear machine (modified to accommodate running 
tests up to 600 F) and Shell four ball EP machine were 
used. On several occasions, an especially good additive 
(as shown by Shell four ball tests) was submitted for 
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four-square small gear tests, but no correlation of the two 
tests was observed. 

The test procedures and screening results of selected 
additives are discussed in the following sections. 


Test apparatus and procedures 
SHELL Four BALL WEAR TEST 


The Shell four ball wear apparatus was modified to 
permit heating the base pot to 600 F. Each additive was 
screened at three temperatures, 167, 400, and 600 F at a 
constant load of 40 kg for 1 hour at speeds of 1265 rpm 
and at 1740rpm. Two types of test balls were used: 
52100 steel balls for the 167 F test and M-10 tool steel 
balls for the 400 and 600 F tests. 


SHELL Four BALL EP Test 


The conventional test, in which the load is increased 
until welding of the four test balls (52100 steel balls) 
occurred, was used throughout the screening tests. 


Discussion 


The results of Shell four ball EP tests, including 
structural formulas, reference number, concentration, 
average wear scar diameters at 167, 400, and 600 F, are 
summarized in Table 2. The numbers in parentheses in 
the last three columns of Table 2 are the wear scar 
diameters of the PPE base stock alone for base-line com- 
parison since variations resulted from changes in speed 
and lot numbers of test balls for a given series of tests. 
The boxing-in of areas in the last three columns of 
Table 2 is for the purpose of focusing attention on the 
best antiwear properties of a given additive for a given 
temperature. No single additive was significantly effec- 
tive at all three test temperatures. 
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Table | 


EFFECT OF ADDITIVES ON THE EXTREME PRESSURE 
PROPERTIES OF POLYPHENYL ETHER LUBRICANT 


No Additive WELD POINT (kg) 
None 140-150 
1 Cl, 1.0 1000 
Sn(occcl,) .0 
©). Ishe 750 
ci 
400 
C10 
v CH,COCH, 1.0 400 
cl 
cl 
vi 1.0 190 
OH 
vil 1.0 180-230 
cl 
cl cl 
vill oO 1. 6 130 
cl rf cl 
snc 0.5 160 
©) CVE 
0 
x (c1,cC0-Ti dy S 0.5 390 
xt +S-Ti-S4+, 1.0 330 
xu -cu 1.0 310 
xu ((n-C, Hy), -NCS) 1.0 220 
XIV | 0.5 180 
més} 
xv (ns Co 0.5 160 
2 2 
xvi NI 1.0 260 
xvil 1.0 250 


No Additive EB kg 
XVII (n-C, Hy), Sn 2.0 250 
H 
XIX 10 240 
xXx ©} P=s 4.0 240 
xXx! Sn ¢ ) 5.0 230 
©): )) 
1.0 230 
x Xt ©) Sb=S 10 210 
XXIV CH,), Lo 210 
bs be 
N (CH,), N(CH), 
(cH), 1.0 150 
XXVI 150 
CH 
(ch, -¢-0), POH 0.5 240 
1.0 270 
CH 
XXIX 10 260 
xxx 10 250 
XXX! 1.0 220 
XXXII 1.0 220 
xX Ni (NCS) 1.0 160 
XIV P| NiBr, 0.5 160 
XXXV (r-C,H, 0), 0.5 160 
0.5 150 


XXXVI —P—MdCO), 
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The Shell four ball EP and wear test results are dis- 
cussed separately in the following paragraphs since it is 
quite apparent from these data that good EP and anti- 
wear properties are not necessarily related. 


EXTREME PRESSURE ADDITIVES 


Table 1 comprises a list of 36 additives that were 
selected for illustrative purposes from over 100 com- 
pounds actually screened. A glance at the chemical struc- 
tures of this list of compounds shows that the major 
portion is limited to halogen-, sulfur-, and phosphorus- 
containing derivatives, many of which contain multiple 
heteroatoms so that any exact classification as to type is 
quite difficult. However, the additives are delineated 
roughly according to the predominant heteroatom such 
as halogen-, sulfur-, and phosphorus-containing com- 
pounds. 

The halogen-containing additives represent chemical 
structures in which the lability and bonding of the halo- 
gen varies from highly covalent to partially ionic. The 
best load carrying capacities were observed with com- 
pounds containing the trichloroacetoxy functional group. 
Thus, trichloroacetic acid (1) gave a weld point greater 
than 1000 kg, and, in fact, no welding occurred up to the 
load limit of the testing apparatus. Unfortunately, tri- 
chloroacetic acid is a very strong acid and is highly 
corrosive to steel. In order to remove this objectionable 
property, the acidic nature of the trichloroacetic acid 
was eliminated by converting it to an ester or salt. This 
resulted in some sacrifice of EP properties as illustrated 
by diphenyltin bis(trichloroacetate) (II), with a weld 
point of 750kg. 

The next highest EP weld points (490 and 400 kg) 
were obtained with halogenated derivatives of p-benzo- 
quinone, (III) and (IV), respectively. Additives (V) 
through (VIII) represent halogenated compounds in 
which the halogens are bound directly to aromatic nuclei, 
and these show very modest EP properties. Diphenyltin 
dichloride-bipyridyl complex (IX) represents an example 
of an ionic halogen-containing compound. It was not 
found to improve the EP weld point (160 kg with addi- 
tive) of the PPE base stock (140-150 kg for base stock). 

Additives (X) through (XXVI) in Table 1 represent 
a series of sulfur-containing compounds. The most effec- 
tive example in this series was the trichloroacetate deriva- 
tive of bis(cyclopentadienyl) titanium sulfide (X) with 
a weld point of 390 kg. The next highest weld point in 
the sulfur-containing compounds was obtained with a 
homolog of (X), the polysulfide of bis(cyclopentadieny] ) 
titanium, (XI). The effect of the trichloroacetoxy group 
in the former was apparent. Additives (XII) through 
(XV) are examples of metal salts of N,N-disubstituted 
dithiocarbamic acid that also show antioxidant activities 
in PPE base stock. The only derivative of these salts for 
which significant EP properties were found was nickel 
N,N-di-n-butyl-dithiocarbamate with a weld point of 
310 kg. One metal xanthate, nickel 2-ethylhexylxanthate 
(XVI), was also included in this series of sulfur-con- 


taining compounds, but its EP property was only mod- 
erate (260 kg). The remaining examples, (XVI) through 
(XXVI), also showed only moderate EP properties. 

A third group of additives represent a selection of 
phosphorus-containing compounds that includes phos- 
phites, dithiophosphate esters, phosphines, and metal 
salts. Diisopropyl hydrogen phosphite (XXVII), at 5% 
concentration, gave the highest weld point (400 kg) of 
the phosphorus derivatives. At 0.5% concentration of 
this same additive, the weld point dropped to 240 kg. 
Like trichloroacetic acid, diisopropyl hydrogen phosphite 
is highly corrosive to steel. The dithiophosphate deriva- 
tives, (XXVIII) through (XXXI), showed only mod- 
erate EP properties (270kg). The triphenylphosphine 
nickel compounds, (XX XIII) and (XXXIV), as well as 
the phosphate compound, (XXXV) and the phosphine 
molybdenum pentacarbonyl compound, (XXXVI), 
showed no improvement in the EP properties of the PPE 
base stock. 


ANTIWEAR ADDITIVES 


The list of antiwear additives tested, summarized in 
Table 3, includes many of the good EP compounds. 
Temperature had a profound effect on the antiwear prop- 
erties of a given additive and no one additive was found 
to exert its beneficial effects at all three temperatures. 
However, several were effective at two temperatures. 

From Table 2 it is apparent that a very limited number 
of compounds were effective at 167 F. Wear scar di- 
ameters of less than 1 mm were obtained with only four 
additives: diphenyltin bis(trichloroacetate), (II), copper 
N,N-di[butyldithiocarbamate] (XII), (both of which 
were also good EP additives), O-O-di-n-butylmorpholino 
phosphate, (XXXV) (which was a relatively poor EP 
additive), and diisopropyl hydrogen phosphite, (XX VII) 
(a relatively good EP additive). 

The only effective additive at both 400 and 167 F was 
diisopropyl hydrogen phosphite (XXVII). Sulfur-con- 
taining compounds represent the largest group of addi- 
tives effective at 400 F. The most effective of the sulfur 
compounds were bis(cyclopentadienyl) titanium poly- 
sulfide (XI), triphenylantimony sulfide (XXIII), 
dibutyltin bis(mercaptobenzothiazole) (XVIII), and 
triphenylphosphine sulfide (XX). The best halogen- 
containing additives effective at 400 F were hexa- 
chlorobenzene (VIII), bis(cyclopentadienyl) titanium 
bis(trichloroacetate) (XX XVIII), and dibutyl bis(pen- 
tachlorothiophenoxy) tin (XLIII). 

At 600 F three of the sulfur-compounds, triphenyl- 
antimony sulfide (XXIII), dibutyltin bis (mercaptobenzo- 
thiazole) (XVIII), triphenylphosphine sulfide (XX), and 
dibutyl bis(pentachlorothiophenoxy) tin (XLIII) were 
the only effective additives that were also effective at 
400 F. The other additives that were effective at 600 F 
represent a distribution of phosphorus, sulfur, and halo- 
gen compounds. The best phosphorus and sulfur com- 
pounds were  methylphenylthiophosphine chloride 
(XXVIII), and nickel O,O-di-2-ethylhexyldithiophos- 
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TABLE 3 
Additives Used in the Discussion 
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Trichloroacetic acid 


II Diphenyltin bis(trichloroacetate) 

Ill 2,5-Dianilino-3,6-dichloro-p-benzoquinone 

IV 2,5-Di-(5-nitrophthamido) -3 ,6-dichloro-p-benzoquinone 
V 
VI Pentachlorophenol 

VII 2,4,6-Triiodophenol 

VIII Hexachlorobenzene 

IX Diphenyltin dichloride bipyridy] complex 

x Bis[ (dicyclopentadienyl) titanium bis(trichloroacetate) ] sulfide 
XI Bis(cyclopentadieny]) titanium sulfide 

XII Copper N,N-dibutyldithiocarbamate 

XIII Nickel N,N-dibutyldithiocarbamate 

XIV Copper N,N-diphenyldithiocarbamate 

XV Cobalt N,N-diphenyldithiocarbamate 

XVI Nickel 2-ethylhexylxanthate 

XVII Diphenylthiocarbazone 

XVIII Dibutyltin bis(mercaptobenzothiazole) 

XIX 3-Benzylideneaminorhodanine 

xX Triphenylphosphine sulfide 

XXI Diphenyltin dithiophenol 

XXII m-Thioxylyl polymer 

XXIII Triphenylantimony sulfide 

XXIV bis(dimethyldithiocarbamate) 
XXV N,N-Dimethyl-2-benzothiazolesulfenamide 

XXVI Dibenzothiophene 

XXVII Diisopropyl hydrogen phosphite 

XXVIII Methylphenylchlorophosphine sulfide 

XXIX Bis(methylphenylphosphine sulfide) 

XXX Nickel bis(2-ethylhexyldithiophosphate ) 

XXXI Triphenylstanny]-di-n-decyldithiophosphate 

XXXII Benzenephosphonous acid 

XXXIII Bis(triphenylphosphine) nickel dithiocyanate 

XXXIV ___ Bis(triphenylphosphine) nickel dibromide 

XXXV O,0-Di-n-butylmorpholinophosphate 

XXXVI Pheny] dithiophene phosphine molybdenum pentacarbonyl 
XXXVII Copper (II) acetylacetonate 

XXXVIII Bis(cyclopentadieny]) titanium bis(trichloroacetate) 
XXXIX N,N-Diethyl-2-benzothiazole sulfinamide 

XL Hexahydro-1,4-diazepine-1,4-bis(2-methylenethiobenzothiazole) 
XLI Nickel O,O-di-2-ethylhexyldithiophosphate 

XLII Bis(pentachlorophenoxydiphenyltin) oxide 

XLII Bis(pentachlorothiopheny]) dibutyltin 

XLIV Bis(cyclopentadienyl) titanium thiocyanato chloride 


phate (XLI). The most outstanding sulfur compounds 
were N.N-diethyl-2-benzothiazole sulfinamide (XX XIX), 
and 3-benzylideneaminorhodanine (XIX). The best 
halogenated derivatives contained the pentachlorophenyl 
nucleus such as (XLII) and (XLIII). 


Summary 


From Shell four ball tests of a wide variety of com- 
pounds, a number of additives emerged that were mark- 
edly effective in improving the antiwear and EP prop- 
erties of PPE base stocks. The most exceptional types 
contain phosphorus, halogen, and sulfur or combinations 
of these heteroatoms. The only outstanding structural 
pattern that was evident from these data is the trichloro- 


acetoxy functional group which gave the highest EP weld 
points. It was also noted that no one compound is an 
effective antiwear additive at all three test temperatures 
(167, 400, and 600 F) although several were effective at 
two temperatures. 
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The Mechanism of Anti-Squawk Additive Behavior in 
Automatic Transmission Fluids 


By CLARENCE E. ALBERTSON! 


A theory of anti-squawk additive behavior was formulated based on diffusion and boundary 
lubricant formation. Friction tests were run with common friction surfaces in a common white 
paraffin oil containing different anti-squawk additives and analyzed to evaluate the theory. The 
data indicated dynamic friction was diffusion controlled at any given additive concentration. 
Additive concentration was not controlling dynamic friction at concentrations above 0.001 moles 
per liter, although the dynamic frictional behavior of zinc stearate solutions was concentration 
dependent below 0.001 moles per liter. Static friction data for zinc stearate followed a Langmuir 
adsorption isotherm. Friction changes caused by boundary lubricants were proportional to the 
amount of boundary lubricant on the friction surfaces. 


Nomenclature 
nm number of reaction rate controlling mol- 
ecules adsorbed on the friction surface 
per cm? 

nm, Maximum number of reaction rate con- 
trolling molecules which will adsorb or 

react with one cm? of friction surface 
n/n, _ traction of friction surface covered by 


adsorbed and reacted molecules (bound- 

ary lubricant) 
t time in seconds 
diffusion coefficient, cm? sec.—* 
thickness of immobile surface layer, cm 
c¢ concentration of reaction rate controlling 
molecules, moles per liter 
Avogadro’s number 
absolute temperature 
Fahrenheit temperature 
viscosity, poise 
diffusion constant for each reactant-oil 
solution, [3] 
average effective asperity distance, cm 

v _ sliding velocity, cm per second 

(Afa),,r change in coefficient of friction due to 
adsorbed boundary lubricant at tem- 
perature T and velocity v 
change in static coefficient of friction at 
temperature T due to maximum bound- 
ary lubricant adsorption 
coefficient of friction of the base oil at 
temperature 7 and velocity v 


Contributed by the ASLE Technical Committee on Lubrication 
Fundamentals and presented at the Annual Meeting of the 
American Society of Lubrication Engineers held in New York, 
April, 1963. 

1 Staff Scientist, Borg-Warner Research Center, Wolf and 
Algonquin Roads, Des Plaines, Illinois. 
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(Af.)r 


(fo) 


300 


(fora)v.r coefficient of friction of the additive and 

oil solution at velocity v and tempera- 
ture T 

d__ an assumed constant for a given sliding 
system, d = (N/b) /(1000 Ahn.,) 

J grouped exponential factors for a test 
condition, J = (cTd)/(u) 

K proportionality constant in n, = — 

x,0,A marks indicating static friction values 
at increasing temperatures, in the 
graphs, or static friction change 

Introduction 


Tuis paper describes the formulation and evaluation of a 
theory of anti-squawk additive behavior in automatic 
transmission fluids. Anti-squawk additives depress the 
static friction and tend to impart a positive slope to the 
change of friction with sliding speed, (df)/(dv) = +. 
when wet clutches are engaged. Several authors (1-3) 
have offered qualitative explanations of this effect based 
on the destruction and repair of adsorbed or chemically 
formed boundary lubricating films. As sliding speed de- 
creases, repair of the lubricating film exceeds destruction 
and friction decreases. Boundary lubricants such as fatty 
acids (1), soaps (4), dialkyl phosphates (3), and dialkyl 
phosphites (3) are anti-squawk agents. This information 
and other hypotheses were used to develop a more 
quantitative explanation of anti-squawk additive action. 

Surface reactions leading to the formation of boundary 
lubricating films were reviewed. Equations were derived 
which describe the amount of oxygen, water or additive 
which diffuses to the sliding friction surfaces. The amount 
of boundary lubricant on the surface was related to the 
change in coefficient of friction at different sliding speeds. 
Friction vs. sliding speed data were measured in different 
additive-oil solutions and analyzed to evaluate the theory. 
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Discussion 
BouNDARY LUBRICANT FORMATION 


Boundary lubricating films are formed on metal or 
oxidized metal surfaces either by adsorption of boundary 
lubricants or by the reaction of boundary lubricant form- 
ing materials with the surfaces (5). For example, adsorbed 
zinc laurate and other soaps were shown to be effective 
boundary lubricants on either reactive or inert metal 
surfaces. Desorption into oil reduced their effectiveness. 
The formation of boundary lubricating films on reactive 
metal surfaces can occur directly, as the case of steel and 
chlorine gas (6), or indirectly, as the formation of soaps 
by fatty acids. In this case, film formation occurs in 
several steps, oxidation of the metal and reaction of the 
metal oxide with acid to form the boundary lubricating 
soap (7, 8). Water accelerates metal oxidation and soap 
formation (5). 

The anti-squawk boundary lubricants selected for test- 
ing included an adsorbing type, zinc stearate, and several 
reactive types, a fatty acid, stearic acid, a reactive, acidic 
lauryl phosphate and dioctyl phosphite, which forms 
reactive mono-octyl acid phosphite in the presence of 
moisture (3). 

Anti-squawk additive behavior has been related to the 
boundary lubrication of sliding surfaces in automatic 
transmission clutches (1-3). It is proposed that the metal 
surface is oxidized by oxygen and water dissolved in the 
lubricant, and anti-squawk additives in the lubricant 
adsorb on or react with the metal and/or oxides to form 
the boundary lubricating films. The amount of film on the 
surface, either oxide and/or boundary lubricant, is 
dependent on the relative rates of formation and abrasion 
of the film in sliding. 

Boundary lubricant formation is rapid as evidenced by 
the anti-squawk additive behavior reported in this work 
and by others (2, 3). Data on zinc stearate, Fig. 1, and on 
stearic acid, Fig. 2, show they form boundary lubricating 
films rapidly. The ordinate, Afa, is the additive induced 
friction change. A degree of mechanical activation (9) 
may be contributing to boundary lubricant formation be- 
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Fic. 1. Frictional behavior of 0.0005 moles per liter zinc 
stearate in white oil. 
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Fic. 2. Analysis of 0.030 moles per liter stearic acid frictional 
behavior. 


cause the steel friction specimen, Fig. 3, is always 
abraided and polished in testing. Since boundary lubri- 
cant formation is rapid, it would be diffusion controlled 
and the reactant or adsorbate which diffuses more slowly 
would control the amount of lubricant on the friction 
surface. The relative rates of formation and desorption 
(10, 11) also influence the amount of boundary lubricant 
on the friction surfaces. 

The relative concentrations of oxygen, water and 
additive in the oil were examined. The solubility of air in 
petroleum oils is about 7—-11% by volume (12). Oxygen 
concentration would be about 2% or 0.0009 moles per 
liter. The white oil used to make the experimental lubri- 
cants contained 18 ppm water, by Carl Fisher analysis. 
This is approximately 0.0008 moles of water per liter of 
oil. The oxygen and water concentrations in the lubricants 
were about the same. Oxygen and water diffusion would 
be a little greater than that of larger additive molecules 
such as stearic acid. Additive concentrations below 0.001 
moles per liter would be expected to control boundary 
lubricant film formation. Oxygen and water concentration 
would be expected to limit boundary lubricant formation 
by additives more concentrated than 0.002 moles per liter. 
Friction tests were run at anti-squawk additive con- 
centrations below 0.001 moles per liter and above 0.002 
moles per liter to test these qualitative predictions. 


SURFACE DIFFUSION 


The diffusion rate of either oxygen, water or anti- 
squawk additive to a friction surface, in a well stirred 
lubricant, can be expressed by Eq. [1] (13), if the 
concentration is not changed appreciably by adsorption 


and reaction. 
dn DcN (1 n ) [1] 
dt Ah 1000 


Nz 
Equation [1] was integrated to determine the amount 
of oxygen, water or additive adsorbed or reacted in time ¢ 
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Fic. 3. 


and transposed to derive an expression [2] for the frac- 
tion of the available surface covered by boundary lubri- 
cant. Boundary lubricant formation is fast, as evidenced 
by the experimental data. 


n ( —cDNt ) [2] 
— = l —e 
nN, my, 1000 Akn,, 


The diffusion coefficient is known to vary directly with 
temperature and inversely with viscosity (14). 


bT 
D = — [3| 
Substituting Eq. [3] in Eq. |2] one obtains Eq. |4| 
n —cNTtb [4] 
— = 1 — exp| ————_- 
Nn, 1000 Akn,u 


This equation relates the basic factors which influence 
the amount of boundary lubricant on a friction surface. 
It is assumed that material which diffuses to the surface 
is adsorbed and/or reacted. This relation holds for either 
oxygen-, water- or additive-controlled boundary lubricant 
formation. 

FRICTION THEORY 


Paper-based automatic transmission clutch linings can 
dissipate oil films and run to sliding speeds above 75 cm 
per second in boundary lubrication. Figure 4 depicts the 


Friction specimens, tumbled 1035 steel, and paper-based lining PB-1 


0.35 


s00F 


x 


20. 40 60 
SLIDING SPEED,v, (cm/sec) 
Fic. 4. Frictional behavior of white, refined, paraffin oil 


COEFFICIENT OF FRICTION 


80 


friction vs. sliding speed behavior, at different tempera- 
tures, of a paper-based lining PB-1 run against a tumbled 
1035 steel disc in a refined, white, paraffin oil. The co- 
efficient of friction at all temperatures is much higher 
than hydrodynamic friction, and friction does not rise at 
increasing sliding speeds. The porous, paper matrix is a 
blotterlike reservoir for transmission fluid. Additives, 
water, and oxygen in the fluid are able to diffuse to the 
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separator plate. The paper fibers are so small and rela- 
tively widely spaced that they are generally not sup- 
ported by a hydrodynamic fluid film and rub off part or 
all of the oxide and or boundary lubricant molecules as 
they slide along. Paper-based linings are a unique 
friction material for evaluating boundary friction be- 
havior. 

If the lining slides faster, less time is available between 
fiber contacts for oxygen or additive adsorption and fewer 
boundary lubricant molecules will be formed at the sur- 
face to lubricate the next fiber. Adsorption time is 
related to sliding speed and average effective asperity 
distance by | 5]. 


v= — [5] 


The frictional effect of the boundary lubricant was 
formulated by Eq. [6], assuming boundary friction 
decreases in proportion to the amount of adsorbed 
boundary lubricant. Previous work indicated this assump- 
tion was valid under equilibrium conditions (10). A fully 
adsorbed surface at infinite time would have the lowest 
friction. Surfaces sliding at infinite rubbing speed, or zero 
adsorption time between fiber contacts, would have the 
coefficient of friction of the base oil. 


(Afa) = (Af.)r(n/nx ) | 6] 


Substituting [4] and [5| in |6] yields [7], the ex- 
pression for the effect of diffusion on friction: 


E ( —cNTlb )] 
(Meer = (Afe)e| 1—emp\ 
[7] 


For purposes of calculation: 
(Afa)..7 = (fo + a) (fo) |8| 


Assuming d is a constant for a given system, |7| can be 
simplified to [9]. 


—cTd 
(Afa)..7 = 1 )] [9] 


For curve fitting over a range of sliding speeds [9] was 
simplified to [10]. 
(Afa),,r = (Afz)r|1—exp (—J/v)] [10] 


An equation can be derived for calculating the coefficient 
of friction of an additive fluid by transposing [8] and 
|9]: 


—cTd 
(fo + = (fo). + — exp( )] 


[11] 


Equations |7] and [11] relate the factors influencing 
the change of dynamic friction with sliding speed due to a 
boundary lubricant. They do not predict the coefficient of 
friction measured in pure oil, fo, the static coefficient 


depression due to the boundary lubricant, Af,, or the 
changes of these basic factors with temperature. They 
correlate the effects of temperature, fluid viscosity and 
sliding speed on dynamic friction. They give an idea of 
the effects of surface asperity configuration and diffusion 
parameters on dynamic friction. : 


Test PROCEDURE 


Friction vs. sliding speed data were taken with a low 
velocity friction apparatus (LVFA), Fig. 5 (2, 3). The 
specimens were immersed in 15 ml of lubricant at 127 psi 
load. They were run-in overnight at room temperature 
and 20.3 cm per second sliding speed. In friction vs. 
sliding speed tests, temperatures were held within + 5 F 
and sliding speed was varied uniformly with a variable- 
speed drive, from 0 to 76.2 cm per second to O in 2 
minutes. Lining surface temperatures were about 5 F 
above holding temperatures at higher sliding speeds. Each 
friction trace was repeated after holding at temperature 
for about 5 minutes at 25.4 cm per second slip speed to 
be sure the trace was stable. Temperature conditions 
during a test are discussed in Appendix 1. 

The friction specimens, Fig. 3, were paper-based linings 
PB-1, cut into %-inch wide annuli, 1-inch mean diameter, 
and driven discs of tumble finished 1035 steel. The same 
friction surfaces were used with the same additive, with 
progressively higher additive content lubricants. When 
increasing additive concentration, several changes and 
runs in fresh lubricant were used to work out the 
previous lubricant. 


The frictional stability of the friction materials was 
checked following additive runs by cleaning them 
thoroughly and repeating the white oil tests. The linings 
were cleaned until the original static friction values were 
duplicated because static friction is the most sensitive 
indicator of additive contamination. Static friction is 
relatively insensitive to changes in lining structure, 
whereas dynamic friction is structure sensitive. Friction 
surfaces were generally stable within the experimental 
variation. Tests with zinc stearate showed the largest 
deviations. During these tests the dynamic friction level 
at all temperatures climbed from 0.01 at 76 cm per second 
to a maximum of 0.03 at 25.4 cm per second. These 
white oil data, measured before and after the additive 
runs, were averaged for the best (fo),.7 data. Testing 
was limited to 300F to reduce side effects due to 
lubricant degradation. Friction traces were reproducible 
and the fluid was still water-white when most of the runs 
were terminated. 


The dynamic friction data were usually reproducible 
within + 0.01 and static within + 0.02 coefficient of 
friction for different specimens of the same materials. 


Experimental data 


The friction vs. sliding speed behavior of common 
paper-based lining PB-1 and 1035 steel specimens was 
measured in a common white paraffin oil containing 
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Fic. 5. Low velocity friction apparatus 


various anti-squawk additives. The white, refined par- 
affin oil data are depicted in Fig. 4. Static friction values 
at increasing temperatures were indicated at the ordinate 
by an x, 0, A. The dashed portions of the curves were ex- 
trapolated because the machine was chattering. Repre- 
sentative data for a white oil solution containing 0.0268 
moles per liter dioctyl phosphite is given in Fig. 6. 

Data for other additive-containing lubricants were 
plotted as Afa vs. sliding speed, in Figs. 6-13. The addi- 
tive-induced friction changes (Afa),,7 were calculated 
[8] from friction vs. sliding speed data. A negative Afa 
indicates the friction of the same lining specimens 
measured in an additive-white oil solution is lower 
than that measured in pure white oil, Afa = fa. 
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Fic. 6. Frictional behavior of white oil containing 0.0268 moles 
per liter dioctyl phosphite. 


COEFFICIENT OF FRICTION 


— f,. The addition of the additive has depressed 
the friction level of the system. Additives which form 
boundary lubricants on the friction surfaces would cause 
Afa to be negative, especially at low sliding speeds. Con- 
versely, a positive Afa indicates the additive has raised 
the friction level of the system above that measured in 
pure white oil. An additive which reduces the viscosity 
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Fic. 7. Frictional behavior of low concentrations of zinc 
stearate in white oil. 
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TABLE 1 
Lubricant Viscosity-Density Data 
100 F 200 F 300 F 
Lubricant Cpe Density? cP Density CP Density 

White oil containing: 
Pure oil 27.3 0.849 4.95 0.812 1.91 0.780 
0.0020 moles per liter 

zinc stearate 26.4 0.850 4.58 0.812 1.86 0.780 
0.0268 moles per liter 

dioctyl phosphite 26.7 0.850 4.94 0.815 1.93 0.778 
0.0803 moles per liter 

dioctyl phosphite 24.6 0.850 4.45 0.819 1.81 0.783 
0.0200 moles per liter 

stearic acid 25.7 0.852 4.62 0.815 1.80 0.780 


@ CP = viscosity (cp). 
» In gm/cc. 


of the white oil could cause Afa to be positive at higher 
sliding speeds, where boundary lubricants do not have 
time to form. 

The mixed lauryl phosphate additive was a mixture of 
mono and di lauryl phosphates. Specimens in zinc stearate 
lubricants were run at 200F and cooled prior to the 
regular test series to activate the surfaces, because in- 
duction periods in boundary lubricant formation have 
been noted. Evidently, some metal specimens do not 
adsorb zinc stearate until their original surface has been 
worn away. The zinc stearate solution was clear to 
0.00005 moles per liter. Zinc stearate and stearic acid 
solutions of higher concentration were cloudy at 100 F 
and clear at 200 F. All other additives were completely 
soluble. 

Lubricant viscosity and density data in Table 1 show 
the decrease in viscosity at increasing temperatures and 
the viscosity changes caused by the additives. 


Analysis 
BEHAVIOR OF ZINC STEARATE AT LOW CONCENTRATION 


Zinc stearate data at 0.00005 moles per liter (Fig. 7; 
Table 2) indicate static friction is controlled by adsorp- 
TABLE 2 


Afa Values for Static Friction of Zinc Stearate Solutions, 
from Figs. 1,7, and 8 


— Afa values at different temperatures 
(moles/liter) 100 F 200 F 300 F 
0.00005 —0.15 —0.09 0 
0.0001 6 —0.115 —0.09 
0.0005 —0.31 —0.25 
0.0020 bad —0.30 —0.27 


@ Zinc stearate was not completely soluble at this temperature 
and concentration. 


tion at this concentration. The static friction change de- 
creased significantly as temperatures increased, indicating 
the amount of boundary lubricant on the friction surfaces 


was decreasing. The zinc stearate was evidently desorbed 
at 300 F. The 300 F dynamic curve has a peculiar shape, 
however it was within + 0.01 of the zero line to a sliding 
speed of 30 cm per second and within + 0.02 at higher 
speeds. The maximum variation of the Afa data would be 
+ 0.02, the sum of the variation of the dynamic friction 
curves used to calculate Afa. Zinc stearate data above 
0.00005 moles per liter, in Table 2, taken from Figs. 1, 
7, and 8, are similar to the more dilute data but are not 
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Fic. 8. Analysis of 0.002015 moles per liter zinc stearate fric- 
tional behavior. 


definitive because the static friction differences at in- 
creasing temperatures usually fall within the maximum 
frictional variation. The friction changes at 100F for 
solutions above 0.00005 moles per liter are less than 
they should be because the zinc stearate was not com- 
pletely soluble at these temperatures. 

Assuming Af, was proportional to the fraction of the 
surface covered by a boundary lubricant, Miller and 
Anderson (10) modified the Langmuir adsorption iso- 
therm to relate the change in friction with additive con- 
centration. The friction and additive concentration data 
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were shown to follow a straight line when c/Af, was 
plotted against c. 

Static friction data for zinc stearate were analyzed to 
determine whether they follow a Langmuir adsorption 
isotherm. The c/Af, values are given in Table 1A (Ap- 
pendix 2) and graphed in Fig. 9. The data followed 
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Fic. 9. Analysis of isothermal zinc stearate adsorption data 


straight lines, within the experimental variation. This 
straight line fit indicates the amount of boundary lubri- 
cant adsorbed under equilibrium conditions is dependent 
on the zinc stearate concentration. It does not differen- 
tiate whether the boundary lubricant is zinc stearate or 
a zinc stearate product. 


METHOD OF ANALYZING DYNAMIC FRICTION DATA 


The diffusion equations derived earlier were used to 
determine whether dynamic frictional behavior is diffusion 
controlled. The exponential in [7], (cNTIb)/(1000 
Ahn,,1v), contains the fluid and system variables which 
influence diffusion-controlled dynamic frictional behavior. 
The / term is the average effective asperity distance, i.e., 
the average sliding distance required to rub off an ad- 
sorbed boundary lubricant molecule. It may be longer 
than the distance between asperities if several wipes are 
required. The n,, term is the maximum number of bound- 
ary lubricant molecules which can adsorb on the friction 
surface. The term increases as a function of concentration 
and is proportional to the decrease in static friction, 
n, — — KAf,. Of these factors /, b, Ahk and n,, have not 
been measured directly. The other factors c, T, wu, and v 
have been determined. The unevaluated terms, /, 5, n., 
and AA were lumped into factor, d, and determined by 
fitting [9] and [10] to one of the experimental curves. 
Using this empirical d factor and measured changes in 
c, T, wu, v and Af,, the related curves at other concentra- 
tions and temperatures were calculated and compared 
with the related experimental curves to determine whether 
the diffusion equation would approximate the experi- 
mental friction vs. sliding speed data. 

The 200 F dynamic friction curves for low concentra- 
tions of zinc stearate were grouped for comparison in 
Fig. 10. Equation [10] fitted to the 0.0005 moles per 
liter curve took the form: 


Afa =—o315|1 exp(— | [12] 
v 


and yielded the dashed curve drawn in Fig. 10. Curves 
were calculated for the 0.00005 and 0.0001 moles per 
liter data, Table 2A, using c = additive concentration, 
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Fic. 10. Analysis of the frictional behavior at 200F of low 
concentrations of zinc stearate. ; 


n, = — KAf, and a common constant for grouped 
factors |13|. The calculated values plotted as broken 
lines in Fig. 10 reproduce the laboratory data within ex- 


c-—3390 


perimental variation. The relative fit of the calculated 
and experimental curves at different additive concentra- 
tions is evidence for the control of boundary lubricating 
film formation by diffusion at additive concentrations of 
0.0005 moles per liter and less. 


FRICTIONAL BEHAVIOR AT HIGHER ADDITIVE 
CONCENTRATIONS 


Data obtained at higher additive concentrations were 
examined. Figures 1, 8, 11, 12, 13, and 14 show that 
three- to fourfold increases in the concentration of zinc 
stearate, dioctyl phosphite, mixed lauryl phosphates, and 
stearic acid did not appreciably change dynamic friction 
vs. slip speed behavior. These tests indicate the dynamic 
friction change, Afa, was not significantly influenced by 
additive concentration changes, at concentrations above 
0.0005 moles per liter for zinc stearate and above 0.002 
moles per liter for the other additives. 

An intermediate factor appears to be limiting the rate 
of boundary lubricant formation by the diffusing addi- 
tives. The limiting step may be the rate of surface oxide 
formation by oxygen and water dissolved in the lubricant 
(5, 7). 

Exploratory tests to detect oxygen and moisture effects, 
by this laboratory and others (3), were not successful 
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because the atmosphere cannot be controlled well enough 
in the LVFA. Tests run in nitrogen filled transmissions 
(3) gave evidence of decreased lubricity, indicating air 
was necessary for anti-squawk additives to be effective. 
Although circumstantial evidence points to an_inter- 
mediate oxide step in boundary lubricant formation, no 
concrete evidence for this effect was found. 

The static friction data at 200 F and 300F, for the 
more acidic additives, stearic acid and mixed lauryl 
phosphate, and data at 300F for 0.08 moles per liter 
dioctyl phosphite, were abnormally high compared to 
dynamic friction behavior (see Figs. 11, 13, and 14). 
These high static values were time dependent. Dynamic 
friction decreased with sliding speed to $ to 1 cm per 
second, indicating the boundary lubricant film was 
forming normally to nearly zero sliding speed. This 
behavior appears to be caused by a lining-additive in- 
teraction. These static friction data on the more reactive 
additives are not stable enough for interpretation. 

The magnitude of the static friction change and the 
lubricity of the boundary lubricant film are determined by 
the type of additive. Stearic acid did not depress friction 
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Fic. 11. Analysis of 0.0803 moles per liter dioctyl phosphite 
frictional behavior. 
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Fic. 12. Analysis of 0.0268 moles per liter dioctyl phosphite 
irictiona! behavior. 


nearly as much as dioctyl phosphite and zinc stearate. 
Lubricity changes are evidenced by the variation in the 
slopes of the Afa vs. sliding speed data of the different 
additives. 

Dynamic friction data at a given additive concentration 
and different temperatures were analyzed to determine 
whether the influence of temperature and viscosity on 
diffusion would account for the shape and relative posi- 
tions of the experimental curves. Equation | 10] took the 
form |14] when fitted to the 200 F data for 0.08 moles 
per liter dioctyl phosphite in Fig. 11. The calculated Afa 
data are given in Table 3A and plotted as dashed lines 
with the experimental data in Fig. 11. The empirical cd 
factor for 0.08 moles 
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Fic. 13. Effect of mixed lauryl phosphate concentration on 
frictional behavior. 
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Fic. 14. Effect of stearic acid concentration on frictional be- 
havior. 
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per liter dioctyl phosphite is 0.000258 (see Table 3A). 
The Afa values at 100 and 300 F were calculated using 
[9], the empirical cd and viscosity data. The cd value 
calculated from J in [14] (see [10]) does not distinguish 
the diffusing species. These values in Table 4A were 
plotted as dashed lines with the experimental data in 
Fig. 11. The calculated curves fit the data within experi- 
mental variations, with the exception of the 0—7.62 cm 
per second portion of the 300 F curve. Other data were 
analyzed with similar results (see Fig. 12 and Table 5A 
for 0.0268 moles per liter dioctyl phosphite; Fig. 8 and 
Table 6A for 0.002 moles per liter zinc stearate, and 
Fig. 2 and Table 7A for 0.03 moles per liter stearic 
acid). The general curve fit, curve spacings at different 
temperatures and similar behavior of the different anti- 
squawk additives is evidence that dynamic anti-squawk 
additive frictional behavior is controlled by diffusion. 
The poor fit of the 300 F curves at 0-7.62 cm per 
second sliding speed is common. It appears to be due 
either to an increase in effective asperity distance /, or to 
a change in Ah, or to both, as the degree of boundary 
lubrication decreases with increasing sliding speed (note 
Fig. 4). 


Summary 


Simplified equations were derived relating dynamic 
friction to the diffusion of a boundary lubricant or 
lubricant-forming species. Using these equations, the 
dynamic frictional behavior of zinc stearate solutions was 
shown to be concentration dependent up to 0.0005 moles 
per liter. The static friction of zinc stearate solutions 
followed Langmuir adsorption isotherms, indicating this 
additive acts by adsorption. Dynamic frictional behavior 
was not concentration dependent at additive concen- 
trations above 0.001 moles per liter for zinc stearate and 
the other additives tested. No concrete evidence was found 
to explain this behavior, although circumstantial evidence 
indicated an intermediate step, such as oxide formation, 
is limiting the additive effect. The diffusion equations 
correlated the effects of temperature, viscosity, and 
sliding speed on the dynamic frictional behavior of these 
more concentrated additive solutions at any particular 
concentration. Agreement of experimental and calculated 
data indicated the additive-induced friction change, Afa, 
is proportional to the amount of boundary lubricant on 
the sliding surfaces. 
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Appendix 1 
LINING SURFACE TEMPERATURES 


In six runs, the average increase in lubricant tem- 
perature per run was 3.3 F at the center of the test 
specimens. This was the monitoring thermocouple. A 
thermocouple positioned through the lining to the lining 
surface, on the mean radius, showed an average increase 
in temperature of 8 F at the higher slip speeds, dropping 
to 3.8 F above the starting temperature at the end of the 
runs. The maximum difference between the center monitor 
thermocouple and the lining surface thermocouple was 
about 5 F. At the beginning of a run the difference was 
zero and at the end of the run the difference was 0.5 F. 

The linings have a high proportion of interconnected 
porosity. The lubricant carries surface heat rapidly into 
the lining and through the lining. The surface thermo- 
couple, bathed in lubricant, should be close to the actual 
surface temperature. 

The maximum lining surface temperature would be 
+ 8F above the starting temperature, in the higher 
sliding speed range of the test. In this range Afa values 
are a minimum. The + 8F surface temperature would 
increase Afa approximately 0.001 at the higher sliding 
speeds, which is well within the experimental variation. 
It is felt that the lining surface temperature increase does. 
not significantly influence the test data and their analysis. 
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Appendix 2 
TaBLe 1A TABLE 2A 
Zinc Stearate Adsorption Data Calculated Friction Changes Caused by Zinc Stearate at 200 F 
Data graphed in Fig. 9. Factors used in Calculation: 366 T, absolute; ~ = 0.0495 poise; 
= a 
Temperature (F) ct Afa® (c)/(Afa + 0.04) (NTIb) /(—K 1000 = — 3390°. 
200 ‘a 5 0.09 0.00038-0.0010 Zinc stearate, moles per liter 0.00005 0.0001 0.0005 
Af,, see nomenclature —0.09 —0.115 —0.315 
0.0003 0.29 0.0009 1-0.0012 
0.0005 0.315 0.0014 -0.0018 Sliding-speed (v) 
0.0010 0.32 0.0028 0.0036 (cm/sec) 
0.0015 0.325 0.0041 -0.0053 
0.0020 0.325 0.0055 -0.0070 2.54 
300 0.00005 0 -- 5.08 —0.0278 —0.051 —0.205 
0.0001 0.10 0.0007 2-0.00167 7.62 —0.0197 —0.037 —0.160 
0.0003 0.20 0.00125-0.0019 12.7 —0.0124 —0.024 —0.109 
0.0005 0.265 0.0016 ~+:.0022 25.4 —0.0064 —0.0126 —0.060 
0.0010 0.30 0.0029 -0.0038 50.8 —0.0033 —0.0065 —0.0316 
0.0015 0.305 0.0044 -0.0057 76.2 —0.0022 —0.0043 —0.0214 
0.0020 0.305 0.0058 -0.0075 @ When ty = —Kaf,,. 
@ Concentration of zinc stearate, moles per liter. » J calculated from additive concentration, J = (cx — 3390)/ 
> Average of 3 different runs: at 200 F, the 3 runs were within (4f,,)- 
+ 0.03; at 300 F, the 3 runs were within + 0.04. ¢ Calculated by Eq. [13]. 
TABLE 3A 


Calculation of the Friction Change at 200 F Caused by 0.0803 Moles per Liter of Dioctyl Phosphite 


2.12 
(Afa) o99 » = —0.285 | 1 — exp| — —) ] [12] 
v 


Sliding speed 
ft/min cm/sec (v) J/v exp(—J/v) 1—exp(—J/v) —exp(—J/v)] 
0 0 1 —0.285 
§ 2.54 0.834 0.434 0.566 —0.161 
10 5.08 0.417 0.659 0.341 —0.0972 
15 7.62 0.278 0.757 0.242 —0.0690 
25 12.7 0.167 0.846 0.154 —0.0438 
50 25.4 0.083 0.920 0.080 —0.0228 
100 50.8 0.0417 0.959 0.041 —0.0117 
150 76.2 0.0278 0.973 0.0274 —0.0078 
Af, = —0.285 cd = (2.12 p)/(T) 
Jz 232 cd = (2.12 X 0.0445) /(366) 


cd = 0.000258 
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TaBLe 4A 
Calculated Friction Change Caused by 0.0803 Moles per Liter 
Dioctyl Phosphite in White Oil 


Factors used in calculation 


F, Fahrenheit 100 200 300 
T, absolute 311 366 422 
cd, see nomenclature 0.000258 0.000258 0.000258 
Ll, poise 0.246 0.0445 0.0181 
J, see nomenclature 0.326 2:42 6.01 
Af, see nomenclature —0.235 —0.285 —0.3104 
Sliding speed (v) 
(cm/sec) Calculated Afa values? 
0 —0.235 —0.285 —0.310 
2.54 —0.0283 —0.161 —0.281 
5.08 —0.0146 —0.097 —0.215 
7.62 —0.0098 —0.069 —0.169 
12.7 —0.0059 —0.0437 —0.117 
25.4 —0.0030 —0.0228 —0.065 
50.8 —0.0015 —0.0116 —0.034 
76.2 —0.0010 —0.0078 —0.023 


“ Af,, = —0.310 was extrapolated from Figs. 5 and 6 because 
of abnormally high static friction. 
» Calculated by Eqs. [9] and [10], from a common cd value. 


TABLE 5A 
Calculated Friction Change Caused by 0.0268 Moles per Liter Dioctyl Phosphite in White Oil 


Factors used in calculation 


F, Fahrenheit 100 200 300 
T, absolute 311 366 422 
cd, see nomenclature 0.000282 0.000282 0.000282 
Lt, poise 0.267 0.0494 0.0193 
J, see nomenclature 0.328 2.09 6.14 
Af,,, see nomenclature —0.235 —0.270 —0.295 
Sliding speed (v) 
(cm/sec) Calculated Afa values” 
0 —0.235 —0.270 —0.295 
2.54 —0.0285 —0.152 —0.269 
5.08 —0.0147 —0.0912 —0.207 
7.62 —0.0098 —0.0647 —0.163 
12.7 —0.0060 —0.0409 —0.113 
25.4 —0.0030 —0.0213 —0.0635 
50.8 —0.0015 —0.0109 —0.034 
76.2 —0.0010 —0.0073 —0.023 


“ Calculated by Eqs. [9] and [10], from a common cd value. 
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TABLE 6A 
Calculated Friction Change Caused by 0.002015 Moles per Liter Zinc Stearate in White Oil 


Factors used in calculation 


F, Fahrenheit 100 200 300 
T, absolute 311 366 422 
cd, see nomenclature 0.000350 0.000350 0.000350 
ll, poise 0.264 0.0458 0.0186 
J, see nomenclature 0.413 2.80 7.95 
Af, see nomenclature —0.200 —0.303 —0.270 
Sliding speed (v) 
(cm/sec) Calculated Afa values@ 
0 —0.200 —0.303 —0.270 
2.54 —0.0300 —0.206 —0.258 
5.08 —0.0136 —0.129 —0.213 
7.62 —0.0085 —0.093 —0.175 
12.7 —0.0064 —0.061 -—0.126 
25.4 —0.0032 —0.032 —0.072 
50.8 —0.0016 —0.016 —0.039 
76.2 —0.0011 —0.011 —0.026 


@ Calculated by Eqs. [9] and [10], from a common cd value. 


TABLE 7A 
Calculated Friction Change Caused by 0.0300 Moles per Liter Stearic Acid in White Oil 


Factors used in calculation 


F, Fahrenheit 100 200 300 
T, absolute 31L 366 422 
cd, see nomenclature 0.000618 0.000618 0.000618 
Ll, poise 0.257 0.0462 0.0180 
J, see nomenclature 0.75 4.50 14.05 
Af,, see nomenclature —0.206 —0.150 —0.190" 
Sliding speed (v) 
(cm/sec) Calculated Afa values? 
0 —0.206 —0.150 —0.190 
2.54 —0.0525 —0.1282 —0.189 
5.08 —0.0282 —0.0928 —0.178 
7.62 —0.0193 —0.0710 —0.160 
12.7 —0.0118 —0.0480 —0.127 
25.4 —0.0060 —0.0263 —0.081 
50.8 —0.0030 —0.0138 —0.046 
76.2 —0.0020 —0.0093 —0.032 


@ Extrapolated due to high static friction. 
» Calculated by Eqs. [9] and [10], from a common cd value. 
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DISCUSSION 


J. J. Ropcers anp M. L. Haviranp (General Motors Research 
Laboratories, General Motors Technical Center, Warren, Michi- 
gan): 


The author of this paper has presented a theory in which dif- 
fusion is the factor controlling boundary film formation with 
particular reference to “the formulation and evaluation of a theory 
of anti-squawk additive behavior in automatic transmission fluids.” 
Accordingly, the author presented relationships equating the rate 
of diffusion of molecules to the rubbing surfaces with various 
parameters. These relationships were then very ingeniously ap- 
plied to determine friction between lubricated surfaces. 

It is without question that if diffusion controls the movement 
of molecules to the rubbing surfaces, the equations set forth 
would be extremely valuable in clarifying boundary film forma- 
tion. Having realized the significance of the idea, we next wondered 
what present information we have, and what information we 
could obtain quickly to support or refute the ideas presented 
by the author. 

Our first thoughts were related to the development of the 
equations. In this respect, assumptions were both presented and 
implied by the author. One of these was that the additives re- 
acting with the surface lower the static friction. This was shown 
to be true by the author with stearic acid and di-n-octyl phosphite, 
and we have evidence which agrees with his findings. However, 
what would happen if an additive increased the static friction? 
Why do the equations relate only to the diffusion of static-friction- 
reducing additives? 

Another assumption was that, like the additive, oxygen dif- 
fuses to the surface and oxygen is necessary in the reaction of the 
additive with the surface. In fact, the oxygen concentration may 
well be the limiting factor in the additive reaction. The oxygen 
concentration of 0.0009 moles per liter is in a way suggested as a 
controlling limit in the film formation. However, the author has 
apparently ignored the oxide layer on the steel surface. Did the 
author attempt to remove this oxide layer from the steel surface? 
What about the paper surface? Is there also assumed to be a re- 
action at the paper surface? 

Although the assumptions can be questioned further, it may be 
better to proceed to a review of the equations as well as the 
results. The diffusion relationship, Eq. [9] in the paper, 


(Afa),, = (1 — exp (— CTd/pv)] 


is the most important one in this discussion. It takes the form 
of an exponential function which is asymptotic to a zero change 
in friction. To determine the change in friction (Afa) by this 
equation, one must measure static friction (Af,), the least re- 
peatable of the friction measurements. 

Now let us study the variables in the exponential term, (Ctd) / 
(uv). We would like to independently vary C, 7, and p and deter- 
mine their effect on the change in friction, but T and wp are not 
independent. However, it does appear that we can independently 
vary C, the concentration of “reaction rate controlling mole- 
cules” in moles per liter. 

The author has stated that in his work, oxygen is the reaction 
rate-controlling molecule for additive concentrations above 0.0009 
moles per liter, which is the estimated amount of molecular oxygen 
dissolved in the lubricant. Decreasing the concentration C by a 
factor of 10 (to 0.00009 moles per liter) and applying this value 
to the relationship for the 200 F data in Table 7A, the calculated 
Afa is increased over the entire sliding speed range, as shown in 
Fig. Al. 
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Fic. Al. Calculated and experimental concentration effects on 
friction change. 


In some of our studies, with 0.03 moles of stearic acid per liter 
of white oil, the oxygen content in the surrounding atmosphere 
was reduced so that it was estimated that no more than 0.00009 
moles of oxygen were present per liter of fluid. The experimental 
results are also shown in Fig. Al, and it is observed that the de- 
crease in oxygen concentration did not appreciably change Afa. 
These results are contrary to what was expected from Dr. Albert- 
son’s diffusion relationship. 

On the basis of our experimental results, the calculated changes 
in Afa shown in Fig. Al are incorrect. Of course, the author could 
say that the exponential term d varies to counterbalance any ef- 
fect produced by a change in the oxygen concentration, but, if so, 
how does he explain this? 

Now let us consider the term d, in the exponent. This term is 
defined as a constant equal to (N16)/(1000Ahn,,). It is not easy 
to determine each of the parameters in d, and as a consequence, 
the author apparently empirically determined d from Eq. [9] 
by substituting measured values of Afa and Af, and solved for 
values of d at various sliding speeds. 

Then, an average value for d was used to calculate the curves 
represented in many of the figures as dashed lines. Further, this 
same value for d was used in the determinations made at the 
other two temperatures. We wonder how the author decided 
which temperature to use for his calculation of d? 

Let us again consider the results the author obtained with 0.03 
moles per liter stearic acid in Fig. 10 and Table 5A. It appears 
that the term d was calculated from the 300F data to be 0.687. 
Using the author’s data we calculated d at each of the sliding 
speeds tabulated, and found that at 300 F, d increased with increas- 
ing sliding speed, as shown in Fig. A2. 

Although we calculated an average value of d which approxi- 
mates Dr. Albertson’s average value, we see that d increased from 
about 0.215 at a sliding speed of 5 fpm to about 0.92 at a sliding 
speed of 150 fpm. Rather than being constant, then, d increased 
fourfold. In fact, the relationship between d and sliding speed 
is very similar to that between Afa and sliding speed. 
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Fic. A2. Calculated factor, d, as a function of sliding speed 


After finding such a variation in d at 300 F, we became skeptical, 
so we calculated this “constant” at temperatures of 200 and 100F. 
At 200F, d varied from 0.65 to 1.50, giving an average value 
higher than 0.687. At 100 F, d varied from 0.569 to —5.03 and no 
averaging technique could possibly produce 0.687. Despite this 
variation in d, the author’s calculated data closely fit the measured 
values of friction. We submit that the reason this is so is that 
the equation is somewhat insensitive to changes in the exponential 
term. 

In Fig. A3, we have shown how sensitive Afa is to changes in 
the exponential term in Eq. [9]. Dr. Albertson’s experimental data 
as well as his calculated results are included. If we decrease the 
exponential term, (CTd)/p, by a factor of 2, we obtain the dotted 
curve; and if we increase the exponential term by a factor of 2, 
we obtain the dashed curve. The variation in the calculated data 
in Fig. A3 is no greater than the variation in some of Dr. Albert- 
son’s data. 

This last discussion brings up another point. According to Eq. 
[9], Afa values cannot exceed zero; that is, Afa cannot be posi- 
tive. However, with all of the compounds investigated by the 
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Fic. A4. Comparison of author’s and discussors’ data 


author, the data do exceed zero. Of course, it is argued that the 
repeatability accounts for this fact. 

However, in much of our work with our low velocity friction 
apparatus we have observed increases in friction with the addition 
of various compounds. In fact, using a 0.03 mole concentration of 
stearic acid per liter of white oil, we obtained the results presented 
in Fig. A4. We show the author’s experimental and calculated 
results for comparison. 

Some of the author’s and our experimental values of Afa exceed 
zero. As previously mentioned, the fact that the author’s results 
exceed zero can be explained on the basis of experimental error. 
Certainly his measured and calculated results are within repeat- 
ability limits. However, since our experimental results are signifi- 
cantly above zero, it is impossible to use Eq. [9] and arrive at 
results that are within experimental error. 

When using dimethylmethoxymethylphosphonate in the white 
oil, we obtained the results at 100 F shown in Fig. A5. Instead of 
reducing static friction, this particular additive seems to have a 
will of its own, producing an increase in static friction much 
higher than that obtained with the white oil. This is another situa- 
tion where Dr. Albertson’s diffusion relationship falls short; it is 
effective only for those additives which reduce static friction. 
How do additives such as dimethylmethoxymethylphosphonate 
reach the rubbing surfaces if not by diffusion? 

Another point that we want to discuss is presented in Fig. A6. 
We show both our data and the author’s data which were ob- 
tained at a lubricant temperature of 300F for a concentration of 
0.03 moles stearic acid per liter of white oil. In both cases, a 
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Fic. A5. Friction change obtained with addition of dimethyl- 
methoxymethylphosphonate to white oil. 
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Fic. A6. Friction change obtained with stearic acid 


high static friction was observed with a rapid decrease in friction 
on noticeable movement between the surfaces. Then, from about 
1 to 100 fpm the change in friction increased with increasing 
sliding speed. The minimum in both curves cannot be obtained 
using Eq. [9]. How does the author rationalize this fact? The 
fact that static friction is higher than kinetic friction is contrary 
to the author’s theory. According to the thedry, friction at zero 
sliding speed is lowest, because the film has infinite time to repair 
itself. 

Our conclusions regarding this paper can be stated rather 
simply. It has been proposed that diffusion controls friction, and 
the author has advanced relationships to this purpose. Although 
we believe that friction must be controlled, in part, by diffusion, 
we have shown that the equations are not adequate to explain 
the friction results obtained by the author as well as some of 
those obtained in sliding systems we have investigated. We think 
that the diffusion approach to frictional relationships is original 
and should be pursued further in order to further the understand- 
ing of frictional effects in lubricated sliding systems. 


W. P. Doyte (Texaco Research Center, Beacon, New York): 


Mr. Albertson’s correlation of the diffusion rates of anti-squawk 
agent and oxygen with the difference of the coefficient of friction 
between additive and nonadditive oil is very interesting. However, 
the function of oxygen in connection with anti-squawk action must 
be regarded as being of broader scope than indicated in the dif- 
fusion mechanism described here. : 

We have recently completed an experimental study of the ef- 
fects of oxidation on anti-squawk action. Some of the results of 
this study may be of interest here. A paper covering our work in 
this field is now under preparation and will be released when 
completed. 

Our work on anti-squawk agents has shown that oxidation 
of the base oil and the anti-squawk additive are extremely im- 
portant in determining the composition of the film absorbed on 
the clutch plates. The per cent monolayer coverage by the anti- 
squawk agent on the clutch plates was determined by means of 
a radioactive tracer method. It was found that significant amounts 
of polar oxidation products are formed in an ATF in only 2 hours 
at 275F. These products were found to absorb on the clutch 
plates, thus decreasing the coverage by the anti-squawk agent. It 
was also shown that the anti-squawk agent is itself oxidized to 
a nonadsorbable product. It should be pointed out that neutraliza- 
tion number data gave no indication of the build-up of polar 
oxidation products in the ATF during the early stage of oxida- 
tion. This means that neutralization number data may be mean- 
ingless in determining the oxidation stability of an ATF. The 
effects of oxidation could be retarded either by the addition of 
an antioxidant or by exhaustive refinement of the base oil. 

These results suggest two questions regarding the work in this 
paper. First, how long did it take to heat the oils to the desired 
temperature and how long were they held at these temperatures? 
Second, how oxidation stable were the blends used in these studies 
and how was this stability determined ? 


AUTHOR’s CLOSURE 


Rodgers’ and Haviland’s careful study and questions were ap- 
preciated. If an additive increased static friction, Af, would be 
positive and [6] and [7] should still be true. These general dif- 
fusion equations should apply equally to static friction increasing 
or decreasing additives. 

No attempt was made to remove the oxide from the steel surface 
because any oxide present was abraided and worn away during the 
wear-in period, see Fig. 3. No boundary lubricant was assumed 
to be on the paper surface in this simplified theory, although a 
reaction may be occurring. Since the paper material is elastic, 
surface fibers maintain contact with the steel surface and are not 
subject to the cyclic wipe-lubricate conditions on the steel surface. 

Similar oxygen insensitive test results were reported in the 
section “Frictional Behavior at Higher Additive Concentrations.” 
It is felt that moisture may be the controlling agent instead of 
oxygen. However, the circumstantial nature of this hypothesis was 
emphasized since no definitive test data are available on the 
lubrication mechanism. This area requires more work. 

The “d” factor was calculated from 300F data at low zinc 
stearate concentrations because diffusion is greater at higher tem- 
peratures. 

That dynamic Afa data in Fig. A4, as well as our own data, 
were greater than zero may be causec by decreased lubricant 
viscosity due to the additive, or to other systemic changes between 
the oil lubricated and oil + additive lubricated runs. The theory 
assumes that no changes of this type are occurring. 
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The diffusion equations are general derivations which should 
apply to the phosphonate data. Since Af, is positive in Fig. AS, 
Equations [6] and [9] would reproduce the phosphonate 
curve. 

Data in Fig. A6 are due to systemic changes not accounted for 
in the simple theory. If the resin covering the paper fibers plastic- 
ally conformed to the steel plate after static engagement, friction 
would increase as a function of engagement time. Since static 
friction was generally time dependent in cases of abnormally high 
friction, it is felt that the additive was plasticizing the resin on 
the lining surface. 


In answer to Mr. Doyle’s questions, it took about 10-15 minutes 
to reach 200F and about 15 minutes to go from 200 to 300F. 
Test duration was usually 4 to 1 hour. If static friction was chang- 
ing or results were other than expected, the lubricant was replaced 
and the run repeated. When equipment was at 300 F, a run could 
be made about 5 minutes after the addition of fresh lubricant. 

A white, refined mineral oil was used as the base fluid. Oxida- 
tion stability was not checked chemically. Stability was checked 
by observing static friction. If static friction varied more than 
+ 0.02, or in an unexpected manner, the run was repeated with 
fresh oil or fresh lubricant. 
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Oxidation Stability of the Fluoroalcohol Esters 
By JONATHAN L. SNEAD! and HENRY GISSER (aste)' - 


Oxidation studies of fluoroalcohol esters indicate that elimination of secondary and tertiary C—H . 
bonds in the acid portion of the molecule considerably improves oxidation stability. 

The trimethylacetates of tetrafluorobutane-1,4-diol, hexafluoropentane-1,5-diol, and octafluoro- E 
hexane-1,6-diol, when subjected to air oxidation at 250 C for 168 hours in the presence of steel, 
copper, and aluminum catalyst strips, had neutralization numbers and viscosity increases of 4.7 H 
and 2.8%, 3.7 and 3.6%, and 5.4 and 1.8%, respectively. Bis(1H,1H,7H-dodecafluoroheptyl) 3- 
methylglutarate gave values of 5:1 and —6.6% after 24 hours, and bis(1H,1H,7H-dodecafluoro- H 
heptyl) dl-camphorate gave values of 8.9 and —15% after 72 hours. Hexafluoropentanediol-1,5 
bis(n-hexanoate) had a neutralization number of 66 and viscosity increase of 38% after 20 hours B 


oxidation at 175 C. The fluoroglycol esters of 2,2- and 3,3-dimethylbutyric acids and bis(1H,1H,7H- 
dodecafluoroheptyl) 3,3-dimethylglutarate were also considerably less stable than the fluoroglycol 


trimethylacetates. 


The physical properties and oxidation and thermal stabilities of the esters having no tertiary or 
secondary C—H bonds in the acid moiety suggest their use as base fluids for lubricants, heat 
transfer media, coolants, and other applications where stability at high temperatures is essential. 


Introduction 


Tue high temperature use of organic fluids is limited by 
oxidation and thermal stability. In the case of esters, 
considerable improvement in these properties has been 
achieved by the use of fluoroalcohols in the esters (/). 
The use of fully fluorinated esters is not feasible because 
esters from fluorinated acids are inherently unstable (2). 
Since both oxidation and thermal stability are dependent 
on molecular configuration, suitable limitations on the 
configuration of the fluoroalcohol esters should result in 
further enhancement of oxidation and thermal stability. 
It has been indicated that in the oxidation of hydro- 
carbons (3), oxidation occurs preferentially at a tertiary 
C—H bond, and evidence is available that the same is 
true of the saturated aliphatic diesters (4). In general, 
the order of increasing resistance to oxidative attack at 
C—H bonds is tertiary, secondary, and primary. This is 
to be expected from the free radical nature of the oxida- 
tion process and the fact that the relative probability of 
attack of alkyl radicals at C—H bonds decreases in the 
order tertiary, secondary, and primary (5). The elimina- 
tion of secondary and tertiary C—H bonds in the acid 
moiety of the fluoroester should result in an increase in 
oxidation and thermal stability. 

Diesters with improved thermal and ‘hydrolytic stabil- 


Contributed by the ASLE Technical Committee on Properties 
of Lubricants and presented at the Annual Meeting of the 
American Society of Lubrication Engineers held in New York, 
April, 1963. 

1 U.S. Army Munitions Command, Frankford Arsenal, Pitman- 
Dunn Laboratories, Philadelphia 37, Pennsylvania. 
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ity were obtained by the introduction of alkyl groups on 
the beta carbon of the alcohol portion, and the alpha 
carbon of the acid portion of the molecule (6). The effect 
on thermal stability was explained in terms of blocking 
the cyclic transition state (described (7) as a step in the 
thermal decomposition mechanism) formed by removal 
of a beta hydrogen by the carbonyl oxygen. The hydro- 
lytic stability improvement was explained to be a result 
of steric hindrance in the acid function. Later, it was 
found that thermal stability also was improved by alkyl 
groups on the alpha (or beta) carbon atoms of the acid 
radical in diesters (8). Further blocking of the next de- 
composition path (probably a free radical chain mecha- 
nism), it was suggested, may be achieved from esters of 
neo-alcohols or glycols and “hindered” acids. Other 
workers have found that esters derived from highly 
hindered acids and phenols, and esters containing the 
neopentyl glycol structure, possess improved stability 
(9). The importance of the neopentyl structure to 
thermal stability was emphasized in earlier work on such 
varied structures as the pentaerythritol and trimethylol- 
propane esters (10, 11) and esters of tert-alkylcarbinols 
and 2,2-difluoroalkanols (12). While these achievements 
in improved thermal stability were considerable, the 
gains made in increased oxidation stability were less sig- 
nificant and left much to be desired. 

Neopentane has no secondary or tertiary C—H bonds, 
and removal of such bonds from organic molecules is 
accomplished by the introduction of the neopentyl (or 
neo) structure. It should be noted however, that the 
effect of the neo-structure is associated with reduction of 
free radical reactions when these pertain (oxidation and 
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TABLE 1 
Dynamic Oxidation Studies of the Esters 
Viscosity 
Neutralization No. Metal Wt. Change 
period at 37.8C (eq mg KOH/gm) (mg/sq cm) Sample 
Compound (hr) (%) Volatile Non-Vol Cu Steel Al appearance 
at 175 C 
Hexafluoropentanediol-1,5 < 20 37.8 39.3 26.7 + 0.07 + 0.01 + 0.04 Dark; 
bis (n-hexanoate) sludge 
Hexafluoropentanediol-1,5 > 168 0 0.06 0 + 0.04 — 0.03 —0.01 Clear 
bis(trimethylacetate) 
Hexafluoropentanediol-1,5 > 168 12.6 2.6 24 —0.5 — 0.02 0 Amber 
bis(2,2-dimethylbutyrate) 
Hexafluoropentanediol-1,5 > 168 9.6 3.5 5.4 — 0.2 0 0 Turbid 
bis(3,3-dimethylbutyrate) 
Bis(1H,1H,5H-octafluoropentyl) > 168 6.9 3.6 2.5 — 0.03 + 0.04 + 0.08 Amber 
3-methylglutarate 
Bis(1H,1H,7H-dodecafluorohepty]) > 168 3.1 1.6 1.2 — 0.05 + 0.05 + 0.09 Amber 
3-methylglutarate 
Bis(1H,1H,7H-dodecafluorohepty]) > 168 2.4 0.5 0.1 —0.1 + 0.01 +- 0.03 Clear@ 
dl-camphorate 
at 200C 
Octafluorohexanediol-1,6 > 168 0.1 0.2 0.6 + 0.05 + 0.04 0 Clear 
bis(trimethylacetate) 
Octafluorohexanediol-1,6 168 14.2 9.5 1.8 —0.7 + 0.07 + 0.08 Turbid 
bis(3,3-dimethylbutyrate) 
at 225 C 
Bis(pentafluoropropy]) 168 11.6 17.8 16.4 —04 +01 +0.1 
tetramethylsuccinate 
Bis(heptafluorobuty]) > 168 1.8 4.1 11.3 + 0.01 + 0.07 
tetramethylsuccinate 
at 250 C 
Tetrafluorobutanediol-1,4 > 168? 2.8 1.9 2.8 — 0.6 0 + 0.04 
bis(trimethylacetate) 
Tetrafluorobutanediol-1,4 < 25 3.9 8.6 + 0.2 +0.4 
bis(3,3-dimethylbutyrate) 
Hexafluoropentanediol-1,5 > 168? 3.6 2.6 1.1 —0.2 +-0.01 + 0.01 
bis(trimethylacetate) 
Hexafluoropentanediol-1,5 < 32 10.7 7.4 —1.3 +0.1 + 0.2 
bis(3,3-dimethylbutyrate) 
Octafluorohexanediol-1,6 > 168? 1.8 4.5 0.9 —- 0.06 + 0.06 + 0.09 
bis (trimethylacetate) 
Octafluorohexanediol-1,6 < 24 42 8.6 4.5 —13 +0.1 + 0.2 
bis(3,3-dimethylbutyrate) 
Bis(1H,1H,7H-dodecafiuoro- < 24 — 6.6 3.7 14 — 0.3 — 0.03 — 0.02 
heptyl) 3-methylglutarate 
Bis(1H,1H,7H-dodecafluoro- 168 —11.1 16.3 5.2 —0.5 + 0.2 +0.1 
heptyl) 3,3-dimethylglutarate 
Bis(1H,1H,7H-dodecafluoro- 72 — 14.9 6.0 2.9 — 0.6 0 — 0.04 


heptyl) dl-camphorate 


@ Light sediment. 


b The induction period was 500-600 hours. 
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probably thermal stability where a cyclic path is not 
available), and in making unavailable a cyclic path for 
thermal decomposition. 

In the work described in this paper, several fluoro- 
glycol esters of neopentanoic and neohexanoic acids, and 
dibasic acid esters of fluoroalcohols with tetramethylsuc- 
cinic acid were prepared and their thermal stability 
studied. The neopentanoic and tetramethylsuccinic acid 
esters contain only primary C—H bonds in the acyl 
group, while the esters of neohexanoic acid contain, in 
addition, several secondary C—H bonds. The relative 
oxidation and thermal stabilities of these compounds 
were determined and compared with those of esters of 
fluoroalcohols containing both secondary and tertiary 
C—H bonds in the acyl group. 


Oxidation and thermal stability 


Preparation and properties of the esters are described 
in the appendix. 

Oxidation studies were made beginning at 175 C and in 
25 C increments up to 250 C. Ten-gram samples of the 
esters, in duplicate, were subjected to oxidation (13) in 
the presence of coupled copper, steel, and aluminum 
catalyst strips. Clean, dry air was bubbled at the rate of 
1.2 liters per hour through the sample in an all glass 
apparatus and the effluent gas passed through a flask 
containing a measured volume of standardized KOH. The 
usual criteria of changes in viscosity, acidity, color, and 
catalyst weight were employed to determine stability. The 
data are summarized in Table 1. 
Hexafluoropentanediol-1,5 bis(trimethylacetate), at 
175 C, showed definite signs of possessing greater stability 
than the other esters at this temperature. Hexafluoro- 
pentanediol-1,5 bis(m-hexanoate), which contains a large 
number of secondary C—H bonds, was severely oxidized 
in a few hours. Beginning at 200 C, differences in stabil- 
ity between the octafluorohexanediol esters of trimethyl- 


acetic acid, containing no secondary or tertiary C—H 
bonds in the acyl group, and 3,3-dimethylbutyric acid, 
containing several secondary bonds became more pro- 
nounced, the former showing greater stability. 


At 250 C, the difference in stability between esters 
containing only primary C—H bonds in the acid portion 
of the molecule, and those containing secondary and 
tertiary C—H bonds is quite significant. Thus, the sta- 
bility of the trimethylacetates is far superior to that of 
the 3,3-dimethylbutyrates; the 3,3-dimethylglutarate is 
superior to the 3-methylglutarate but inferior to the tri- 
methylacetates. Bis(1H,1H,7H-dodecafluoroheptyl)  dl- 
camphorate, which has a tertiary C—H bond, also shows 
lesser stability even though some increased stability is 
derived from its cyclic structure. The trimethylacetates, 
after 168 hours, were practically unchanged, being clear, 
colorless, and without sludge or sediment. The related 
3,3-dimethylbutyrates on the other hand, showed, after 
only 24 hours, high increases in viscosity and acidity as 
well as the formation of insoluble substances. 


Oxidation of bis(heptafluorobutyl) tetramethylsucci- 
nate at 225 C indicated a high degree of stability com- 
parable to that shown by the trimethylacetates at 200 C. 
Oxidation of 
could not be determined at 225 C, however, because of 
its relatively high volatility at that temperature. 

Limited thermal stability studies were made on repre- 
sentative esters. About 5 gm of each ester, in duplicate, 
was sealed under vacuum in glass ampoules previously 
flushed with nitrogen, and heated at 260 C for 72 or 168 


_ hours. Changes in viscosity and acidity were used as 


criteria for stability. Although none of the esters decom- 
posed severely (not more than 0.5% calculated on the 
assumption that the increase in acidity was due solely to 
regeneration of the parent acid), it is evident from 
Table 2 that the trimethylacetates are superior to the 
other esters. 


TABLE 2 
Thermal Stability of the Fluoroalcohol Esters 


heptyl) 3,3-dimethylglutarate 


Viscosity 
increase Decomposi- 
Hours at 37.8 C Neut. No. tion 
Compound at 260 C (%) increase (%) 
Hexafluoropentanediol-1,5 72 0 0 Nil 
bis(trimethylacetate) 
Hexafluoropentanediol-1,5 72 + 0.6 0.7 0.2 
bis (3,3-dimethylbutyrate) 
Tetrafluorobutanediol-1,4 168 0 0 Nil 
bis(trimethylacetate) 
Tetrafluorobutanediol-1,4 168 +2 0.8 0.2 
bis(3,3-dimethylbutyrate) 
Bis(1H,1H,7H-dodecafluoro- 168 — 2.3 0.6 0.4 
heptyl) 3-methylglutarate 
Bis(1H,1H,7H-dodecafluoro- 168 —3.1 0.7 0.5 
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The observed oxidative and thermal stability of organic 
fluids is determined by a series of chemical reactions, 
and the over-all rate reflects the rates of those reaction 
steps which are slowest. The process from initial to final 
products may take place by various paths, that path 
being preferred which requires lowest energy for the 
intermediate configurations (activated complexes). Since 
the energy for formation of the intermediate complexes 
is of thermal origin, the less readily accessible an acti- 
vated complex corresponding to a mode of decomposition, 
i.e., the higher energy required for its formation, the 
higher the temperature required for decomposition. 


In oxidation processes (of organic fluids), it is gen- 
erally accepted that the mechanism involves a radical 
chain reaction. Since the first isolatable product of oxida- 
tion (of a hydrocarbon configuration) is a hydroperoxide, 
the decomposition mode of the hydroperoxide provides 


the path for further oxidation. From earlier work on the 
oxidation of alkylbenzenes (14-16) it may be deduced 
that a probable mode of decomposition of primary and 
secondary hydroperoxides is dehydration with the forma- 
tion of aldehydes and ketones. Tertiary hydroperoxides 
on the other hand, may decompose by scission of the 
ROOH into RO- and -OH, the former then further split- 
ting by rupture of a C—C bond yielding a carbonyl com- 
pound and a free radical of lower molecular weight (17). 
In contrast to primary and secondary hydroperoxide de- 
composition, a free radical is now available for initiation 
of new oxidation chains. The kinetics of the over-all 
process makes it likely that the bulk of the observed 
oxidation is due to initiation of new chains by hydroper- 
oxide decomposition. Thus, in spite of the fact that ter- 
tiary hydroperoxides are more stable than primary or 
secondary hydroperoxides, compounds having tertiary 
C—H bonds would be expected to oxidize at a much 
more rapid rate than compounds which do not. The dif- 


TABLE 3 
Viscometric and Low Temperature Properties 


Viscosity (cs) a 
Compound 98.9C $7.8C —17.8C — 28.9C —40C — 53.9C Ce 

bis(trimethylacetate) 

Tetrafluorobutanediol-1,4 2.00 8.41 8 
bis(3,3-dimethylbutyrate) 

Hexafluoropentanediol-1,5 1.70 6.63 159 509 2152 22,000¢ — 43 
bis(trimethylacetate) 

Hexafluoropentanediol-1,5 2.23 10.03 266 837 3441 — — 42 
bis(3,3-dimethylbutyrate) 

Octafluorohexanediol-1,6 1.69 7.33 182 591 2600 28,000 <— 54¢ 
bis(trimethylacetate) 

Octafluorohexanediol-1,6 2.30 10.80 314 987 4074 39,631 < — 54¢ 
bis(3,3-dimethylbutyrate) 

Hexafluoropentanediol-1,5 1.67 5.53 54.8 118 3004 — —31 
bis(n-hexanoate) 

Hexafluoropentanediol-1,5 1.99 8.30 185 560 2071 15,680 <— 54¢ 
bis(2,2-dimethylbutyrate) 

Bis(pentafluoropropyl) 1.14 3.88 57.6 168 586 5118 < — 54¢ 
tetramethylsuccinate 

Bis(heptafluorobuty]) 1.40 6.10 153 470 1855 17,750 < — 54¢ 
tetramethylsuccinate 

Bis(1H,1H,7H-dodecafluoro- 3.700 26.40 1280° 42700 19,0007.» — 
heptyl) 3-methylglutarate 

Bis(1H,1H,7H-dodecafluoro- 3.83 28.4 1561 5505 24,338 — — 43¢ 
heptyl) 3,3-dimethylglutarate 

Bis(1H,1H,5H-octafluoro- 2.790 15.90 5220 15102 5360? 40,8002 — 59».¢ 
pentyl) 3-methylglutarate 

Bis(1H,1H,7H-dodecafluoro- 8.224 115.54 — 32¢4 


heptyl) dl-camphorate 


@ Determination made on supercooled liquid. 
> Values obtained from Ref. 1. 
¢ Pour point. 


@ Values obtained from manufacturer’s technical information bulletin. 
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ferences in oxidation rate between primary and secondary 
C—H bonds is probably due mainly to ease of formation 
of a free radical at the site (this has been discussed 
earlier). By elimination of tertiary and secondary C—H 
bonds from the structure, one would therefore expect that 
the oxidation stability would be increased and the data 
reported above indicate that such is the case. 

The fluoroalcohol esters possess secondary C—H bonds 
on the carbon atom alpha to the ester group. It has been 
shown in earlier work (4, 18) that oxidation rate at a 
site increases with increasing electron density, and that 
in case of esters, resonance at the ester group reduces the 
electron density at neighboring carbon atoms (the effect 
falling off as distance from the ester group increases). 
Oxidation at a C—H site alpha to the ester group would 
therefore probably contribute little to overall rate. 

Apart from oxidation and thermal stability, the use- 
fulness of the compounds described above depends on 
their physical properties, primarily fluid range and vis- 
cosity. Viscosity and low temperature data are given in 
Table 3, and these appear to be within ranges applicable 
to hydraulic fluid and lubrication applications. 

The physical properties and oxidation and thermal 
stabilities of the fluoroalcohol esters having no tertiary 
or secondary C—H bonds in the acid moiety suggest 
their use as base fluids for lubricants, heat transfer 
media, coolants, and other applications where stability at 
high temperatures is essential. 
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Appendix 
PREPARATION AND PROPERTIES OF THE ESTERS 


Starting Materials 


Hexafluoropentanediol-1,5 was recrystallized from 
chloroform containing a little alcohol. Pentafluoropro- 
panol, heptafluorobutanol, and 1H,1H,7H-dodecafluoro- 
heptanol were purified by distillation through a 24-inch 
Podbielniak Heli-Grid column, at a 10:1 reflux ratio. 
Trimethylacetic acid, n-caproic acid, tert-amyl chloride, 
and 3,3-dimethylbutyric acid were fractionated through 
a 36-inch heated column filled with glass helices at a 4:1 
reflux ratio. 3-Methyl, and 3,3-dimethylglutaric acids 
were recrystallized, the former from benzene-petroleum 
ether (3:2), and the latter from benzene. Ethyl per- 
fluorosuccinate, ethyl perfluoroadipate, and 2,2’-azodi-iso- 
butyronitrile were used without further purification. 

2,2-Dimethylbutyric acid was prepared by the reaction 
of tert-amylmagnesium chloride with carbon dioxide 
(19). Perfluoroglycols were prepared by reduction of 
ethyl perfluorosuccinate and ethyl perfluoroadipate with 
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lithium aluminum hydride (20). Tetramethylsuccinic 
acid was prepared by a 3-step synthesis beginning with 
the thermal decomposition of 2,2’-azodi-iso-butyronitrile 
in toluene (21). The tetramethylsuccinonitrile formed 
was hydrolyzed in alkaline solution to the imide, which 
was converted to the desired compound by acid hydrol- 


TABLE 4 
Intermediates 
Boiling Melting 
point point 
Compound (°C/mm Hg) (°C) np? 
Tetrafluorobutanediol-1,4 85 
Hexafluoropentanediol-1,5 77-78 
Octafluorohexanediol-1,6 68 
Ethyl perfluorosuccinate 89/15 
Ethy] perfluoroadipate 70-71/2.5 1.3541 
Pentafluoropropanol 81-83 /760 
Heptafluorobutanol 95/749 
1H,1H,7H-Dodecafluoro- 82/20 1.3156 
heptanol 
n-Caproic acid 105.7/7 1.4163 
Tetramethylsuccinic acid 200 
Trimethylacetic acid 161/760 34 
2,2-Dimethylbutyric acid 187/760 
3,3-Dimethylbutyric acid 109/49 1.4096 
3-Methylglutaric acid 84-84.5 
3,3-Dimethylglutaric acid 101.5-102.0 
t-Amyl chloride 86/760 
2,2’-azodi-iso-butyronitrile 103 (dec.) 


ysis. The melting and boiling points and refractive index 
(liquids) of the intermediates after purification are shown 
in Table 4. 


Esters of Fluoroglycols 


The fluoroglycol esters were prepared by direct ester- 
ification in the presence of catalytic quantities of con- 
centrated sulfuric acid (22). A detailed description of 
the preparation of hexafluoropentanediol-1,5 bis(tri- 
methylacetate) follows and illustrates the procedure. 
Trimethylacetic acid, 61.3 gm (0.60 mole), hexafluoro- 
pentanediol-1,5, 53 gm (0.25 mole), 2 gm of H.SO,, and 
250 ml of toluene were placed in a liter, two-necked, 
round-bottom flask fitted with a mechanical stirrer, 
condenser, and a Dean and Stark water trap. A 12-inch 
Vigreux column was placed between the boiling flask and 
the water trap. The reaction mixture was stirred and 
refluxed until the amount of water collected was close to 
the theoretical (9-10 hours). The crude ester was diluted 
with 500 ml of ether, washed successively with VY KOH, 
water, 0.01 N HCl, and distilled water until neutral. 
Distillation over 0.5 gm of anhydrous potassium car- 
bonate, fractionation twice through a 12-inch Vigreux 
column, and finally percolation through a column con- 
taining alumina and Fuller’s earth yielded a water-white 
liquid with neutralization number less than 0.01. A yield 
of 85.5 gm (90%) of ester was obtained after purifica- 
tion. 
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Esters of Fluoroalcohols 


Esters of 3-methyl- and 3,3-dimethylglutaric acids 
were prepared directly in a manner similar to above. 
Tetramethylsuccinate esters could not be prepared by 
direct esterification since conversion of the acid to an- 
hydride proceeded at a faster rate on heating. The anhy- 
dride distills with water and as a result, the acid is 
removed from the reaction with the water-solvent aze- 
otrope. These esters were prepared by reaction of tetra- 
methylsuccinyl chloride with sodium fluoroalcoholates. 

The acid chloride was prepared by the dropwise addi- 
tion of excess thionyl chloride to solid tetramethylsuccinic 
acid (molar ratio 6:1) with agitation, at a rate which 
allowed the reaction flask to be maintained at 0-10 C in 
an ice-salt bath. After the addition, stirring was continued 
at room temperature for 3 days until the entire mass was 
liquid. The excess thionyl chloride was distilled at reduced 
pressure leaving a viscous, dark brown liquid, which was 
used without further purification. 

The preparation of the fluoroalcohol esters is illustrated 


by the following procedure for bis(heptafluorobutyl) 
tetramethylsuccinate. Eighty grams (0.4 mole) of hepta- 
fluorobutanol as its sodium salt in 250 ml of ether was 
added with stirring to 21.1 gm (0.1 mole) of tetramethyl- 
succinyl chloride and the mixture refluxed 2—3 hours. The 
crude ester was washed successively with water, 10% 
potassium carbonate, and 0.01 N HCl, then further 
purified in a manner similar to the esters of fluoroglycols 
above. Repeated fractional distillations were necessary 
to remove small amounts of tetramethylsuccinic anhy- 
dride. A yield of 23.7 gm (44%) was obtained after 
purification. 


The boiling points, elemental analyses, and saponifica- 
tion equivalents are given in Table 5. It can be seen that 
the boiling points within a homologous series increases 
by approximately eight degrees for each additional per- 
fluoromethylene group. A similar increment has been 
found for the perfluorobutyl esters of adipic through 
sebacic acids (23). Complete saponification required 24— 
48 hours, indicating that esters of fluoroalcohols, as a 


TABLE 5 
Analytical Data on the Fluoroalcohol Esters 


Boiling Saponification 
point Carbon (%) Hydrogen (%) Fluorine (%) equivalent 
Compound (°C/mm Hg) Calc. Found Calc. Found Calc. Found Calc. Found 

Tetrafluorobutanediol-1,4 78/0.35 50.39 50.99 6.72 6.85 23.01 22.89 165.1 165.3 
bis(trimethylacetate) 

Tetrafluorobutanediol-1,4 100/0.35 53.63 53.77 7.31 7.44 21.21 20.97 179.2 179.2 
bis(3,3-dimethylbutyrate) 

Hexafluoropentanediol-1,5 86/0.35 47.38 47.67 5.83 5.79 29.98 29.68 190.1 190.5 
bis(trimethy!acetate) 

Hexafluoropentanediol-1,5 108/0.35 50.00 50.35 6.42 6.66 27.91 27.33 204.4 204.4 
bis(3,3-dimethylbutyrate) 

Octafluorohexanediol-1,6 94/0.35 44.67 44.75 5.15 5.43 35.32 35.10 215.1 215.0 
bis (trimethylacetate) 

Octafluorohexanediol-1,6 116/0.35 47.16 47.30 5.72 5.83 33.16 32.05 229.2 229.2 
bis(3,3-dimethylbutyrate) 

Hexatluoropentanediol-1,5 137/0.8 50.00 50.43 6.42 6.44 27.91 27.05 204.4 204.2 
bis(-hexanoate) 

Hexafluoropentanediol-1,5 130/0.62 50.00 50.28 6.42 6.69 27.91 27.55 204.4 204.2 
bis(2,2-dimethylbutyrate) 

Bis(pentafluoropropy]) 82/2.15 38.35 38.46 3.68 3.74 43.35 43.18 219.1 219.2 
tetramethylsuccinate 

Bis(heptafluorobuty]) 74/0.55 35.70 35.67 3.00 32a 49.42 48.84 269.1 269.2 
tetramethylsuccinate 

heptyl) 3-methylglutarate 

Bis(1H,1H,7H-dodecafluoro- 143/0.33 31.99 32.07 2.05 2.14 57.84 57.34 394.1 394.6 
heptyl) 3,3-dimethylglutarate 

pentyl) 3-methylglutarate 


heptyl) dl-camphorate 


@ Values obtained from Ref. 1. 


+’ From manufacturer’s technical information bulletin. 
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TABLE 6 
Physical Properties of the Fluoroalcohol Esters 


Molar refraction 


(MR) 
Compound Calc.4 Obsvd. 
Tetrafluorobutanediol-1,4 1.1192 1.3930 70.6 70.4 
bis(trimethylacetate) 
b Tetrafluorobutanediol-1,4 1.0985 1.4047 79.9 79.9 
bis(3,3-dimethylbutyrate) 
Hexafluoropentanediol-1,5 1.1852 1.3854 75.3 75.3 
bis(trimethylacetate) 
Hexafluoropentanediol-1,5 1.1609 1.3973 84.6 84.8 
bis(3,3-dimethylbutyrate) 
Octafluorohexanediol-1,6 1.2512 1.3778 80.0 79.2 
bis(trimethylacetate) 
Octafluorohexanediol-1,6 1.2233 1.3899 89.3 88.8 
bis(3,3-dimethylbutyrate) 
Hexafluoropentanediol-1,5 1.1676 1.3996 84.6 84.7 
bis(n-hexanoate) 
it Hexafluoropentanediol-1,5 1.1647 1.3961 84.6 84.3 
bis(2,2-dimethylbutyrate) 
Bis(pentafluoropropyl) 1.3607 1.3635 70.8 71.7 
tetramethylsuccinate 
Bis(heptafluorobutyl) 1.4512 1.3543 80.2 80.7 
tetramethylsuccinate 
Bis(1H,1H,7H-dodecafluoro- 1.6484? 1.3505? 99.0 101.2 
heptyl) 3-methylglutarate 
Bis(1H,1H,7H-dodecafluoro- 1.6140 1.3528 103.6 105.9 
heptyl) 3,3-dimethylglutarate 
Bis(1H,1H,5H-octafluoro- 1.55702 1.35950 80.2 81.3 
pentyl) 3-methylglutarate 
Bis(1H,1H,7H-dodecafluoro- 1.597 1.37090 


heptyl) dl-camphorate 


@ Bond refractivities of Denbigh (25) were used in the calculations. 
» Literature value for density and refractive index at 20 C. 


class, show much greater hydrolytic stability than esters 

in general. Esters derived from fluoroacids in particular, 

exhibit very poor stability towards alkaline hydrolysis 

(2, 24). 3. 
Table 6 shows that the density increases with increas- 

ing fluorine content within a homologous series while the 

refractive index decreases, as expected. The observed and 

calculated molar refractivities (Table 6) are in excellent 

agreement with fluorine content up to about 40%. With 

more than 40% fluorine content, however, the values be- 

come more divergent. This phenomenon has been attrib- 

uted to variation in the atomic refraction of the C—F 

bond (1). 
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A First Order Solution for the Stress Concentration 
Present at the End of Roller Contact 


By C. A. MOYER! and H. R. NEIFERT? 


An approximate solution of the “critical” maximum compressive stress present in rolling contact 
fatigue is given. Geometries normally used in rolling contact rig tests having contact areas varying 
from rectangular to full elliptical contact are considered. For cylinders in contact it is shown that 
the first order solution for end of contact stress is considerably above the nominal line contact 
Hertz stress usually calculated. Confirmation of the approximate stresses calculated is provided by 
deformation studies just into the plastic regime, and rolling contact fatigue tests on both cylindrical 


and 100-inch radius toroidal specimen rollers. 


Nomenclature 


a = Semi-major axis for a full elliptical contact 
area (see Fig. 3), inches 

bo = Semi-minor axis for an elliptical or a trun- 
cated elliptical contact area (see Fig. 3), 
inches 

bo = One-half width of contact for a rectangular 
contact area (see Fig. 2), inches 

b, = One-half width of contact at the truncated 
edge for a truncated ellipse (see Fig. 3), 


inches 
a’, b’ = Semi-axes of the roller end ellipse (see Fig. 
4b), inches 
L = Contact length between two rollers (see Fig. 
1), inches 


1 = Semi-length of actual contact between two 

rollers (1 = L/2), inches 

P» = Applied load, pounds 

P;,= Applied load giving contact ellipse with 
major axis = /, pounds 

P = Equivalent load for a truncated ellipse, 
pounds 

P’ = Load necessary to make a’ = bo (cylinder 
against cylinder) or a’ = 0, (at least one 
roller a toroid), pounds 

Qo = The maximum compressive stress in psi cal- 
culated using Eqs. [1], [3], or [7] (see 
Figs. 1 or 2) 


Contributed by the ASLE Technical Committee on Bearings 
and Bearing Lubrication and presented at the Annual Meeting 
of the American Society of Lubrication Engineers held in New 
York, April, 1963. 

1 Senior Research Engineer, Physical Laboratories, The Timken 
Roller Bearing Company, Canton, Ohio. 

2 Chief Engineer, Physical Laboratories, The Timken Roller 
Bearing Company, Canton, Ohio. 
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gi: = The maximum compressive stress in psi at 
the truncated end of contact, calculated using 
Eq. [10] (see Fig. 3) 

Qe = The maximum compressive stress concen- 
tration in psi at the end of contact, calculated 
using Eqs. [12] or [13], in psi (see points 
A and C in Fig. 4) 

q = The maximum compressive stress in psi 
related to plastic deformation (taken from 
the line in Fig. 11) 

R,, Re = Radii of cylinders (or toroids) in the direc- 
tion of rolling, inches (see Fig. 1) 
R,’, R»’ = Contour radii of toroids perpendicular to the 
rolling direction, inches (see Fig. 4a) 
~1/R = Summation of the curvatures of the con- 
tacting roller bodies (1/R; + 1/R2 + 
1/Ri’ + 1/R2’) 

r;’ = End radius of shorter roller perpendicular to 
the rolling direction, inches (see Fig. 4a) 

A= Permanent vertical displacement of roller 
surface from its original profile, microinches 

, v = Transcendental functions related to the axes 
of a contact ellipse 
u’, v = For example v/u = bo/ao. Values of the 
functions can be found in the Appendix, Fig. 
14 


Introduction 


In the theory of elastic contact stresses, special attention 
has been given to the problems of line contact between 
two finite length elastic cylinders. The solutions for this 


problem usually start with plane strain assumptions, and 


the work of Hertz (1), Belajef (2), and Thomas and 
Hoersch (3) on line contact are classic examples of this 
approach. All three sources provide equations usable in 
calculating the static stresses present in the contact region. 
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Radzimovsky (4), Poritsky (5), and Smith and Liu (6) 
extended the equations to map the complete internal stress 
field for the dynamic condition of rolling cylinders. They 
consider a “hydrodynamic” analogue with the material 
“flowing” through a calculated stress lattice below the 
contact surface and bounded by the contact region. The 
latter two investigators have included the effect of both 
rolling and sliding on the calculated stresses. 

The original contact stress problem can also be am- 
plified to consider stress concentrations occurring at the 
axial end of contact when one finite length cylinder is 
longer than the other (Fig. 1). Solutions for this problem 
have been proposed (7, 8). For the further problem of 
two contoured finite rollers, or one contoured and one 


Fic. 1. Two cylinders in contact with the contact length de- 
termined by cylinder 1. 


cylindrical roller loaded against each other apparently 
nothing has been written. In this presentation first order 
solutions to both latter problems are offered, with par- 
ticular consideration to the stress concentrations present 
at the ends of the finite roller contact. In addition some 
rolling contact loading investigations just into the plastic 
regime are presented as corroborative evidence supporting 
the approximate solutions. The roller geometries discussed 
are of importance since a great many fatigue tests are 
presently being conducted on rolling bodies with these 
general contours. 


Cylinders in contact 


Following the methods outlined by Thomas and 
Hoersch (3) or Belajef (2), the contact stress problem 
of two finite cylinders in contact will be examined. Under 
no load, line contact is produced. As load is applied to 
press the two cylinders together, the line broadens into 
an assumed rectangular area of contact. 

For this rectangular area the stress distribution within 
either body over the contact region may be represented by 
one-half of an elliptical cylinder, as shown in Fig. 2. The 


Fic. 2. Stress profile over a rectangular contact area produced 
by two cylinders in contact. 


maximum compressive stress go and the one-half width of 
contact 5» can be calculated from the equations: 


Jo = 2290 (po/L)* (1/Ri+ 1/R2)* 


Po 
L(1/R; + 1/R2) 


by = 0.0002779 | 


Toroid rollers in contact 


Whenever a roller is no longer cylindrical but has a 
contour radius along its length L, it is termed a toroid. 
Two toroidal rollers, or a toroidal roller and a cylindrical 
roller, under no load produce point contact. As load is 
applied, an elliptical contact area is formed for which the 
major axis 2a of the ellipse may be at right angles to the 
direction of rolling. For the problems to be considered it 
will be assumed the geometries are such that the major 
axis will always be oriented in this direction. 

The equation for calculating the contact stress go for 
toroidal or elliptical contact is, 


1/3 
= 23,600 (= 1/R)28 [3] 
uv 


The half width of contact, or semi-minor axis bo is 
calculated by, 


Po 1/3 
bo = 0.004498v ( ) [4] 
In the finite length roller, the limit of a full ellipse is 
reached when the major axis 2a equals the length of 
contact (ZL). The ends of the ellipse are truncated as 
shown in Fig. 3, when the applied load is sufficient to 
cause the theoretical ellipse to exceed the finite length 


3 The constants in Eqs. [1] and [2] and those following assume 
the elastic modulus E equal to 30 million psi and Poisson’s ratio 
equal to 0.3. 
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Fic. 3. Stress profile over a truncated elliptical contact area 


of contact L possible. Whenever this condition exists it is 
necessary to calculate a new equivalent load based on the 
applied load Py, and the load P;,, which is just sufficient 
to fill the contact length Z with a full ellipse. 

The equivalent load equation (derived in the Appendix) 


is, 


where Py = applied load in pounds limited by P, < Po. 
By replacing Po with P; in Eq. [3], Pz, can be 
obtained, 


= 1.374 10° Z1/R [6] 


P should be used to calculate go, bo at the center of con- 
tact, using the following modified equations of elliptical 
pressure areas: 


Yo = 23,600 (21/R)? [7] 
by = 0.004498v (P/Z1/R)"3 [8] 


Also from the equation of an ellipse, the half width of 
contact 5, and the stress qg, at the end of contact, shown 
in Fig. 3, are: 


b = bo|1 — = bol 1 — [9] 
= qol1 — = go[1 — [10] 


Therefore, from Eq. [1] it is seen that for cylinders in 
contact the assumption is made that the maximum com- 
pressive stress has the same magnitude over the full 
length of contact Z, while for toroidal contact the maxi- 
mum compressive stress go has a single peak value that 
occurs only at the center of the contact area. From 
Eq. [10] it can be determined that the contact stress q, 
at the end of a truncated ellipse will be less than qo 
from Eq. |7], the magnitude depending on the geometry 
of the toroid(s). In reality, however, there may exist a 
stress concentration at the end of the contact so that the 
maximum stress from Eq. [1] or Eq. [7] may actually be 
exceeded by an edge stress g, that would be the critical 
maximum compressive stress for rolling contact fatigue 


considerations. In the next section this stress at end of 
contact is considered. 


Stress concentration at the end of roller contact 


At present the authors are aware of two works dealing 
with the problem of stress concentration in a finite cylinder 
pressed against a longer cylinder. Because of the manner 
in which the shorter cylinder (roller) tends to embed 
itself in the larger body, Lundberg (7) and Nelson (8) 
have worked out formulas to show that the load distribu- 
tion and consequently the deformation is approximately 
uniform over the roller length except near the axial ends 
of the area of contact, and here the local deformation is 
considerably greater. 

In the following, an approximate solution to this 
problem is offered that consists simply of taking into ac- 
count the end radii on the shorter roller. The two roller 
body radii R,; and R» taken from Fig. 1, plus the end 
radius r;’ of the shorter roller are used as shown in Fig. 


R,' 


(c) 
Fic. 4. Geometric approximation of the contact area between 
rollers when one roller is shorter. 
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4a. An ellipse of contact can be constructed from these 
radii in which the major axis is parallel to the rolling 
direction. Then the ellipse shown in Fig. 4b can be 
imagined “split” and each half joined to the ends of a 
rectangular or, as shown in Fig. 4c, a truncated elliptical 
area of contact. 

It is assumed that a’, the semi-major axis of the end 
radius ellipse is equal to 5,, the one-half width of contact 
at the roller end so that the contact stresses reduce to zero 
around the periphery of the combined contact area and 
no discontinuities are present. 

It must be realized that the transition of stresses along 
the major axis AC in Fig. 4c is not rigorously specified 
but at present the main interest is in a first order “stress” 
at points A, B, and C for an approximate solution. 
Following the above assumptions, the stress equation for 
Qe, the maximum compressive stress concentration at 
points A and C for a cylindrical or toroidal roller is 

Ge = [11] 
2x a’b’ 


For a cylindrical roller (loaded against a cylinder), Eq. 
{11] can be rewritten for calculation purposes as, 


__ (5.246 X 10°) by (1/Ri + 1/R2 + 1/n’) 


e 


[12] 


For a toroidal roller forming a truncated elliptical contact 
area, the corresponding stress equation is 


_ (5.246 X 10%) by (1/Ri + + 1/n’) 
v 


e 
[13] 

The increased length of contact that occurs when the 
end radii are taken into account has been ignored, since 
it is negligible and does not measurably change the load 
per linear inch on the roller. The authors have purposely 
considered only the maximum compressive stress to 
simplify the presentation, but this does not mean the 
importance of static principal shearing stresses and 
dynamic range of shearing stresses, as they pertain to the 
rolling contact fatigue problem, can be ignored. The 
rolling contact fatigue rig tests recently reported by 
Greenert (10) indicate a higher order approximation to 
the solution may be obtained considering the critical 
shear stresses present. However, in Greenert’s tests the 
variations in the range of orthogonal shear stresses for 
the three orientations and shapes of the pressure ellipses 
studied were all in the same plane. 

For the problem presented here the critical maximum 
range of shear stress at the center of contact and the criti- 
cal shear stress range at the end of contact in all probabil- 
ity exist on different planes that are oriented differently. 
The findings of the effect of fiber orientation on rolling con- 
tact fatigue life such as Ref. (11) or on torsional fatigue 
life as in Ref. (12) add to the uncertainty as to what 


plane should be associated with the critical shear stress 
range. 


Results of stress calculations 


To present a concrete example of the problem proposed, 
consider a 0.5-inch diameter roller loaded against a 12- 
inch diameter cylinder. Assume the small roller to be 
shorter (Z = 0.93 inch). If the small roller is a cylinder, 
the maximum compressive stress go (from Eq. [1]) will 
have values varying according to load as shown on the 
log-log plot in Fig. 5. 
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Fic. 5. Comparison of maximum compressive stress in the 0.5- 
inch diameter rollers determined two ways: Hertz stress (gp) ; 
end of contact stress (q,). 


Using Eq. [12] the end of contact compressive stress 
de can be computed for the same roller assuming the 
roller end radius equals 0.050 inch. Again Fig. 5 shows 
how q- will vary and be greater than go over the load 
range given. 

For example, under 4000 pounds load, go equals 305,000 
psi for the cylindrical roller, while g, under the same load 
is calculated to be 560,000 psi. If the small roller were 
a toroid with, for example, a 100-inch contour radius 
loaded against the larger 12-inch cylinder, the maximum 
compressive stress go would be significantly higher than 
for the small cylinder. At 4000 pounds load g» would equal 
380,000 psi for the toroid roller (using Eq. [3]). At 
this load a full ellipse of contact is formed so g, equals 
zero. 


The geometry and loads selected above were deliberate 
since one of the rolling contact fatigue rigs in the authors’ 
laboratory has rolling elements of just such geometries. 
The rig is of the “nut-cracker” type with two 12-inch 
diameter cylinders squeezing the toroidal or cylindrical 
test roller between them. The centers of the rolling 
elements are aligned vertically, with the bottom cylinder 
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driving the top cylinder through the rolling contacts on 
the small middle roller. 

The rig has been used extensively for single element 
rolling fatigue tests. As data were collected on toroidal 
and cylindrical roller tests it was apparent that the usual 
Hertz contact stresses were not sufficient to rationalize 
even approximately the fatigue test results. Therefore, 
based on the usual Hertz stress assumptions, a solution 
for the roller end stress was determined so that an analysis 
of fatigue data could be made on the basis of a more 
critical contact stress. 


Experimental evaluation of roller end stress 


An experimental verification of the preceding mathe- 
matical development for computing roller end stress has 
been obtained from a series of deformation studies made 
on rolling elements of several geometries. By measuring 
the small amount of gross plastic deformation occurring 
within the rolling contact area of toroidal rollers, an index 
relating deformation and maximum compressive stress has 
been established. This index has then been used to 
evaluate the maximum stress at end of roller contact of 
cylindrical rollers once the plastic deformation was 
measured. 


Plastic deformation of rolling elements 


Using the rolling rig previously described, deformation 
studies just into the plastic regime were first made on 0.5- 
inch diameter rollers following several revolutions of the 
rolling elements under load. Two cylindrical rollers and 
one 100-inch radius toroidal roller were selected and care- 
fully measured as to profile and contact length. Each 
roller was marked on the end at three points 120° apart, 
as shown in Fig. 6, and using these marks as location 
guides, full length traces of the roller body profile were 
made at each location with the profilograph.* 

The rollers were carburized and hardened 4620 steel 
having a surface hardness of 62 Rc, retained austenite 
content of approximately 21%, and surface finish of 
4-7 winches. The 12-inch diameter loading cylinders were 
equivalent. 

For these deformation tests a new set of ground 12-inch 
diameter cylinders was placed in the rolling contact 
fatigue rig, the cylinders being carefully gauged before 
assembly to insure that their surfaces were true cylinders. 
Each roller was placed in turn in the test rig and a 7500 
pound or 4000 pound load applied to the roller through 
the cylinders. The driving cylinder was rotated by hand 
until each roller rotated about its own center 24 times 
and traveled over only about one-eighth of each loading 
cylinder surface. A new section of loading cylinder sur- 


4 The profilograph is a machine that magnifies and records the 
true surface roughness and contour of finished parts rather than 
giving an average. Vertical magnification is obtained by an optical 
arrangement and in the traces shown vertical magnification is 
365 the horizontal magnification. 
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RECTILINEAR 
PROFILOGRAPH TRACES 


Fic. 6. Location of roller profiles traced in Figs. 8, 9, and 10 


face was used for each roller. The rolling surfaces were 
lubricated with a spindle oil before rotation. After being 
rotated under load, profilograph traces of roller body 
profile were again made over the same locations, as shown 
in Fig. 6. 

To illustrate the plastic deformation that occurred 
during the brief loading described, superimposed copies 
of the profilograph traces of the rollers taken before 
and after loading are shown in Figs. 7, 8, and 9. It may 
be seen that the cylindrical rollers had permanent de- 
formation occurring at the roller end of contact, the 
amount depending on the load. The 100-inch radius 
toroidal roller had center plastic deformation that is just 
discernible from the superimposed traces. It was now 
apparent that since a quantitative measure of the vertical 
plastic displacement of roller material from the original 
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Fic. 7. Profilograph traces of cylindrical roller No. 17. Traces 
taken before and after loading roller to 4000-pound test load. 
Assuming the whole contact surface is permanently deflected 
0.000005 inch (based on Eq. [14]), then the total deformation 
at end of roller contact ranges from 0.000035 to 0.000045 inch. 
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ROLLER POSITION | 
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AFTER sarten” 
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Fic. 8. Profilograph traces of cylindrical roller No. 16. Traces 
taken before and after loading to 7500-pound test load. Assuming 
the whole contact surface is permanently deflected 0.000015 inch 
(based on Eq. [14]), then total end deformation ranged from 
0.000085 to 0.000125 inch. 
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Fic. 9. Profilograph traces of 100-inch radius toroidal roller 
No. 14. Traces taken before and after loading to 4000-pound test 
load. Only the center % portion of the traces shown. 


surface line had been made, it would be very 
beneficial to have some approximate relation between 
plastic deformation and contact stress so that an estimate 
of the compressive stress necessary to have caused the 
permanent deformation could be made. 


Such relations have been developed by other investiga- 
tors. Reference (13) describes the amount of plastic de- 
formation produced when a loaded carbide ball was 
rolled across SAE 52100 steel plates of various hardnesses 
and various amount of retained austenite content. Figure 
10 is summary of data extracted from Figs. 6 and 7 in 
Ref. (13). The plots show depth of the plastic groove 
developed in the flat plates when either a }-inch radius 
or a }-inch radius ball is rolled across each plate once. 
The plastic deformation obtained is plotted versus cal- 
culated maximum compressive stress (go). Three hard- 
nesses and three retained austenite contents from 0 to 
18.4% are represented. The conclusion from Fig. 10 is 
that no understandable relation exists between plastic 
deformation and contact stress unless the important 
parameters of the material are specified. Since the grooves 
from }-inch and 4-inch radii balls can be fairly well 
represented by a single line, it appears that the effect of 
geometry is of secondary importance. Drutowski and 
Mikus (13) state that stress necessary for initiation of 
plastic deformation is independent of ball diameter. It 
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Fic. 10. Depth of groove developed in flat plates versus cal- 
culated contact stress for 3 structures of 52100 steel. Points shown 
are values taken from test curves in Fig. 6, Ref. (13). Hollow 
points indicate tests made with 34-inch radius carbide ball; black 
points 4-inch radius carbide ball. 
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3 EQUATION 14 
9 (1.637 X 107) A ROLLER NO. 17 
~ 40- 2 400 ROLLER NO.14(q)-—@ 
q 
oe 
w 
S 300 
2 3 4 5 6 7 890 
FIRST ROLLER SET APPLIED LOAD Po IN 1000 POUNDS 
: ° @ SECOND ROLLER SET Fic. 12. Comparison of maximum compressive stresses in the ; 
= a 0.5-inch diameter rollers determined three ways: g, (Eqs. [1], [3], 
or (71); g, (Eqs. [12] or [13]); (Eq. [14]). 
pressive stress for the carburized and hardened 4620 steel ( 
; “a SSW ae | considered here, a test was run in a second rolling contact s 
rig. In this test a 2;4;-inch diameter toroidal roller> with 
8 8 2S & g 2 a }-inch radius profile was rotated once under successively c 
higher loads against a 14-inch diameter cylinder, _using 2 
j (CALCULATE D) the same track. The depth of the plastic groove in the 7 
Fic. 11. Depth of plastic groove (A) versus calculated maxi- yang was measured after each loading. F igure 11 I 
mum compressive stress. 14-inch diameter cylinder loaded against Saws & log-log plot of the measured deformations — . 
a 27-inch diameter 3-inch radius toroid (carburized and hardened calculated maximum compressive stress (qo). A straight I 
4620 rollers). line was drawn by eye through the data points similar to 
Fig. 10; this then provided the approximate relationship 
i also seems reasonable to assume that the depth of the needed for the 4620 roller material: 
3 lastic groove is approximately independent of geometr 
P PP 8 y q = (1.637 X 107) [14] 


if the elastic deformations are considerably greater than 
the plastic deformations encountered. a 


‘“ To determine a similar first order approximation of the 5 The physical properties (material, hardness, retained austenite, c 
, relation between plastic deformation and calculated com- etc.) were equivalent to the 0.5-inch diameter rollers. r 

TABLE 1 e 
Comparison of Contact Stresses le 

3 Roller geometry Stress calculation Max. comp. stress Cc 
(load) Method Equation Symbol Magnitude® 
: Cylinder (4000#) Hertz [1] %% 305 ir 
E End of contact {12] GV 560 ti 
; Plastic deformation [14] q 530-590 
i (measured range) 

Cylinder (7500#) Hertz [1] 415 
End of contact [12] Ve 770 = 
Plastic deformation [14] 720-820 
(measured range) 
ea 100-inch radius toroid (4000#) Hertz [3] WG 380 
Plastic deformation [14] qd 355--400 
(measured range) c 


® Thousands of psi. 1¢ 
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A repeat test on a second. set of rollers at the highest 
stress level indicated that the depth of plastic deformation 
obtained on the cylinder of the first roller set had not 
been appreciably affected by the over-rollings. The de- 
formation obtained on the second cylinder is shown in 
Fig. 11. Table 1 shows the comparison of stress calcula- 
tions using Eqs. [1] or [3], [12] and [14]. It is felt 
the agreement between the maximum compressive stresses 
qg- and q is sufficient to justify this presentation. Figure 12 
shows this agreement graphically. 


Discussion 


A first order solution for the critical maximum com- 
pressive stress present in rolling contact for the variety 
of geometries presently used in rolling contact rig testing 
has been provided. Confirmation based on deformation 
studies just into the plastic region established that the 
end of contact stress for the conditions specified can be 
approximated. 

A number of rolling contact fatigue tests have been 
completed on the 0.5-inch diameter roller discussed in this 
paper. The bulk of this data on a variety of toroidal 
geometry specimen rollers has been presented in Ref. 
(14). However, to illustrate the success of the first order 
solution for rationalizing cylindrical and toroidal rolling 
contact fatigue tests, the fatigue data for a group of 
cylindrical rollers tested under 3500 pound load and a 
group of 100-inch radius toroidal rollers tested under a 
7500 pound load are plotted on Weibull paper® in Fig. 13. 
Lines passed through the data points as shown enable 
estimates to be made of Weibull By) (10% failure) and 
Bso(median life) life values. 

Table 2 compares calculated go and g, contact stresses 
for the respective rollers tested and their corresponding 
Byo and Bso (median life) life values. 

When these specific life estimates are examined in the 
light of the usual Hertz stresses (qo), the fatigue results 
are inexplicable. However, when the life estimates are 
compared with the end of contact stresses (g,) the 
results are clearly what one would expect. Actually 
whether go or g- is the critical stress to be used in life 
evaluation depends on which is higher for the particular 
loads and geometry studied. Table 2 provides additional 
confirmation that a suitable, although approximate, end 


6 The use of the Weibull distribution is so universally accepted 
in rolling contact fatigue investigations that it is felt no explana- 
tion is necessary. See, for example, Ref. (15). 
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SPECIMEN LIFE IN MILLIONS OF STRESS CYCLES 


Fic. i3. Rolling contact fatigue life data on cylindrical and 
100-inch radius toroidal rollers. 


of contact stress has been determined for the cylindrical 
roller examined. 

Other tests considering other conditions (various 
toroids, roller end radii, materials) are in progress in 
order to determine how far the concepts presented may be 
generalized. 

The influence of the lubricant film on contact stresses 
suggested by the recent work of several investigators (16, 
17, or 18) has not been considered. According to Ref. 
(17) the lubricant film thickness is affected more by 
surface velocity than by pressure so that nominal stresses 
calculated, even go (Hertz stress) following classical 
theory, must be only approximate under rolling contact 
conditions. The role of cumulative deformation has not 
been discussed, but if contact stresses change as geometry 
changes during rolling contact until some “shakedown 
limit” (19) is reached then a plasticity correction as sug- 
gested by Moyar (20) or a relation between plastic 
deformation and over-rollings (21) will also be needed to 
determine more closely the contact stresses under actual 
operating conditions. 
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TABLE 2 
Contact Stress and Fatigue Life Comparison for the Cylindrical and 100-inch Radius Toroidal Rollers 


Cale. max. Weibull life 
Calc comp. stress (psi) estimates (10 
Load contact Qo (Hertz q, (end of stress cycles) 
Roller (Ibs) area stress) cont. stress) Byo Bs 
Cylinder 3500 Rectangle 285,000 525,000 41 235 
100-inch R toroid 7500 Truncated ellipse 480,000 490,000 67 262 
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Appendix 
TRANSCENDENTAL FUNCTIONS 


\t and v are transcendental functions useful in calculat- 
ing parameters in elastic contact stress problems. An 
explanation of the functions and their values can be 
found in Ref. (9). Approximate values of these func- 
tions, as they relate to toroidal radii of curvature, are 
shown in Fig. 14. The general relation to toroidal radii 
is by the geometry factor, 


+ 1/Re— (1/Ry’ + 1/R2’) 
1/R, + 1/R2+ 1/Ry’ + 1/R2’ 
For the ellipse in Fig. 4b, that considers the end radii 
of the shorter roller, the geometry factor is, 
+ — 1/ry’ 
1/R; + 1/R2 + 
The values of this geometry factor may be minus simply 
meaning the major axis of the ellipse is in the rolling 


direction. However, only magnitudes are important in 
using Fig. 14. 


DERIVATION OF THE EQUIVALENT Loap EQuATION FOR 
A TRUNCATED ELLIPSE 
When a toroid and toroid, or toroid and cylinder are 
loaded against each other the total load is considered dis- 
tributed over the contact surface in the shape of a semi- 
ellipsoid, as shown in Fig. 15. The sum of the pressures 
over the contact surface can be expressed as 


(ao/bo) (bo? — 
Po =4Be [1 — (x?/ao”) — (y?/bo”) ]** dx dy 
15 
where 


Po = applied load 
4o,69,Co = semi-axes of the ellipsoid 
E = the elastic modulus 
x,y = axes of the pressure area 
2 = Co[1 — (x?/ao?) — (97/07) = height 
from the pressure surface to the ellipsoid sur- 
face, as shown in Fig. 15. 


Integrating Eq. [15] with respect to x and y gives, 
NECobo — ] "1 


dy" 3 0 


= 1 —% (= )} 
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Fic. 14. Plot of geometry factor versus transcendental functions, 
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z When 
= / (when / < ay) then Eq. [16] becomes, 
Po = xbolqo [1 — (1/ao)?] [18] 


Fic. 15. Semi-ellipsoid of pressure 


where 


qo = Ecy = maximum compressive stress 
¥1 = a limit to be defined below. 


When 
¥1 = do Eq. [16] simplifies to, 
Po = (2% aobogo)/3 or Go = (3Po)/(2x aobo) [17] 


From elliptical contact theory it is known that the rela- 
tion between load and semi-minor and major axes for 
the same geometry is 


= (@;/az) = (Pi/P2)** [19] 


When a, = /, then P,; equals P;, the load necessary to 
form an ellipse just sufficient to fill the actual roller 
contact length 2/ = L. If the applied load Py» is greater 
than P,, a truncated ellipse will be formed; the theoretical 
semi-minor and major axes for this ellipse will be a2 = 
a, b2 = be and an equivalent theoretical load P2 — P 
can be determined. Equation [19] becomes, 


a (P,/P)*3 [20] 


Substituting Eq. [20] for 7 and 3P/2xa9bo for go in Eq. 
[18], the equivalent load equation is obtained: 


P is equal to an equivalent load for a full ellipse that will 
give the same maximum compressive stress go as the 
applied load Py produces for the truncated ellipse. 


DISCUSSION 


P. F. Brown (Pratt & Whitney Aircraft Corp., Hartford, Connecti- 
cut): 


The authors have made a comprehensive and practical contribu- 
tion to furthering the understanding of the stress distribution 
picture associated with cylindrical and toroidal roller contact. 
Their very careful analysis and the implications of the results of 
the analysis are substantiated to a large degree by their appropriate 
deformation studies and the correlation of their results with the 
rolling contact fatigue test data obtained on both cylindrical 
and 100-inch toroidal specimen rollers. Studies of this type are 
becoming increasingly necessary as rolling element bearings load 
and environmental conditions become more severe, and long life 
is no longer a luxury but a serious requirement for applications 
such as advanced aircraft and space vehicles. The authors should 
be extended a vote of thanks for their addition to the state of 
the art for which the goal is to arrive at a modification of the 
life formula that will precisely predict performance. 

Although the expressions for stress concentration on the ends 
of cylindrical and toroidal rollers are only a first order of magni- 
tude approximation, as admitted by the authors, it is both in- 
teresting and encouraging that such a simplified approach can 
be correlated with actual fatigue test results obtained from rolling 
element fatigue testers. It would be interesting to hear whether 
the authors have had any experience in actual applications in- 
volving full scale bearings or have conducted full scale bearing 
fatigue tests that would add further testimony to the influence 
of edge loading on fatigue life and whether these results correlate 
with the findings presented in their paper. 


A type of roller geometry that has not been considered by the 
authors and which is of prime interest to many roller bearing 
makers and users is the “crowned” roller design. This type of 
profile is a combination of the cylindrical and toroidal contour 
with the cylindrical section occupying from 30 to 60% of the 
roller length. It is realized that the authors could not consider 
every possible profile in their paper, and this is the obvious ex- 
planation for the omission. In light of this, I would ask the au- 
thors if they are conducting or plan to conduct a “crowned” 
roller study similar to that presented in this paper. If such a 
study were made and it were presented in the same precise fashion 
as this paper has been, it would receive perhaps even more ac- 
claim. 


L. B. Sistey (SKF Industries, Inc., Research Laboratory, King 
of Prussia, Pennsylvania) : 


The problem of stress concentration at the ends of the roller 
race contacts in roller bearings is an important practical design 
consideration that has been studied both experimentally and theo- 
retically for a number of years. The guiding philosophy in such 
investigations has always been to reduce the level of stress at the 
ends of roller-race contacts enough that failures no longer concen- 
trate in these regions but are distributed evenly over the width 
of the raceway in the bearing. 

The authors have presented a highly approximated and probably 
misleading method for calculating the magnitude of the stress 
concentration at the ends of the contact between a cylindrical 
roller of finite length loaded between two cylindrical races whose 
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width is greater than the length of the roller. Essentially, they 
have assumed that the major axis of the half contact ellipse 
formed at the corner radius of the roller is equal in magnitude to 
the width of the contact on the cylindrical portion of the roller 
under load, assuming infinitely long cylinder theory according to 
Hertz, and have thus calculated the maximum compressive stress 
so determined at the roller end, again using Hertzian theory, 
with the further assumption that the cylindrical portion of the 
roller is completely removed and the two end corners are joined 
into a body of toroidal shape. In the single numerical example 
given a roller end radius of 0.050 inch is assumed on a 0.5-inch 
diameter roller loaded against a 12-inch diameter cylinder thus 
producing, according to the above method, a maximum compres- 
sive stress at the roller ends some 80% greater than the maximum 
compressive stress on the cylindrical portion of the roller. This 
method of calculating roller end stresses is approximate, at best, 
and is probably misleading for the range of roller dimensions and 
end radii assumed by the authors, for the following reason: when 
the roller end corners are joined for the calculation of the end 
stresses, it is tacitly assumed that the support provided by half 
of the toroid so formed is similar to the support provided by the 
cylindrical portion of the roller which, of course, produces an 
error, making the stresses so calculated too high by an amount 
which increases as the roller end radius decreases. Conversely, by 
taking the width of the contact at the roller end from Hertz 
theory for infinitely long cylinders, it is tacitly assumed that the 
support provided by the half toroid of the roller end radius is 
similar to the support that would be provided if the roller were 
continued to infinite length which, of course, produces an error 
causing the stress calculated to be too low. It would be highly 
fortuitous and very unlikely that these two effects cancel and 
thus, the end stresses so calculated are probably in error, especially 
for small roller corner radii (incidentally, in practical roller mak- 
ing, the roller corner radius rarely intersects the finished roller OD 
on a tangent due to final grinding operations, thus accentuating 
the error discussed above, by making the actual corner radius 
even smaller). 

The authors go on to discuss a toroid roller with a profile radius 
so calculated, apparently to provide a maximum compressive 
stress at the center of the roller which is just as great as the end 
stresses but somewhat less than the end stresses calculated for the 
cylindrical roller discussed above. The comparison of the rolling 
contact fatigue life of these two roller types shown in the au- 
thors’ Fig. 13 indicates that practically the same life characteristics 
are obtained. There was an insufficient number of bearings tested 
to determine with any confidence that there was any difference 
at all between the life of these two groups of rollers. Figure 13 
does show that no great improvement in life was realized in going 
to the toroidal roller, certainly not as much life as would be ex- 
pected from the calculated stresses indicated in Fig. 12. 

It should be pointed out that this problem of roller end stresses 
was studied extensively by Lundberg in 1939 (the authors’ Ref. 
(7), who derived the theoretical profile required on finite length 
cylindrical rollers in order to produce a constant maximum com- 
pressive stress right out to the end of the roller contact. This 
theoretical profile, incidentally, is not anywhere near toroidal in 
shape. More recently, Kunert (Kunert, K., “Spannungsverteilung 
im Halbraum be; elliptischer Flaechenpressungsverteilung ueber 
einer rechteckigen Druckflaeche,” Forsch. Gebiete Ingenieurw. 27, 
165-174 (1961)) has derived the sub-surface stresses for the 
Lundberg profile and has found that even though the surface 
maximum compressive stress for this profile is the same on the 
ends as it is on the center of the roller, the subsurface shear 
stresses which are apparently critical in rolling-contact fatigue 


failure are substantially higher at the roller ends than they are 
at the center of the roller. Thus, some modification to the Lund- 
berg profile is required for optimum roller endurance character- 
istics. 


GeorcE S. REICHENBACH (Associate Professor of Mechanical Engi- 
neering, M.I.T., Cambridge, Massachusetts) : 


The authors have presented a simple approximate method for 
estimating stresses at edge discontinuities in Hertz-type contacts. 
The correlation between these calculations and the inception of 
plastic flow is encouraging and, in view of the reasonable re- 
producibility of yield data, is probably well confirmed enough 
for engineering purposes. 

The correlation with fatigue life, however, is not very well 
supported by the limited data presented. The scatter in rolling 
contact fatigue data can often be extreme, and one should be very 
careful in drawing conclusions from Weibull plots containing fewer 
than twenty or thirty points. In this writer’s experience repeat 
sets of test containing only eight points, as the plots of Fig. 13, 
will often show substantial changes in B-10 and B-50 lives. A 
variation of more than a factor of two is not unusual. When 
viewed in this light, the reported change of 41 108 to 67 « 108 
cycles does not seem very significant. 

One can use the statistical techniques of Johnson (C1) to place 
confidence limits on the data presented. For a Weibull slope of 1.0 
and a sample of 10, approximately the same values as the data 
reported, the following confidence limits result: 


Cylinder 7.8 X 108 < B-10 life << 213 x 106 —90% confidence 
Toroid 12.7 X 108 < B-10 life < 349 x 106 —90% confidence 


These are very wide limits indicating that changes as small as 
30-50% in life are not very significant when so few test points 
are involved. Similar calculations of confidence limits on the mean 
life show a four to one spread possible at the 90% level. 

This is not meant to indicate that the approach presented is of 
no use but rather that for fatigue purposes more extensive ex- 
perimental correlation is required. 


REFERENCE 


C1. Jounson, L. G., “Ball Bearing Statistical Guide Book,” Gen- 
eral Motors Research Staff, New Departure Division, General 
Motors Corporation, Bristol, Connecticut (no date). 


AUTHORS’ CLOSURE 


The authors express appreciation to the discussers for their 
comments. In reply to Mr. Brown, similar work has been done 
on full-scale tapered roller bearings; however, not all of the 
techniques described are directly applicable. Because of this, exact 
correlation of results cannot be expected. 

Work has also been done in the authors’ laboratory on 
“crowned” or cylindrical-toroidal combination rollers. Actually the 
“crowned” roller can be produced by making the roller end 
radius r,’ larger than is shown in Fig. 4a in the paper. How far 
such “contouring” can profitably be explored is illustrated in Figs. 
7 and 8. Since the plastic deformation shown extends in from the 
ends of contact roughly 0.040 inch or less than 10% of the contact 
length, it would serve no useful purpose to contour the ends of 
the cylindrical roller any more than necessary to relieve the end 
stresses. Equations [1], [2], and [12] can be used to determine an 
approximate “end radius” geometry Mr. Brown alludes to in 
“crowned” roller design. 

Mr. Sibley has reviewed the approximations in the stress 
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equations, as previously acknowledged by the authors, very com- 
mendably. Since we have not attempted to generalize beyond the 
roller geometries discussed in the paper, we do not see how the 
method of approach presented can be misleading. 

Professor Reichenbach’s remarks concerning the statistical nature 
of fatigue are well taken; however, significant statements, 
especially at the Bs, level can be made with small test samples. 
The point of Table 2 is that assuming Hertz stresses (gg) are 
valid, then with the toroid roller B;, life of 262 million stress 
cycles, one can calculate an expected B;, life for the cylindrical 
roller using a stress-life exponent of 20/3 or Bs, (cylinder) = 
(262 X 106) (285,000 psi/480,000 psi)—2°/3 stress cycles = 
8450 X 108 stress cycles (expected assuming Hertz or gq, valid). 


Making the same calculation with q, stresses from Table 2 one 
obtains: Bs, (cylinder) = (262 X 106) (525,000 psi/490,000 
psi) —2/3 stress cycles = 165 X 108 stress cycles (expected as- 
suming q,). 

As given in Table 2, B;, for the cylindrical roller test sample 
actually equalled 235 x 10® stress cycles, which certainly agrees 
with the expected value of 165 x 106 stress cycles using g, as 
compared to the expected value of 8450 x 108 stress cycles as- 
suming g,. Therefore, even considering the confidence limits 
Professor Reichenbach suggests, the conclusion remains that 
meaningful correlation of the rolling contact fatigue data can be 
realized only when critical maximum compressive stresses (for this 
example q,) are compared among various tests. 
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Theoretical Investigation of Externally Pressurized 4 
Thrust Collar Gas-Bearings i 


By HARUO MORI’ and HIROSHI YABE! 


Externally pressurized thrust collar gas-bearings with several supply holes are investigated theoret- 


ically. The gas supply holes are spaced circumferentially and communicate directly with the film 


without feeding into recesses. 


A complex potential function which satisfies the boundary conditions of the bearing is found, and 
the pressure distribution, quantity of flow, and load capacity are calculated for the two cases 


where there are or are not radial grooves between the holes. 
The theory is applicable for any number of gas supply holes. 


Nomenclature 


A, A’, B,, = functions defined in text 
a — radius to the supply hole from 
the bearing center 
h = film thickness 
i = unit of imaginary number 
j= number of supply holes ar- 
ranged in radial direction 
k = number of supply holes ar- 
ranged in circumferential di- 
rection 
Ky, Ky. Ki, Ke, 
= constants 
m = integer 
n = polytropic exponent 
P= patel» 
= pressure, absolute 
= pressure at gas supply hole 
= ambient pressure 
Q = volume rate of flow 


q = intensity of source 
R = radius in Z-plane 
r — radius in €-plane 


r, — radius of circular thrust bear- 
ing 

r, — radius of supply hole 

U = sliding velocity 

w — complex potential function 
x,y = orthogonal coordinates 

Z = complex number 

© — argument in Z-plane 


Contributed by the ASLE Technical Committee on Bearings 
and Bearing Lubrication and presented at the Annual Meeting of 
the American Society of Lubrication Engineers held in New York, 
April, 1963. 

1 Department of Mechanical Engineering, Kyoto University, 
Kyoto, Japan. 


= argument in C-plane 
¢ = complex number 
ut = absolute viscosity 
= function 


SUBSCRIPTS 

inner edge of thrust collar 
outer edge of thrust collar 
supply 


0 
1 


Introduction 


It has been reported that the externally pressurized gas 
bearing with large recesses has good load capacity but 
poor stability (7). It may be better to use a thrust collar 
bearing having gas supply holes which are arranged 
circumferentially and communicate directly with the film, 
thereby avoiding recesses. 

Accordingly, in this paper the thrust collar bearing 
without recess is investigated by using the complex 
potential theory in order to obtain pressure distribution, 
quantity of flow and load capacity. 


Pressure distribution 
CoMPLEX POTENTIAL 


Using the generally employed assumptions of the : 
hydrodynamic theory for the gas lubricated bearing, 3 
Reynolds equation is as follows for the compressible, - 
viscous, laminar film without inertia effects: 


(1+n)/n (1+n) /n 


Ox Ox oy \ Oy 
l+n O(hp'") 
= ow ( ) | + 4 ] [1] 
n Ox oy 


the gas flow is assumed to be polytropic. 
When the sliding velocity U is small in comparison with 
the gas flow velocity induced by film pressure gradients, 


ad 

: 

‘ 

if 
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the right-hand side of Eq. [1] can be neglected. If U = 0, 
the following Laplacian equation is obtained as the film 
thickness / is taken to be constant in this type of thrust 
bearing: 


=0 [2] 
ox" oy 
where 
P = [3] 


For an incompressible fluid, P = p so that the actual 
pressure p satisfies Eq. [2] (2). 

The problem of solving Eq. [2] can be transformed 
into the problem of finding a suitable complex potential 
in the plane of complex number Z = x + iy. If a com- 
plex potential is obtained, the real part of the potential 
function gives the pressure distribution for this case and 
the imaginary part gives the stream function. The pres- 
sure distribution must satisfy the boundary conditions of 
the actual bearing. 


Consider an externally pressurized thrust collar gas- 
bearing as shown in Fig. 1, and take a conformal trans- 
formation given by 


where & is the number of the gas supply holes in Fig. 1. 
A point on the Z-plane defined by Z = Re'® corresponds 
now to a point on the C-plane 


C= = 


The hatched region of the Z-plane in Fig. 1 is transformed 
to the whole of (-plane in Fig. 2. Thus, if the thrust 
collar bearing with only one gas supply hole as shown in 
Fig. 2 is analyzed, other bearings with several gas supply 
holes can be investigated by using the above transforma- 
tion. 


Now consider the following analytic function in the 
{-plane; 


Fic. 1. Configuration of thrust collar bearing on Z-plane 


Fic. 2. Thrust collar bearing on ¢-plane 


q { | 


— In sin eS In —— 
ro" 


Substituting Co (= roe'@), which represents the inside 


-circle of the bearing edge, the analytic function vanishes 


its real part and only has an imaginary part as follows: 


q 


Wo = Wo [6] 
2x 
where 
mln (7o/a) 
| 6a | 
2In(79/r1) 
+ 2mx 
0, +1, +2,..... ) [6b] 
2In (70/11) 


as shown in Appendix 1. 

Similarly substituting €; (= r,e'@), which designates 
the outside circle of the bearing edge, into Eq. [5], 
w becomes 


w, = w(rie) = i(g/2x) [7] 
where 
] 
A, | 7a] 


as shown in Appendix 2. 


Thus the real parts of wo and w,; become zero so that 
the analytic function given by Eq. [5] can satisfy the 
boundary conditions at both edges of the bearing, and 
gives a suitable complex potential function for the bear- 
ing. 

Dividing w into two parts of w; and wy;;, and substitut- 
ing C = re” into Eq. [5] to obtain the real part of w 
which gives the pressure distribution, 


the 


R( 


wh 


th 


n= 
| 
wl 
fli 
Supply Hole 
/ 
7 
k 
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= — _ In — 
2x 2In(10/r1) a 


| + 2mzx) 
= sin —————— cosh 
2x 21n(r0/r1) 
an(r/a) 2mz) 
+ icos inh [8] 
2In(ro/r1) 2In(ro/r1) 
the real part of wy is 
R(w;) = In | sin?A cosh*B,, + cos*A | 
where 
r 
A = —— ln — 
2In(ro/r1) a 
Now 


q ac 


the real part of w;; becomes 


q 


R(wy) = an 


In [sin?A’ + sinh*B,, | [11] 


where 


ar 


, 


= ——— ln — 
2In(ro/r1) 


Therefore, the pressure distribution for an incompressible 
fluid can be expressed as 


p= Ky, [In(sin?A’ + sinh*B,,) — In(sin?A 
4+ sinh?B,,)] + Ke [12] 


where K, and K» are constants to be determined by the 
following boundary conditions: 


at r—a—r, 


at Orr, [13] 


where r, is the radius of the gas supply hole. 


Putting the boundary conditions into Eq. [12], Ki 
and K» are given by 


Ps— Pa 


Ko = 


In[ (sin?A,’ + sinh*B,,) /(sin*A, + sinh?B,,) | 


Ks = Pa [14] 


where 
a—r, 
A,= In 
2In(ro/r1) a 
a(a—T,) 


, 
8 


— n 
2In(r0/r1) fo 


It should be noted that the boundary conditions are 
satisfied only at one point of the edge of the gas supply 
hole. However, if € is near to a in Eq. [5], the curve of 
constant pressure is given by a circular form so that 
this pressure distribution may satisfy the condition of the 
gas supply hole, if 7, is small enough compared to the 
bearing dimensions. 

In the summation of Eq. [12], since the terms of 
m = 0 are negligibly small in comparison with the term 
of m = 0, the pressure distribution can be expressed by 


p = [In (sin*A’ + sinh*Bo) 
— In (sin?A + sinh?By)| + Kz [15] 
Ps— Pa 


In (sin?A,’/sin?A,) 


CASE OF A COMPRESSIBLE FLUID 


The solution for a compressible fluid can be easily 
obtained by transformation through Eq. [3]. If the letter 
P is substituted for p of the preceding solution, P satisfies 
Eq. [2] and gives the pressure distribution in the form 
of p("+1)/" where p is the real pressure for the com- 
pressible fluid case. But K, and Ky must be changed to 
satisfy the boundary conditions for the real pressure. Con- 
sidering the relation between p and P in Eq. [3], p, and 
pa in the K, and Ky must be changed to P, and P,, 
respectively, through the relations 


n 
= n 


where ~, and p, correspond to the boundary pressures for 
the compressible fluid case. 

Thus the solution for a compressible fluid can be ex- 
pressed by 


PCTS [16] 


where P is the corresponding solution obtained for the in- 
compressible fluid case using the correct boundary con- 
ditions. 


MULTIPLE Gas SupPLY HOLES 


When there are several gas supply holes in each sector, 
which are equally spaced circumferentially, as shown in 
Fig. 3, the pressure distribution is obtained by super- 
posing the solution which corresponds to the case of one 
supply hole in each sector. 


If there are & arrays of supply holes circumferentially 


i 

Al 

re: 

A" 4 

x 

* 


Fic. 3. Configuration of the bearing with multiple gas supply 
holes. 


and j supply holes in each array in radial direction, the 
pressure distribution is given by 


j k 


m=— 


+ —In (n+ } 


1 


+ sinh*B,’) — In (sin?A, + sinh*Bg’}] + 


where 
r 
2In(70/r1) dq 
, 1,2 
B 6+2 1) 6 
= + 2mx— — 1) 6x] 
B,’ 6 (p’ 1) 6 
[@ — (p’—1) 4x] 
2n 
(angular interval) 


The boundary conditions to determine the constants 
Ki(a,g) and are 


at 


P= Pa at 
(a=1,2,...,j; B=1,2,...,2) 


6= (B— 
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Fic. 4. Thrust collar bearing with radial grooves between 
supply holes. 


Sink 


Fic. 5. Arrangement of sources and sinks in the case of 6, = 
60°. 


BEARING WITH RADIAL GROOVES 


If the bearing has radial grooves between the gas 
supply holes, as shown in Fig. 4, the pressure distribution 
on one sector can be obtained by taking alternate sources 
and sinks spaced 6; as shown in Fig. 5. 

The angle 6; can be taken as the angle of the real 
sector in a particular bearing, so that & may be an 
arbitrary real number. 

The pressure distribution obtained similarly as Eq. [17] 
is 


P= ¢a,8) [In (sin?A,’ 


a=1 
+ sinh?Bg’) — In (sin?A, + sinh*Bg’)]} + 
[18] 
Equation [18] is valid for a sector bounded by rp < r < 


r, and — (6%)/2 < @ < (6,)/2. The boundary con- 
ditions to determine the constants K1,4,8) and Ko,q,g) are 


(Q 


(Q 
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P= Pa at 


Load capacity 


The load capacity can be obtained by integrating the 
pressure distribution over the whole bearing surface. It 
should be noted, however, that the pressure at the supply 
hole remains constant at p,. 


Quantity of flow 


The quantity of flow can be calculated from the out- 
flow from the bearing edges. 


For the incompressible fluid, it becomes 


{ Op 
necomp — r 
(Q); 6u er dé 
ikl 9 
0 or 
i 
<i [19] 
3u a=1 
in the case of the bearing without radial grooves, 
(QO) meomp = 
kh’ Op ‘Kl Op 
{ f —— | r=r9 dé r dé 
6u 0 or 0 or 


j 
a=1 


[20] 


1 OP 


for the case of the bearing with radial grooves. 

In the neighborhood of the supply hole, the pressure 
can be expressed as follows for the incompressible fluid, 
when the supply hole located at the origin 


1 
p = 2K, ln—-+ Ke 
r 


For the case of the compressible fluid, the pressure 
distribution in the neighborhood of the supply hole is 
similarly obtained and if isothermal flow (m = 1) is 
assumed we get 


°—2K,In—+ Ke 
r 


where 
(sin2A,’/sin2A,) + Pe) 
= Ef 


Therefore, the volume rate of compressible flow from 
one supply hole can be given by 


OP K, 
(Qs) comp r = 
12 0 or 6u Ps 
s+ Pa 
6u Ps 
Ps + Pa 
= (Q) incomp x 


Since volume rate of flow (Qs) comp corresponds to the 
quantity at the pressure p,, it can be transformed to the 
volume rate of flow corresponding to the atmospheric 
pressure Pq by Qa = Qs X (Ps/Pa). 

Thus, the volume rate of compressible flow cor- 
responding to the pressure p, can be given by 


Ps + Pa 
(Qa) comp = (Q)incomp X [21] 


It is obvious from Eq. [21] that one can get the 
volume rate of flow for a compressible fluid from the 
volume rate of flow for an incompressible fluid which is 
given by Eq. [19] or Eq. [20]. 

Results of calculation for some examples 

The bearing dimensions chosen for the calculations 
were as follows: 

Lal / %o = 3.0 
a/ % = 2D 


0.04 

ps = 1 kg/cm? 

The number of supply holes is taken to be one, as 
shown in Fig. 2, and six, as shown in Figs. 1 or 4. The 


case of one supply hole is fundamental since the case 
of multiple supply holes is obtained by superposition of it 


= 
—— Incompressible Fluid 
—-— Compressible Fluid 
0.8 ...... Conventional Circular 
an Thrust Bear 
| 
r/Yo= 
04 
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Fic. 6. Pressure distribution in the circumferential direction 


for the case of one supply hole. 
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Figure 6 shows the pressure distribution calculated on 
the circle of radii with r/ro = 1.5, 2.0, and 2.5; and 
Fig. 7 shows the pressure distribution on the radii locating 
on @ = 0°, 5°, and 30° for the case of one supply hole. 

One of the practical methods of estimating the load 
capacity of this type of bearing has been to substitute for 
: sectors a conventional circular pad which has its center 
os at the gas supply hole, as shown in Fig. 8. The load capa- 
+ city for each circular pad is summed to give the load capa- 
city of the thrust collar. 

The pressure distribution for this conventional substi- 
tution is also shown by a broken line in these figures. 
There is considerable difference in the pressure distribu- 
tion in the circumferential direction so that the con- 
ventional calculation of load capacity can yield an under- 
estimation. 

Figures 9 and 10 show the pressure distribution for 
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Thrust Bear j 
1.0 ! 
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Fic. 7. Pressure distribution in the radial direction for the case 
of one supply hole. 
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the case of six supply holes spaced equally in the cir- 
cumferential direction. 

The load capacity obtained by numerical integration 
and volume rate of flow calculated for the same given 


Conventional 
Circular 


Fic. 8. Conventional substitution of circular thrust bearing for 
thrust collar bearing. 
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Fic. 9. Pressure distribution in the circumferential direction 
for the case of six supply holes. 


TABLE 1 
Load Capacity and Volume Rate of Flow 
Incompressible fluid Compressible fluid 
No. Vol. rate of Vol. rate of 
supply Radial Load flow at p, Load flow at p, 
holes grooves (kg) (cm3/sec) (kg) (cm3/sec) 
1 no 0.930 0.151 h3/po 1.236 0.227 h3/p 
6 no 0.885 0.144 h3/p 1.119 0.220 h3/p for 
6 exist 0.577 0.158 h3/p 0.743 0.237 h3/p each 
Conventional circular pad 
thrust bearing 0.478 0.163 h3/p 0.622 0.244 h3/u 


@ hin cm, p in kg sec/cm?. 
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——— Incompressible Fluid | 

Compressible Fluid | 
Conventional 

10 Compressible) _ 


p 


Pressure, 


Fic. 10. Pressure distribution in the radial direction for the 
case of six supply holes. 


conditions, in which rp = 10 mm, a = 20 mm, 7; = 30 
mm, 7, = 0.4 mm and p, = 1 kg/cm’, are shown in 
Table 1. The radial grooves reduce the load capacity and 
increase the quantity of flow. The conventional circular 
thrust bearing has lower load capacity than true one, but 
it might be a good approximation for the case with radial 
grooves. 

Concerning the number of supply holes, the bearing 
seems to be better having at least five or six without radial 
grooves in order to keep the pressure between them 
higher. 


Summary 


In this study, the authors have investigated the ex- 
ternally pressurized thrust collar gas bearing having 
supply holes in both the radial and the circumferential 
direction and the effect of radial grooves. An appropriate 
complex potential function has been found to yield the 
pressure distribution for each case. The quantity of flow 
and load capacity are then calculated from the pressure 
distribution. 

Experimental investigations may be necessary to assure 
the theoretical results. The authors have obtained solu- 
tions using an electrical analogy which coincide with the 
theoretical results given here. 
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Appendix 1 


Logarithmic function of complex number is calculated 
as follows: 


2 2 cos—1 
In (a + ib) = $1n (a? 4+ +i cos 
[22a] 
or 
In (re®?) = Inr+i(0+2mx) m=0,+1,+2,... 
|22b] 
Using Eq. |22b], 


Wo = w(roe'*) 


q | ain (ro/a) + 2mz) 
= —In |} sin cos. 
2In (70/r1) 2In (70/71) 
ain (ro/a) + 2mzx) ] 
+ i cos sin 
2In (70/r1) 2In (r0/r1) 
q | mln (a/ro) + 2mz) 
——lIn | sin cos. 
2x 2In (70/11) 2In 
alin (a/ro) + 
+ icos sin 
2In (70/71) 2In (r0/r1) 
q sinAg coshB,, + i cosAg sinhB,, 


sin (— Ao) coshB,, + i cos (— Ao) sinhB,, 
q (— sin*Ay cosh*B,, + cos*Ao sinh?B,, ) 
sin*A, cosh*B,, + cos*Ao sinh*B,, 
2i sinAy cosAg sinhB,, coshB,, 
sin*Ag cosh*B,, + cos*Ao sinh*B,, 


where Ay and B,, are defined by Eqs. [6a] and [6b] 
in text. 


Using Eq. [22a], the real part of the function wo is 
—-sin*Ay cosh?B,, + cos*A9 sinh*B,, \ 
4x (= cosh*B,, cos*Ao sinh*B,, ) 
( 2sinAy cosAy sinhB,, coshB,, | 


sin?A, cosh?B,, + cos*Ao sinh*B,, 


q 
=—In = 
4x (1) 


Thus wo has only an imaginary part determined by Ao 
and By. 


| l 0° 

| | 
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Appendix 2 
Similarly 


w= 
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_ (r:/a) 

= sin 
2In (ro/r1) 

an (ar;/ro”) 
2In (70/r1) 


ain (7;/a) 
os 
2In(ro/r1) 


q _ (r:/a) + 2mzx) 
(10/11) 2In (ro/r1) Using Eqs. [22a], [7a], and the above, one obtains 

nln (r;/a) + 2mz) | q coshB,, + i cosA, sinhB,, 
Zin (70/r1) 2In (10/11) sinA, coshB,, — i cosA, sinhB,, 

q ain + 2mzx) q sin?A, cosh*B,, — cos*A, sinh*B,, 
—— In} sin cos. = —In— - 

2x 2In (70/71) 2In (70/11) 2x sin?A, cosh?B,, + cos7A, sinh*B,, 
mn (ar;/ro*) 4+ ] 2i cosA, sinA, coshB,, sinhB,, 

2In(r0/r1) 2In (r0/r1) sin?A, cosh?B,, + cos?A, sinh?B,, 

mln (ar;/ro") xin (a/r;) ] q y 

sin ——————— = sin | — =—i— Bm 
2In (ro/r1) 2In(7o/r1) 2x 
DISCUSSION 


J.S. AusMan (Litton Industries, Woodland Hills, California) : 


The authors present a very elegant solution applicable to several 
types of thrust bearings lubricated either with liquids or with 
gases. I have two minor points to raise, one a question to the 
authors and the other a clarification to the reader applying the 
author’s results to a practical problem. 

First, what is the magnitude of the error introduced by 
truncating all terms in the series in Eq. [12] except for the 
m = 0 term, and over what range of parameters is this approxima- 
tion acceptable? 

Second, lest the reader be lulled into complacency by the seeming 
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Fic. Al. Technique for determining p, 


simplicity of the authors’ results, it should be pointed out that in 
most practical applications, p, is not really a constant but 
differs from the actual supply or reservoir pressure p, by the 
pressure drop across the inlet orifice and the piping system 
supplying the inlet orifice. This pressure drop and hence p, can 
be determined as indicated in Fig. Al. The flow rate Q is plotted 
as a function of p, for various clearances h using the authors’ 
Eqs. [19], [20], or [21] to produce the dashed lines shown in 
Fig. Al. The Q vs. p, characteristic for the orifice and its piping 
supply system can be plotted on this same graph. Intersections 
of this orifice characteristic with the dashed lines give possible 
operating values of p, as a function of the clearance h. 

It is this variation of p, with clearance h which provides the 
necessary bearing behavior of an increasing load with decreasing 
clearance. Note that if , was in fact held constant, the pressure 
distribution given by Eq. [12] would be independent of the 
clearance h, and the bearing would have zero stiffness. 


W. A. Gross (Ampex Corp., 2385 Bay Road, Redwood City, 
California): 


This paper is an extension of a previous article co-authored by 
one of the authors (B1). The configuration analyzed by the 
authors have been frequently used, and design engineers should find 
the equations for pressure, load, and flow to be of value. Since 
this paper does not take into account external restrictors, results 
may be used only if the inlet and outlet pressure at the lubricating 
film are known. In order to obtain an approximate idea of the 
characteristics of bearings such as shown in Figs. 1-4, the designer 
may find it convenient to use this paper in conjunction with a 
design paper (B2) which does take external restrictors into ac- 
count. 

Balcerzak and Raynor (B3) also solved the Laplace equation 
for a variety of gas bearing shapes. Their studies, which included 
bearing shapes with square corners, were obtained for incom- 
pressible lubrication films. Their results, therefore, may be applied 
only directly to the pressure in a gas bearing film. Gas flow, 
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stiffness, and load in such bearings can be approximated by the 
present paper and reference (B3). 
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The authors wish to thank Dr. Ausman and Dr. Gross for their 
valuable discussions. 

For the first question of Dr. Ausman, the magnitude of the 
error is very small when the ratios of r,/ry and r,/a are not so 
large as is common for thrust collar bearings. For example, the 
values of the pressures at 6 = x and 2x are shown in Table C1. 
These are for an incompressible fluid with one supply hole located 
at (a, 0). 

TaBLe Cl 
p at 0 = wand 2z in kg/cm? when r,/r, = 3, a/ry = 2, 
r,/@ = 0.02 and p, = 1 kg/cm? 


o=<x O= 2x 
15 6.1 10-5 7.6 X 10-9 
2.0 6.1 x 10-5 7.6 X 10-9 
25 33 «x 4.1 xX 10-9 


The convergence in Eq. [12] is very rapid, almost geometric. 
Therefore the values in Table C1 can be considered to show the 
order of the error from which it may be concluded that it is less 
than 0.01% when all terms are truncated except the term of 
m = 0 in Eq. [12]. 

But it increases to a certain extent with increase of r,/ro, r,/a@ 
and number of supply holes when the other conditions are not 
changed. These variations are shown in Table C2 and Table C3 
for several values of r,/rg and r,/a respectively for an incom- 
pressible fluid with one supply hole; these show the increase of 


the order of the maximum error when we ignore the terms except 
the term of m = 0 in Eq. [12]. 


TaBLe C2 
p at 6 = xandr = ain kg/cm? when r,/a = 0.02 
and p, = 1 kg/cm? 


1 


a= (ror,)? @= 
3 0.00007 0.00006 
0.00110 0.00078 
10 0.00643 0.00345 
TaBLe C3 


pat @ = xandr = ain kg/cm? when r,/r, = 3, 
a/r, = 2 and p, = 1 kg/cm? 


r,/a= 0.01 5.1 X 10-5 
0.02 6.1 X 10-5 
0.03 6.9 x 10-5 
0.05 8.2 x 10-5 


The variation of the error with respect to the number of the 
supply holes can be estimated by the superposition of the case with 
single supply hole. 

In any case, we may have to take several terms in Eq. [12] 
for a precise calculation of the pressure distribution when the 
error is presumed to be fairly large. 

As to the second discussion by Dr. Ausman, he is quite right 
and the authors thank him for his reference. This is the same case 
that Mori (2) has analyzed in which fp» corresponds to p, in 
this paper. So one should take the variation of p,/p, into con- 
sideration as Mori has done. But the authors did not mention 
this in the paper because they intend to obtain the pressure 
distribution as an ideal case providing that the pressure p, is 
known. On a practical problem, one must adopt the method of 
Dr. Ausman. 

One may have to take into account the variation of the pres- 
sure at the inlet hole in practical design of bearing. It is 
noticed that the inlet pressure is determined by the bearing load 
and does not concern the bearing clearance. Our results may be 
useful to obtain gas flow, bearing clearance and stiffness in con- 
nection with the external restrictor as Dr. Gross mentioned. 

The authors wish to thank Dr. Gross for his comments on 
external restrictors and for his reference of the application to the 
different bearing shapes. 


+ 
d 
ex 
- 
Shy 


ASLE TRANSACTIONS 6, 346-357 (1963) 


Theoretical and Experimental Investigation of 
Gas-Lubricated, Pivoted-Pad Journal Bearings 


By E. J. GUNTER, JR.,! V. CASTELLI and D. D. FULLER (astE)* 


This paper presents some of the theoretical predictions and experimental results for the steady- 
state characteristics of a gas-lubricated, pivoted-pad journal bearing of finite length. An analytical 
expression for the prediction of load-carrying capacity for these bearings is developed from 
numerical computer solutions. Sample calculations are presented to show how the theoretical data 
may be utilized in the design of an actual pivoted-pad bearing. Improvements over previous ap- 


proximate theories are discussed. 


Nomenclature 


B= shoe length in the circumferential direc- 
tion = Ra, inches 

C = radial clearance between shoe and shaft 
= Ranoe — Renate; inches 

C, = load coefficient = W/(P,RL), dimension- 
less 

C’ = pivot circle radial clearance = Ryivor — 
Rgnatt, inches 

e= eccentricity between shoe and _ shaft 

center, inches 


e’ = eccentricity between pivot circle center 
and shaft center, inches 
fa = side leakage factor, dimensionless 


h» = minimum film thickness, inches 
= pivot film thickness for shoe i, i = 
1,2,3, inches 
h, = leading edge film thickness, inches 
hy = trailing edge film thickness, inches 
H = dimensionless film thickness = h/c 
L = shoe width, inches 
L/B= aspect ratio, dimensionless 
= rpm 
P,, = ambient pressure, psi 
p = local lubricant pressure, psi 
P = p/P, = dimensionless pressure 
P, = ?’th pivot location; i = 1,2,3 
R = shaft radius 
= shoe load capacity, pounds 


Contributed by the ASLE Technical Committee on Bearings 
and Bearing Lubrication and presented at the Annual Meeting of 
the American Society of Lubrication Engineers held in New 
York, April, 1963. 

1 Research Engineer, Friction and Lubrication Laboratory, The 
Franklin Institute Laboratories, Philadelphia 3, Pennsylvania. 

2 Assistant Professor, Department of Mechanical Engineering, 
Columbia University, New York. 

3 Professor, Department of Mechanical Engineering, Columbia 
University, New York. 
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W; = load on ?’th shoe; i = 1,2,3, pounds 
W, = shaft load 
O, = pivot circle center 
O; = journal center 
x = horizontal displacement of the journal 
center from the pivot circle center 
y = vertical displacement of the journal center 
from the pivot circle center, inches 
\n,, |My = orthogonal cartesian unit vectors 
n; = unit vector set directed from O, to P; ; 


|r; = unit vector set directed from O; to P; ; 
423 


a = shoe arc length, degrees 

6B = pivot location, degrees 

€ = eccentricity ratio of shoe and shaft = 
e/C, dimensionless 

e’ = eccentricity ratio of pivot circle and shaft 
e’/C’, dimensionless 

1 = dimensionless axial coordinate 

= bearing compressibility parameter, dimen- 
sionless 


\t = viscosity, Ib sec/in* 
— = angle between line of centers and shoe 


leading edge, degrees 

= angle between shoe leading edge and pivot 
location, degrees 

= dimensionless pivot location 

6’ = attitude angle of journal center with 

respect to pivot circle center, degrees 
6 = angular coordinate, degrees 
® = angular velocity, rad/sec 


Introduction 


IN recent years considerable interest has been devoted to 
the development of design criteria for gas-lubricated bear- 
ings appropriate for high-speed rotating machinery. The 
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numerous advantages of gas bearings are widely known 
(1, 2) and need not be discussed here. While many of 
the associated design problems have been successfully 
solved or circumvented, the extreme susceptibility of gas 
bearings to the phenomenon of self-excited half-frequency 
whirl instability under the conditions of high speed and 
light load has limited the application of the more 
conventional bearing configurations. On the basis of ex- 
perimental (3) and approximate analytical results (4), 
the tilting-pad bearing configuration has demonstrated 
favorable stability characteristics. This feature, together 
with the desirable characteristic of self-alignment, sug- 
gested that detailed investigation of this type of bearing 
might prove to be very worthwhile. 

Until now, the best available design technique for 
tilting-pad gas bearings was to treat the bearing as a flat 
inclined slider. Such a technique was the basis of the 
early analysis performed on pivoted-pad gas bearings 
by Snell (3), who utilized the work of Muskat e¢ al. (5), 
and of Raimondi and Boyd (6). This technique consists 
of replacing the sinusoidal film distribution due to shoe 
curvature, with a linear function and then considering 
that the lubricating film is incompressible. This leads to 
considerable error in the prediction of the minimum film 
thickness and the selection of the optimum pivot location 
at high speeds and high loads due to the neglect of com- 
pressibility and shoe curvature. 

A computer program for the tilting-pad bearing has 
been developed by Castelli and Stevenson (7) which 
includes shoe curvature and compressibility of the gaseous 
lubricant. The present paper will compare the results of 
the computer solution to those obtained by using the flat 
slider and the assumption of an incompressible lubricant. 
Experimental verification of theoretical predictions for 
the static-load case is provided by data obtained from a 
pivoted-pad research rig developed and operated at the 
authors’ laboratories (8a, 6, and c). 


Bearing geometry 


A schematic representation of a pivoted-pad journal 
bearing with three pads is shown in Fig. 1. For the 
purpose of analysis, the pivot points of each pad are taken 
to be located on the pad surface. These pivot points are 
also located on a circle, called the “pivot circle,” of radius 
R + C’, where R is the radius of the shaft and C’ is the 
radial clearance of the pivot circle. 

The pivot circle clearance C’ is an important design 
parameter and may be expressed in terms of the film 
thicknesses under each pivot by means of the following 
relationship: 


+ Apo + cos 
2(1 + cos B) 


The shaft center eccentricity measured with respect to 
the pivot circle center is: 


ef = (a? + y*)% [2] 


[1] 


22 
22 
2 
— 


JOURNAL BEARING 
Fic. 1. Schematic diagram of a pivoted-pad bearing con- 
figuration with three shoes. 


where 


— hye 
[3] 
2 sin B 
and 


y= [4] 
2(1 + cos B) 
The shaft attitude angle with respect to the pivot 
circle and the load line is: 


= [5] 
x 
The geometry of an individual shoe and shaft is shown 
in Fig. 2. To describe the film distribution of a given 
partial-arc bearing, it is necessary to specify three 
quantities. These quantities consist of the clearance 
parameter C, which is assumed as fixed and the two 
variables of eccentricity ratio « and lead angle E, which 
may assume various values. In the analysis of a composite 
pivoted-pad bearing configuration, the total load-carrying 
capacity of the bearing is obtained by computing the 
performance of each shoe and combining the individual 
results to determine the net bearing capacity. The 
variables of ¢ and & were used in the computer program 
to analyze the individual shoe characteristics. 
For a multi-pad bearing, there would be a set of ¢ and 
— which would correspond to each shoe. As an example, 
for a three-pad bearing, the film thickness distribution 
between the 7’th shoe and shaft is given by 


PIVOT CIRCLE 

Za SZ a 
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Fic. 2. Individual pivoted-pad bearing and shaft geometry 


h, = C[1 + & cos 4]; i = 1,2,3 [6] 


where C is the difference between the radii of the shoe 
and the shaft surfaces as manufactured, and ¢; is the ec- 
centricity ratio = e;/C. 

To analyze the total bearing performance, the pivot 
film thickness of each pad must be specified. In this sys- 
tem the pivotal film thickness is given by 


hy = C[1 + ecos (E+ ¢)];#=1,2,3 [7] 


Theoretical solution 


The steady-state pressure distribution in the bearing 
is assumed to be governed by the two-dimensional 
Reynolds equation for an isothermal film (2) which can 
be expressed as 


P R? P 


00 00 on on 
O(PH 
086 
where 
6um R\? 
Acs bearing compressibility parameter 


yn = Z/L = dimensionless axial coordinate 
P = p/P, = dimensionless pressure 
H = h/C = dimensionless film thickness 


The boundary conditions require that the dimension- 
less pressure assumes a value of unity at the periphery of 
the bearing as follows: 


Ey) = P(O= + = 
P(6y=+4)=1 


Equations [8] and [9] were expressed in finite dif- 
ference form and a computer program was developed for 


[9] 


their solution by an iterative technique (7). This program 
was used to produce field maps (to be discussed later) 
determining the operating characteristics for a partial 
bearing with an arc of 94.5° and a shoe aspect ratio of 
1. 

From Fig. 1 it is evident that, if the pads are sym- 
metrically located with respect to the load line, they will 
each have a maximum arc length equal to 26. The 
optimum value of the angle 6 and consequently of the 
angular extent of the bearing pads is obtained as a com- 
promise between load-carrying capacity and resulting 
load on an individual shoe or pad. 

The average load per bottom shoe is given by 


W3+W, 
was 2 cos B [10] 

From these considerations an optimum value of 6 
equal to 50° was chosen and the corresponding maximum 
allowable arc length of 100° was reduced to 94.5° for 
practical design reasons. 

Figure 3 represents a typical computer field map for 
the partial bearing under consideration for a particular 
value of the parameter A = 3.5. (This corresponds to a 
shaft rotational speed of approximately 18,000 rpm, which 
is the design point for the experimental unit.) 

For a given value of the compressibility bearing 
parameter A, a computer field map is obtained by 
selecting numerous combinations of shoe eccentricity 
ratio e and lead angle &. For each set of values of ¢ and E 
there is obtained a corresponding dimensionless load and 
the position at which it acts. The location of the resultant 
load vector will determine the pivot position for static 
equilibrium. Therefore, for a pivoted-pad shoe, the 
performance characteristics are determined by the cor- 
responding values located along the vertical line where 
¢/a equals a constant. 

An important factor in the design of a tilting-pad bear- 
ing is the selection of the optimum pivot location. This is 
defined as the position for which the maximum load-carry- 
ing capacity is obtained for a prescribed minimum film 
thickness. 

Lines of constant minimum film thickness can be 
plotted directly on the computer field map (Fig. 3). For 
example, the line H,, = 0.40 represents a minimum film 
thickness Hy», = A»/C = 0.40. To the right of the line 
—E = 85.5°, the fluid film is converging, and the minimum 
thickness occurs at the shoe trailing edge; to the left, the 
film is converging-diverging and the minimum thickness 
occurs before the trailing edge. 

From the study of flat sliders with incompressible 
lubrication, it is seen that the optimum pivot location is 
independent of load, but varies with the aspect ratio of 
the shoe ranging from 0.60B for an aspect ratio of unity 
to 0.579B for an infinitely wide bearing (6) (Fig. 4). 
From the computer field map of Fig. 3 for the com- 
pressible case, it can be seen that the optimum pivot loca- 
tion for a shoe with aspect ratio of 1.0 is not constant 
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Fic. 3. Theoretical performance of a partial journal bearing for A = 3.5 
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0.7 08 0.9 10 
PIVOT POSITION FACTOR, z (DIM.) 


Fic. 4. Theoretical performance of a flat slider with incompressible lubricant 


but varies with load and speed and ranges from 0.62 to 5 

0.64. With this partial bearing using a compressible ost 

lubricant, an increase in load causes the optimum pivot _ OT} 

position to move closer to the trailing edge. However, for z 06 

very light shoe loads the effect of compressibility becomes pales 

negligible and the optimum pivot position approaches se 

0.60B, the value predicted for bearings with incom- CONDIT 

Since the choice of an optimum pivot location is 5 03+ 
dependent on both load and speed, it should be selected to = Guu (Ry2 

correspond to the design conditions. Examination of the i Pa 
slope of the minimum film thickness curves reveals, how- 8 on i a=94.5° 
ever, that the magnitude of the minimum film thickness t $.2 
is fortunately not sensitive to small deviations in pivot ? 
location. 

Results from computer field maps for values of A 0. | | ar Vie a a 
ranging from 1.5 to 4.0 and ¢/a = *% were used to “Ol 02 03 04 05 06 O08 10 
develop the theoretical load deflection curves of Fig. 5. MINIMUM FILM THICKNESS, Hy= (DIM) 

(The pivot location of the experimental shoes was also at Fic. 5. Theoretical load coefficient C,, versus minimum film 


¢/a = %%.) The linearity of these curves on a logarithmic thickness for various values of A. 
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plot indicated that an equation could be fitted to them. 
Such an equation was developed relating the basic 
performance parameters as shown in the following form: 


2 (R/L, ¢/a, A 
Ww =— (R/L, ¢/a, A) 
[11 


Equation [11] makes possible a direct comparison with 
the analytical results for flat sliders using incompressible 
lubricants. The incompressible case is equivalent to 
taking m = 2 in Eq. [11] and the side leakage factor f, 
is a function only of the aspect ratio and pivot location. 

The factor f, was used by Raimondi and Boyd (6), 
and Muskat e¢ al. (5). For the test conditions of speed 
and load encountered with the authors’ gas-lubricated, 
tilting pad shoes, the f, factor is a constant which is not 
speed dependent. This value was computed to be about 
0.235 for ¢/a = %4, as compared to the values of f, = 
0.245 for the flat slider with an aspect ratio L/B = 1 as 
given by Raimondi and Boyd (6). Thus, it can be seen 
that one effect of compressibility is to slightly reduce the 
fa factor with a corresponding reduction in load as com- 
pared to the incompressible case. 

The other important variable of the above equation is 
the exponent “n” of the film thickness term. The value of 
“n” is slightly less than 2 for the compressible case and 
decreases slowly with increasing A. The value of “n” 
may be given by the expression: 


n = 2.08 — 0.052A [12] 


for A varying from 1.5 to 4.0. 

Thus, the two theories predict that, for the same 
minimum film thickness, the compressible film will carry 
less load than the incompressible film with a flat slider. 


Sample calculation 
EXAMPLE I 


An an illustration of the use of Eq. [11], the load- 
carrying capacity for a single pivoted pad will be com- 
puted from knowledge of the speed and minimum film 
thickness. The results will then be checked against the 
computer field map to test the accuracy of Eq. [11]. 


Assuming the following input data: 


R = 2.00 inches 
C = 0.0015 inch 
B= 3.3 inches 
L = 3.3 inches 
P, = 14.7 psia 


N = 17,650 rpm 
u = 2.61 reyns 
Re, = 0.0007 inch 
a = 94.5° 
¢/a= % 


Then 
= 0.235 
A= 3.5 


n = 2.08 — 0.052 (3.5) = 1.898 
Equation (11) yields 
14.7 ( Cc (3.5) (3.3)?(0.235)? 
6 0.606 
C \ 1.898 
= 36 lb 


(It is understood that a design pad load could have been 
specified for which the minimum operating film thickness 
h,, could have been calculated instead of specifying a film 
thickness and computing a load as in the example above. ) 


To check the results of Eq. [11] through the use of the 
computer field maps we can compute the load coefficient 
C_, which equals (W)/(P,RL), to be 


Cy = 36/(14.7) (2) (3.3) = 0.371 


From this load coefficient and ¢/a — 74 the computer 
field map of Fig. 3 yields ¢ = 0.59 and —E = 60°. 

The corresponding leading (4,,) and trailing edge (7) 
film thicknesses are 


h;, = cos = 0.00194 [13] 
hy = C[1+ cos (E+ a)] = 0.000701 [14] 


Since E is smaller than 85.5° the film thickness occurs at 
the trailing edge. The calculated value of A, is in ex- 
tremely good agreement with the one specified as a design 
parameter for this problem, thus confirming the good fit 
of Eq. [11] to the computer results. 


Design of complete bearing 


Figure 6 is obtained as a crossplot of the computer field 
map and is extremely useful in the design of complete 
bearings. 

One of the crucial parameters governing the operation 
of a complete bearing is the setting of the pivot circle 
clearance c’. This is caused by the necessity for maintain- 
ing a stable hydrodynamic film under the top or preload 
shoe at very low speeds. With low speeds, the shaft will 
operate at very high eccentricity ratios with respect to the 
pivot circle setting, and the pivot film thickness of the top 
shoe may become larger than the machined-in radial shoe 
clearance. 


If the top shoe should vibrate or be disturbed in any 
way from its equilibrium position when operating under 
the above conditions, the hydrodynamic restoring moment 
becomes negative, causing the shoe to rotate about its 
pivot until the leading edge of the shoe touches the shaft. 
In this “locked-in” position, the shoe forms a diverging 
fluid film clearance area which in turn creates a low pres- 
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PIVOT FILM THICKNESS, Hp (DIM.) 


PIVOT LOCATION, £ (DIM.) 


Fic. 6. Pivot film thickness versus pivot location for A = 3.5 


sure region keeping the shoe in contact with the shaft 
(see Fig. 7). 

To avoid this phenomenon, two recourses are available: 
one is to stabilize the shoe by the application of external 
moments and thus physically raise the leading edge and 
maintain the normal shoe attitude. It is this attitude that 
creates the load-carrying converging clearance space in 
the bearing. Second, by a proper choice of c’/C such that 
at the lowest predicted operating speed, a stable con- 
verging clearance is maintained under the top shoe. The 
proper choice of c’/C is dictated by the requirement that 
Ing C. 


EXAMPLE II 


As an illustration, a complete three-shoe bearing could 
be designed in the following manner: 


Input data: Same individual shoe geometry as previous 
example. 
Shaft weight = 86.5 lb 
50° 
A= 3.5 
Selecting a dimensionless pivot film thickness for the 
top shoe such that A,3/C < 1.0 or assume H,3; = 0.9. 
Figure 6 yields ¢ = 0.65 and E = 36°. 


Entering the computer field map with these values, the 
top shoe load coefficient is determined to be C3 = 0.25. 


Then the load coefficient of the bottom shoes is: 


] 57 
= 0.574 
2cosBL 
[15] 
The computer field map (Fig. 3) gives the operating 


068 70 o72 O74 O76 


LOW PRESSURE 
REGION 


Fic. 7. Unstable shoe condition when (a) shoe position before 
“lock-up” and (b) shoe position at “Lock-Up.” 
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Fic. 8. Test rig with self-acting bearings, SA-2 


conditions of the bottom pads as €;,2 = 0.62 and &1. = 

From Fig. 6, Hy1,.2 = 0.53. 

The resulting pivot circle clearance is obtained from 
Eq. [1] rearranged as 
Hy: + Apo + 2H ps cos 


2[1 + cos 


0.53 + 0.9 X 0.643 
= 0.675 [16] 
1.643 


Therefore 
C’ = 0.675 X C = 1.012 X 10-3 inches 


or the pivot circle radius should be approximately one mil 
larger than the shaft. The machined-in clearance C is of 
course 14 mils larger than the shaft. 


Experimental investigation 


The experimental rig is shown in Fig. 8. It consists of 
a plain cylindrical rotor driven by an air impulse turbine 
and supported by two symmetrically located three-pad 
bearings. The unit is provided with two externally pres- 
surized lifters to independently support the shaft at low 
rotational speeds. At approximately 10,000 rpm the air 
supply to the lifters is shut off so that the total shaft 
load is carried hydrodynamically by the tilting-pad bear- 
ings. Capacitance probe mounts are provided so that the 
motion of the shaft and all shoes may be monitored with 
respect to ground. The film thickness is then obtained by 
differences. One shoe is also equipped with probes for the 
direct measurement of leading and trailing edge film thick- 
nesses (Fig. 9). A close up of the bearing-shaft assembly 
without the top shoe is shown in Fig. 10. 

Testing was performed without the top shoe in each 


bearing so that the bearing loads were due to shaft weight 
only and were thus accurately known. Under the above 
mentioned conditions the operating leading- and trailing- 
edge film thicknesses were measured and are shown in 
Fig. 11 along with a comparison of the theoretical predic- 
tions by the exact computer solution (8) and also Snell’s 
approximation (3). 

It will be remembered that Snell treated the pad surface 
as equivalent to a flat inclined slider and the fluid as be- 
ing incompressible. As expected the latter assumption 
causes his solution to diverge from the exact one as the 
value of A increases. The deviation caused by the neglect 
of the sinusoidal variation of the clearance distribution 
(shoe curvature) will be considerable at high eccentrici- 
ties. 

The agreement between the computed and the experi- 
mental results is very satisfactory. Part of the small devia- 


Fic. 9. SA-2 shoe with small probes attached 
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tion could be ascribed to squeeze film effects that are 
caused by the departure of the shaft surface from a per- 
fect cylinder. Another source of error could be the un- 
certainty in the value of C, the machined-in shoe clear- 


ance, used in converting from ¢ and € to 4, and Ay. The 
original value of C = 0.0015 inch, employed in all cal- 
culations, is now probably slightly higher due to repolish- 
ing of the shaft surface. 


Fic. 10. SA-2 shoe installed on shaft 
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Fic. 11. Shoe film thickness measurements versus theoretical predictions 


| 
i 

| 

| 

| 


354 


It was considered desirable to test the tilting-pad bear- 
ing at various pivot locations in order to experimentally 
ascertain the optimum pivot position. 

In order to shift the effective pivot location, external 
moments were applied to the existing shoes. Whereas the 
resultant hydrodynamic force developed remains relatively 
unchanged, the location of the resultant load vector shifts 
from the physical pivot position to counteract the ap- 
plied moment. This then is equivalent to having the shoe 
pivoted at a new position corresponding to the intersec- 
tion of the shifted load vector with the bearing surface. 


Figure 12 shows the effect of external moments on the 
leading and trailing edge film thicknesses and the cor- 
responding effective pivot location. It can be seen that 
shifting the pivot location beyond the two thirds point 
toward the trailing edge of the pad causes a slight reduc- 
tion in the trailing edge film thickness and a relatively 
large increase in the leading edge film thickness. In order 
to carry the equivalent load at increasing values of ¢/a, 
the shoe inclination must increase, thus demonstrating that 
the optimum pivot location must be equal to or less than 
two thirds for the A condition of test. This is in agree- 
ment with the conclusions drawn by examining the com- 
puter field map, Fig. 3, which predicts values ranging 
from 0.62 to 0.65 depending upon speed and load. This 
trend is the same as for the flat slider bearing using an 
incompressible lubricant. 
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Fic. 12. Pivoted-pad film thicknesses versus A for various ap- 
parent pivot locations. 
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Conclusions 


The good correlation obtained with the experimental re- 
sults confirms the validity of the numerical computer 
solutions of Ref. (7). The usefulness of the field maps in 
the design of pivoted-pad gas bearings is enhanced by the 
establishment of a closed form equation which is ap- 
plicable for the prediction of load-carrying capacity over 
a moderate range of the compressibility parameter A. The 
optimum pivot location has been determined to be in the 
range 0.62 to 0.65 for the range of A’s investigated in 
contrast to the value of 0.5 previously proposed in Ref. 
(3). 
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Appendix 
DERIVATION OF Pivot CIRCLE CLEARANCE 


For a three-pad bearing configuration the shoe pivot 
points form a unique circle referred to as the pivot circle. 
An important design parameter is the value of C’, the 
pivot circle clearance setting. To completely specify the 
orientation of the journal with respect to the pivot circle 
three quantities must be known. These quantities are: 


(a) R+ C’ = pivot circle radius 
(b) e’ = pivot circle eccentricity 
(c) 6 = attitude angle 


The above quantities will be expressed in terms of the 
shoe pivot film thicknesses /,;. 


Irs 
Ing 
Ps CIRCLE 
hp. 
4 


Fic. 13. Pivot circle and shaft geometry 


From Fig. 13 it can be seen that the position vectors 
from the pivot circle center O, to the shoe pivot points 
P; (i = 1, 2, 3) can be expressed as 


OpPPi — (R |m = 1,2,3 [17] 
where 


|n; = unit vector set directed from O, 


The above position vectors can also be expressed as 
OppPi Oppo; O;pP; 
+R, 


where 


|r; = unit vector set directed from O; 
123 

e’ = eccentricity vector = x |n, + y |n, 

|n,, |\ny = orthogonal cartesian unit vectors 


Equating the two sets of position vectors 


(R+C’) =x + y |my + (R + 
[19] 
The equations of transformation between the |; unit 


vector set and the orthogonal cartesian unit vectors is 
given by 


—sin 6 cos \n 
| = sin cos [20] 
Nz 0 —1 


Take the dot product of Eq. [19] by the unit vectors 
im 


[((R+C’) =x \ny + (R + - in 
[21] 

fa 

if e’/R << 1 then 


|r, + = |ro = |mg ~1 
The above results in three scalar equations which are 


(a)R+C’=-xsinbB+ yoosB +4,,+ R;i=1 
(b)R+C’=xsinB+ yeoosB +42.+R;i=2 
(c) R+C’ = — y+ R51 =3 [22] 


Solving for C’ 


+ + 2 cos B Ays 123] 
2[1 + cos B] 


Solving for y 
2hys (Api hy2) 


yr=ky—C= | 24] 
[22.c] 2[1 + cos B] 
Solving for x 
—h 
[25] 
2 sin B 


Solving for the attitude angle 6’ of the journal center 
with respect to the pivot circle 


6’ = tan-! — 
x 


= tan~| ( sin )] 
1 + cos — 


[26] 
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DISCUSSION 


A. A. Rarmonpi (Westinghouse Research Lab, Pittsburgh 335, 
Pennsylvania) : 


A pivoted pad journal bearing with only three pads would 
ostensibly have an undesirable characteristic, particularly for a 
gas bearing. The variation in eccentricity with the direction of 
the load (between the pads as compared say to “on” the pivot) 
is appreciable resulting possibly in large changes in stiffness. This 
effect is mitigated as the number of pads is increased (see Ref. 
A1). It would seem that at least four and preferably six or more 
should be used. Viewed in this light, the effect of pad curvature is 
lessened as the number of pads increases (since the pad arcs 
become considerably shorter) so that flat pad theory should give 
better agreement. The agreement for the present case appears to 
the writer to be tolerable and might have been better had com- 
pressibility been accounted for. 

The treatments of pivoted pad journal bearings available 
presently are not necessarily restricted to flat pad approximations. 
Reference (A2) points out that curvature is important and ac- 
counts for curvature by using a crowned pad technique. 

The subject of lock-up has already been treated in some detail 
and brought to attention in a recent paper (A2). There it is 
called “spragging” and means of controlling it shown. The top 
(unloaded) pad of a three pad bearing would be almost certain 
to sprag. 

Spragging can be a real problem in gas bearings and we have 
found that when the unloaded pad does sprag, it can interfere with 
the proper operation of the loaded pad(s). 
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P. C. WarNnER (Mechanical Technology, Inc., Latham, N. Y.): 


The authors’ interesting paper and invaluable data will be 
particularly appreciated by those designers of high-speed rotating 
equipment who are faced with the most difficult gas-lubricated 
journal bearing design problems. In view of the strong assets of 
the tilting pad type of bearing and its widespread use where con- 
ventional lubricants are employed, it is difficult to understand 
the almost complete lack of theoretical treatment accorded it. 
Fortunately, the authors’ paper goes a long way toward filling 
the obvious gap. 

Inasmuch as the chief reason for employing this pad type of 
bearing is to eliminate fluid film whirl, it would appear that perhaps 
less attention should be spent on the precise details of minimum 
film thickness and other similar quantities and more on that 
property which makes a properly designed pad bearing unable to 
support fluid film whirl. In Hagg’s analysis (4), which visualized 
incompressible lubricants, and perhaps is restricted to reasonably 
small A values, the salient point would appear to this discussor 
to be the statement that, “Because the bearing pad is pivoted 
and free to rotate, not only is the direction of P from pivot 
point to journal center, but the magnitude of P is a single valued 
function of the distance separating the pivot center and journal 
center.” Certainly if pad inertia is so great that the pads cannot 
follow the vibrating shaft, then the pads are not free to rotate, 
and the force need not go through the pad pivot point; hence, 


the stability properties of the bearing must degenerate to those of 
fixed partial arc bearings. While adequate film thickness at each 
pad is certainly a necessary condition for proper bearing design, 
it is not sufficient in most cases, and several design decisions 
should take cognizance of the extreme desirability of low inertia 
pads and high angular stiffness of the lubricant film. It is felt 
by this discussor that it is very desirable to keep the pitching and 
rolling natural frequencies of all the pads at least twice as high 
as operating speed. Such an accomplishment assures proper ability 
of the pads to track the journal, hence freedom from fluid film 
whirl. 

Further, the fluid film properties required to estimate the pitch- 
ing and rolling frequencies can be calculated by numerical dif- 
ferentiation of the computer map, Fig. 3, provided that the 
frequency of motion involved isn’t too high. 


AuTHorsS’ CLOSURE: 


Mr. Raimondi, in his discussion, points out that a three-pad 
bearing arrangement would have an appreciable variation in 
bearing stiffness in transverse directions. This is true and has been 
experimentally observed. One aspect of the experimental phase of 
the pivoted-pad gas bearing investigation has been to obtain quasi- 
static bearing stiffness values at various speeds. These values have 
been obtained from load-deflection data obtained by a foil bearing 
device which preloads the shaft by various amounts. The dif- 
ference in bearing stiffness characteristics has manifested itself 
in the behavior of the rotor rigid body resonance frequencies 
(“critical speeds”). Several resonance frequencies are monitored 
on the oscilloscope in which the rotor motion is predominantly a 
horizontal mode at low speed and changes to essentially a vertical 
mode at higher speeds. All resonance frequencies encountered 
occur well below the design speed of 18,000 rpm for the unit, in 
the 5000-10,000 rpm range. No difficulty has ever been en- 
countered in passing through these resonance speeds, and as 
such the authors feel that unsymmetrical bearing stiffness in 
itself will not preclude the successful operation of a three-pad 
gas-lubricated bearing arrangement. 

Of paramount importance, however, to the successful operation 
of a gas bearing rotor of this size, is the high degree of rotor bal- 
ancing that must be attained. Since the viscosity of air is very low, 
approximately 1/1000 that of oil, the damping characteristic of a 
gas-lubricated bearing is substantially lower than that obtained 
with a conventional oil-lubricated bearing. As such, in order to 
maintain permissible amplitudes of motion while passing through 
the rotor resonance frequencies, the rotor must be extremely well 
balanced. Considerable difficulty was encountered in trying to 
obtain the necessary rotor balance by commercial techniques. A 
method was finally devised whereby the rotor was accurately and 
rapidly balanced in place by monitoring the rotor motion by 
means of the capacitance probes (8b, 9). 

It was felt that the three-pad bearing arrangement is preferable 
because this configuration yields the maximum load-carrying ca- 
pacity and also for the consideration of aligning the shoes. For a 
three-pad arrangement the three pivot points will determine a 
unique pivot circle setting. If more than three shoes were em- 
ployed, it was felt that some difficulty might be encountered in 
trying to align the shoes so that the pivots would fall on a 
common pivot circle. 

The authors wish to iterate here that the analysis does include 
the effect of compressibility of the lubricant as well as shoe curva- 
ture since the developed computer program takes into considera- 
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tion the exact film distribution between the shaft and the shoe. 
Hence Eq. [11] presented in the text includes the effects of shoe 
curvature as well as lubricant compressibility. The film distribu- 
tion corresponding to the above case is not linear (shoe arc 
length = 94.5°) but is a slowly varying sinusoidal function. As 
Mr. Raimondi mentions, the use of more shoes (reduction in shoe 
arc length) will reduce pad curvature and hence the shoes will 
be better approximated by flat slider theory. 

The authors appreciate the viewpoint taken by Mr. Warner in 
that it would be desirable to have more information regarding 
the stability characteristics of the pivoted-pad bearing arrange- 
ment. It is because of the highly stable characteristics with re- 
spect to half-frequency whirl that this system has been so in- 
tensively investigated. The main purpose of this paper, however, 
has been to present some of the steady-state characteristics of the 
pivoted pad bearing such as shoe load capacity, minimum film 
thickness, and optimum shoe pivot location. The authors hope 
at a later date to present information on such system dynamic 
characteristics as half-frequency whirl, shoe flutter, and rotor 
resonance frequencies. 

There are several operating characteristics that the authors wish 
to comment on at this time. One is the possible occurrence of 
shoe flutter with tilting-pad shoes. As Mr. Warner has stated, it 
seems desirable to construct the pads with low moments of inertia 
so as to keep the pad rolling and pitching flutter frequencies well 
above the operating speed of the unit. In the design of the au- 
thors’ unit, the self-acting shoes were designed with large inertia 


properties in order to obtain accurate film thickness data. The 
natural shoe flutter frequency in the direction of motion appears 
to be less than one-half of the design speed of the unit, i.e., 18,000 
rpm. No operating difficulties have ever been encountered due to 
shoe motion. 

Another factor which seems to have a strong influence on the 
dynamic behavior of a tilting-pad gas-lubricated bearing system 
is the design of the shoe pivots. 

Several different sets of pivots have been used in the testing 
and have resulted in a considerable variation in both shoe and 
rotor performance. One set of pivots employed had a cone and 
socket design which permitted an accurate alignment of the shoe 
with respect to the shaft for precise film thickness measurements. 
With this type of pivot, an unstable rotor orbit was monitored in 
which half-frequency whirl appeared to be present. This phe- 
nomena occurred with two different sets of shoes having different 
aspect ratios. In both cases, the apparent half-frequency rotor 
orbits occurred at a speed of approximately twice the rotor trans- 
latory resonance frequency. It was felt that sufficient friction 
present in this pivot design could have caused the shoes to behave 
more as fixed partial arc bearings rather than as pivoted pad 
bearings. When a second pivot design of a crown pivot on a flat 
bushing was employed, the unstable whirl orbiting was eliminated, 
but the observed shoe motion was three times the previous value. 
This might indicate that some external damping in a pivoted pad 
shoe is desirable to limit shoe motion but not so much as to cause 
the bearing to behave as a fixed pad and hence induce whirling. 
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